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ANNOTATION 


In  the  book  the  theory  of  a  turb¬ 
ulent  boundary  layer  of  a  compressible 
gas  Is  set  forth  or  based  on  a  study  of 
the  relative  change  of  the  coefficients 
of  friction  and  heat  transfer  with  growth 
of  the  M  number,  the  heat  transfer  factor 
tof/  and  the  wall  permeability  factor  b. 

The  existence  is  shown  of  a  limiting  law, 
corresponding  to  very  large  Re  numbers 
and  almost  full  of  self-simulation  of 
relative  changes  of  the  coefficients  of 
friction  and  heat  transfer.  On  this 
basis  simple  engineering  methods  are  pre¬ 
sented  for  solving  the  basic  problems  of 
friction  and  heat  transfer  in  turbulent 
flow  of  a  gas  past  a  solid  body. 

The  theoretical  conclusions  sre  com¬ 
pared  with  experimental  data. 

The  book  is  designed  for  scientists, 
engineers  of  aerodynamics  and  thermophysics, 
students  of  senior  courses  in  these  speci¬ 
alties  and  may  be  used  as  a  guide  for  prac¬ 
tical  calculations  in  design  offices. 
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PREFACE 

In  nature  turbulent  flows  are  the  most  wide-spread,  starting  with 
the  flow  of  water  in  a  stream  and  finishing  with  the  motion  of  inter¬ 
stellar  gas.  Flow  of  operating  media  in  machines  and  apparatuses  of 
contemporary  technology  in  the  overwhelming  majority  of  cases  is  turb¬ 
ulent.  During  high  speed  flight,  in  combustion  chambers  and  the 
nozzles  of  motors,  in  nuclear  reactors,  in  high-forced  gas  equipment 
and  so  forth  streamlining  of  solid  bodies  of  different  shapes  by  a 
turbulent  stream  of  compressible  gas  takes  place.  In  addition,  flow 
is  accompanied  by  Intense  heat  transfer  or  mass  transfer  of  some  sort 
of  matter. 

As  is  known,  aerodynamic  drag  and  heat  and  mass  transfer  are  de¬ 
termined  by  the  conditions  of  the  transfer  of  momentum,  heat,  and 
matter  in  the  boundary  layer  of  gas  or  liquid  forming  near  a  stream¬ 
lined  surface.  In  a  stream  of  compressible  gas,  all  phenomena  are 
complicated  by  innate  compressibility  and  by  changes  of  temperature, 
connected  with  high  glow  velocities.  In  this  connection,  the  great 
attention,  which  has  been  allotted  in  the  past  few  years  to  the  problem 
of  the  boundary  layer  of  a  compressible  gas  is  understandable.  How¬ 
ever,  if  the  theory  of  a  laminar  boundary  layer  can  be  considered  as 
basically  completed,  then  with  respect  to  the  turbulent  boundary  layer 
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the  position  to  now,  remains  far  from  satisfactory. 

The  semi-empirical  theories  of  near-wall  turbulence  of  Prandt.L- 
Karman  and  Taylor  to  a  certain  degree  permitted,  explanation  of  'hp 
existence  of  a  logarithmic  velocity  profile  in  an  isothermal  fluid 
current  with  weak  pressure  gradients  and  impenetrable  surface,  i’xtei  - 
sion  of  this  theory  to  noniso thermal  flow  with  zero  pressure  gradients 
was  carried  out  in  the  works  of  F.  I.  Frank*  and  V.  V.  Vcyshel’, 

A.  A.  Dorodnitsyn,  R.  Deisler,  L.  Ye.  ho  ikhmar.  7.  M  iyevlev,  £.  Von 

Driest,  W.  Dorrance  and  F.  Dore,  V.  P.  Mctulevich,  Yu.  V.  Lapin,  et  al. 
In  the  works  of  K.  K.  Fedyaevskiy,  V.  Stsablevskiy ,  L.  G.  Loytsyanskiy , 

et  al.  Certain  problems  are  considered  of  the  fl:  w  of  an  incompressible 

liquid  in  the  presence  of  a  pressure  gradient  in  the  flow  region,  not 
close  to  the  point  of  separation  of  the  boundary  layer.  In  addition, 
it  is  necessary  to  accept  a  series  of  additional,  physically  insuffi¬ 
ciently  valid,  and  sometimes  contradictory  assumptions.  As  a  result 
calculations  by  these  methods  substantially  differ  among  themselves, 
which  was  noted  repeatedly  in  the  literature.  Therefore,  it  is  not 
surprising,  that  in  engineering  practice  purely  empirical  methods  of 
calculation  of  friction  and  heat  transfer  in  the  turbulent  boundary 
layer  of  a  compressible  gas,  based  on  the  introduction  of  some  "gov¬ 
erning"  temperature,  have  obtained  wide  circulation.  It  is  necessary 
to  note  that  in  the  semi-empirical  methods  of  calculation  it  is  also 
necessary  to  introduce  a  "governing"  temperature  for  calculation  of 
physical  properties  of  gas  in  a  viscous  sublayer.  The  problem  of  the 
parameters  of  separation  of  the  turbulent  boundary  layer  remains  theo¬ 
retically  nearly  unexplored. 

This  monograph  is  devoted  to  discussion  of  the  limiting  law  of 
change  of  the  coefficient  of  friction  in  a  turbulent  boundary  layer 
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under  nonisothermal  conditions,  transverse  flow  of  matter,  and  pressure 
gradient  as  established  by  the  authors.  This  law,  derived  for  Re  -*•  oo, 
in  the  general  form  does  not  depend  on  empirical  constants  of  turb¬ 
ulence  and  is  not  connected  with  any  special  type  of  semi-emoiricnl 
theories • 

The  known  fact  of  the  weak  influence  of  the  Reynolds  number  on  the 
relative  change  of  the  coefficients  of  friction  and  heat  transfer  in 
connection  with  nonisothermal  conditions  and  transverse  flow  of  matter 
permits  extending  the  limiting  law  to  flows  with  finite  Re  numbers, 
with  good  accuracy.  As  a  result  we  managed  to  construct  relatively 
simple  methods  of  solving  the  integral  equations  for  the  momentum  and 
energy  of  a  turbulent  bounday  layer. 

Theoretical  results  are  compared  with  a  large  number  of  diverse 
experimental  data. 

It  is  possible  to  propose  that  the  possibilities  of  the  new  method 
are  by  no  means  exhausted  by  the  problems,  considered  in  this  mono¬ 
graph  and  will  be  still  developed,  supplementing  the  existing  theory 
of  a  turbulent  boundary  layer. 

It  is  assumed  that  the  reader  is  acquainted  with  the  princinles 
of  aerodynamics  and  the  theory  of  heat  and  mass  transfer. 

The  authors  will  be  very  grateful  for  all  critical  remarks, 
which  are  raised  by  readers  in  connection  with  the  material  considered 
below. 
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Basic  Symbols 

thermal  equivalent  of  work; 

relative  speed  of  blowing  (In  pumping)  of  matter  Into 
the  boundary  layer  through  the  wall; 

thermal  diffusivity  coefficient; 

speed  of  sound; 


A=-Jfca*£S —  wall  permeability  factor; 
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wall  permeability  factor,  referred  to  the  reel  value 
of  the  friction  coefficient; 
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local  friction  coefficient 
average  friction  coefficient; 

friction  coefficients  during  streamlining  of  a  flat, 
impenetrable  plate  to  an  unlimited  isothermal  flow; 

specific  heat  at  constant  pressure; 

specific  heat  at  constant  volume; 

diameter; 

diffusion  coefficient; 

surface  area; 
shape  parameter; 

mass  flow  rate; 

acceleration  of  gravity; 

shape  parameter,  the  ratio  of  the  depths  of  dlsplace- 
ment  to  loss  of  momentum; 

specific  enthalpy; 

pulsational  component  of  mass  velocity  (flow  +  rate); 

mass  velocity  of  transverse  flow  of  matter  through  a 
permeable  wall; 

criterion  of  phase  transition; 
index  of  the  adiabatic  curve; 
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/M  -  linear  dimension; 
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wj 
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f  “ 


s/=  —  - 
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length  of  the  path  of  hydrodynamic  mixing; 
length  of  the  path  of  thermal  mixing; 
full  length  of  the  body; 
exponent; 

Nusselt  number; 
exponent; 

Mach  —  Mayevskiy  number; 

Prandtl  number; 
pressure; 
heat  flow; 
heat  flux; 
gas  constant; 
radius; 

Reynolds  number; 

characteristic  Reynolds  number  of  boundary  layer 
temperature  recovery  coefficient; 
latent  heat  of  flow  transition; 

Stanton  number; 
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ri°K|  - 
rTKI  - 

r;  I  ku 

T^K)  - 

t  rci  - 


sr, 


absolute  temperature; 
lhhlbltlon  temperature; 
adiabatic  wall  temperature; 

calculated  temperature,  determined  by  formula  (2*45); 

temperature  on  the  centigrade  scale; 

ratio  speed  of  undisturbed  flow  to  the  maximum  pos¬ 
sible  speed  of  flow; 


X  - 


y  - 


the  same  for  the  local  value  of  velocity  component  w  ; 

X 

volume ; 

volume  flow  +  rate; 

pulsational  vector  component  of  velocity; 
specific  volume; 
dynamic  velocity; 
flow  velocity; 

velocity  outside  the  dynamic  boundary  layer; 

velocity  on  the  calculated  boundary  of  turbulent  nu¬ 
cleus  and  viscous  underlayer; 

longitudinal  gradient  of  velocity  outside  the  boundary 
layer; 

coordinate,  directed  along  the  flow  according  to  the 
outline  of  the  contour; 

coordinate,  normal  to  the  streamlined  surface; 
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KKil.l  I 
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r, 
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*  W  - 
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0*M  - 
A-M  - 

CM- 

f  - 

to  - 

*•  _ 

V--^- 


calculated  thickness  of  the  viscous  underlayer; 

heat  transfer  coefficient; 

angle; 

compressibility  factor  in  the  expression  for  turb  - 
lent  tangential  stresses; 

form  parameter; 

density; 

sign  of  difference; 
heat  transfer  factor; 
pulsational  component  of  density; 
thickness  of  the  dynamic  boundary  layer; 
thickness  of  the  thermal  boundary  layer; 
thickness  of  displacement; 
thickness  of  loss  of  momentum; 
thickness  of  loss  of  energy; 

coefficient  nonsimilarity  temperature  and  velocity 
fields ; 

coefficient  nonsimilarity  of  concentration  and  velocity 
fields; 

coefficient  of  aerodynamic  drag  during  flow  in  a  pipe; 

dimensionless  distance  from  wall,  expressed  in  the  form 
of  "local  Reynolds  number;" 
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Vt 
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dimensionless  thickness  of  the  viscous  underlayer; 
pulsational  component  of  temperature; 
dimensionless  temperature  (see  formula  [2.46]); 
form  parameter; 

form  parameter,  referred  to  c^q; 
coefficient  of  thermal  conductiuity; 
coefficient  of  dynamic  viscosity; 
coefficient  kinematic  viscosity; 
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density  of  the  medium; 
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relative  magnitude  of  tangential  stress  in  the  thick¬ 
ness  of  the  boundary  layer; 
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relative  change  of  coefficient  of  friction  at  Re**  = 

=  idem; 
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kinetic  temperature  factor; 

cross  sectional  area; 
dimensionless  velocity; 
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CHAPTER  I 

BASIC  EQUATIONS  OF  A  TURBULENT  BOUNDARY  LAYER 

1.1.  Equations  of  Motion  and  Thermal  Conductivity  of  a 
Plane  boundary  Layer  of  Gas 

List  of  Designations  Appearing  In  Cyrillic 
ct  =  wall  *  wall 

Kr  =  kg 

M  =  m 

A  stream  of  fluid  forms  near  a  streamlined  surface  a  dynamic 
boundary  layer,  i.e.,  a  region,  in  which  the  velocity  of  the  fluid 
changes  from  the  velocity  at  the  wall  (for  a  nonrarefied  gas  it  is 
equal  to  zero)  to  a  magnitude  very  close  to  the  speed  of  an  undistrubed 
flow  —  wQ.  In  the  presence  of  heat  transfer  and  diffusion  thermal  and 
diffusion  boundary  layers  appear.  In  a  thermal  boundary  layer  the 
temperature  in  practic  .  changes  from  Twall  to  TQ.  in  a  diffusion 
boundary  layer  the  concentration  at  the  diffusing  substance  changes 

n  ti 

from  Pwall  to  pQ  Strictly  speaking,  for  a  dynamic  boundary  layer 
the  boundary  conditions  take  the  form; 

jy  =  0,  wt=0;  y  -  cot  «,  =  ».!.  (1.1) 
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Due  to  a  sharp  change  of  speed  in  direct  proximitv  to  the  wall,  the 
exact  conditions  of  (1.1)  can  be  replaced  by  the  approximate; 

|y  =  0,  if,  =  0;  y  =  v,  =  .  I  ~  a\,|.  / *  „  \ 

\  <»*  •  <  ' 

where  e  is  a  prescribed  small  magnitude. 

In  that  sense  we  are  speaking  of  a  boundary  layer  of  finite  thickness 

5. 

In  experiments  the  value  of  e  can  coincide  with  the  sensitivity 

of  the  measuring  instrument. 

In  Pig.  1  is  given  a  diagram  of  a  boundary  layer  on  a  certain 

curvilinear  surface. 

In  a  plane  boundary  layer  in  the  absence  of  significant  transverse 

forces  (for  instance,  centrifugal)  the  following  conditions  are  ful¬ 
filled! 


*>  ' 


it 


< 


*Ii  +U 

dxt  if  ’  > 


(1-3) 


where  for  f.  are  understood  w  T,  and  p. 

X  X  | 


Fig.  1.  Diagram  of  a  boundary 
layer  on  a  curvilinear  surface. 
In  front  of  the  body  is  a 
shock  wave. 


In  connection  with  this, 
the  equations  of  thermal  conduc¬ 
tivity,  motion,  and  continuity 
for  e.  stationary,  plane  boundary 
layer  of  compressible  gas  on  an 
impenetrable  surface  have  the  form. 


(1.4) 
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t» t 
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J*£-  .£. 

H.t 


7r(^)  =  -iw’Jrr+w’J7r): 


dy  V  dy  I  '  '  " 

*(**»)'•',  _n 
d.t  *•  .  rfy  ~U‘  ‘ 


(1-5) 

(1.6) 


At  c  const 
P 


Si 

Sx 


=  c. 


St 

-3^—  ,  what  is  practically  true  for  a 


uniform  gas. 

To  these  equations  one  should  join  the  equation  of  state,  determ 
J.ning  the  magnitude  of  the  density  of  the  gas  p,  and  temperature  func 
Lions  of  the  coefficients  of  thermal  conductivity  and  viscosity. 

For  a  gas,  obeying  the  Clapeyron  —  Mendleyev  equation,  we  have 


*  =  -£- 


gRT 


(1.7) 


The  coefficients  of  thermal  conductivity  and  viscosity  are  re¬ 
lated  to  the  specific  heat  via  the  Prandtl  number. 


Pr=-ErL-- 


(1.8) 


for  undissociated  gases  the  values  of  Pr  and  c  change  little  with 

P 

temperature  and  pressure.  Therefore,  for  such  a  medium  practically 
one  may  assume  that 

-const.  (1.9) 

The  magnitude  of  Pr  depends  chiefly  on  the  number  of  atoms  in  the 
molecules  of  the  gases,  and  its  order  is  given  in  Table  1.1. 


Table  1.1.  Order  of  Pr  number  for  undissociated  gases. 


Atoms  1  2  3  *4 

Pr  number . .  0.66  0.75  0.84  1 


Outside  the  boundary  layer  frictional  forces  do  not  appear 

cHv  „ 

3^: —  =  0  )  and  for  steady-state  flow 
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-  df'  -  ft  7*1  dW» 

,tx  dx  "  (1.10) 
Substituting  in  equation  (1.4)  the  value  of  — from  equation  (1.5), 
we  obtain  the  thermal  conductivity  equation  for  a  plane  boundary  layer 
in  the  form  of  M.  P.  Shirokov  (at  c  *  const)  i 

It 


Here 


d 


T*  -r  (Pr-  I  t 

» 


=  [™x 


d  T * 
d.x 


+  Wy 


dT « 
dy  • 


(1.11) 


r*=.  t — 


Aw] 

2c  tp 


(1.12) 


is  the  temperature  of  stagnation. 


1.2.  Turbulent  Viscosity  and  Thermal  Conductivity 
in  a  Plane  Boundary  Layer 

The  real  (actual)  instantaneous  characteristics  of  turbulent  flow 
effect  at  every  point  disordered  oscillations  around  a  certain  mean 
value.  Thus,  for  the  flow  velocity  we  have 


tr  —  tr  ’■  V ,  (l  •  15  ) 

where  w  —  vector  of  the  averaged  velocity  at  a  given  point  of  flow, 

?  —  vector  of  pulsational  component,  giving  the  deviation  of 
the  true  velocity,  at  a  given  moment  of  time,  from  the  averaged  value. 

In  a  pompressible  gas  the  flow  velocity,  pressure,  temperature, 
density,  and  flow  +  rate  of  the  medium  pulsate.  The  corresponding 
equation  of  motion  of  a  plane  steady  state  boundary  layer  has  the  form 


it  x  0  y  ' 1  d  v 


-j,  K.)  = 


-  dw,  i  ~~  dwx 


(1.14) 
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Here  a  line  above  a  letter  signifies  averaging  during  a  time  fairly 
large  as  compared  to  the  period  of  pulsation. 

In  turn,  the  averaged  product  of  the  vector  components  of  pulsa¬ 
tion  of  the  flow  +  rate  J  and  vector  velocity  pulsation  Vx  are 

T^m  =  tVF,  «,TT*  ^  (1.15) 


where  6  is  the  density  pulsation. 

By  comparing  in  equation  (1. 14 )  the  members  and  jJVx.  ^ 

is  possible  to  arrive  at  the  conclusion  that  the  latter  value  can  be 
considered  as  some  tangential  stress,  appearing  in  the  averaged  stream 
under  the  influence  of  turbulent  pulsations. 

In  a  gas,  obeying  equation  (1.7). 


8=P - P = 


(1.16) 


and  correspondingly 


— +-£■ 

?  T 


(1.17) 


where  8  is  the  temperature  pulsation. 
The  value 


is  = 


s 

~*T 


(1.18) 


is  called  the  coefficient  of  turbulent  viscosity.  This  value  is  a 
complex  function  of  velocity  and  temperature. 

The  idea  of  the  coefficient  of  thermal  conductivity  is  developed 
in  an  analogous  manner. 

In  the  method  proposed  the  values  M-t  and  Xt  are  not  used  directly 


and  only  their  ratio  is  important 


1 


called  the  turbulent  Prandtl  number. 

This  value  near  a  solid  wall  is  close  to  unity. 


1.3.  Integral  Equations  of  a  Boundary  Layer 

The  equation  of  motion  can  be  written  recorded  ir  tie  form 


4Wn 
4  x 


+ 


d  "t 

*y 


-  v « 


X 


dtrv 

d.T 


4- 


d «  r 
d  v 


(1.20) 


true  for  laminar  and  turbulent  flows,  if  corresponding  velues  of 


variables  are  introduced  into  it. 


Integration  of  this  equation  with  respect  to  y  from  0  to  5  gives 

*  *’• -Jr  ‘ + ■»  *”■'»  "  = 


i 


On  the  outer  boundary  of  the  layer,  by  definition  t  -  0  and  wx  »  wQ. 
On  the  surface  of  a  streamlined  body 


Here 


"  ~  ”<■?«  -  0,  pK*y  =  /,. 

jy  =  Prr  O’ct 


(1.22) 


is  mass  flow  velocity  of  matter  through  a  surface.  This  flow  may  be 
as  a  result  of  a  change  of  the  aggregate  state  (evaporation,  conden¬ 
sation),  porosity  of  the  wall  (blast  of  gas  into  the  boundary  layer 
or  pumping  from  the  boundary  layer),  chemical  reactions. 

Subsequently,  a  surface,  on  which  /  0,  we  will  call  permeable. 

Taking  into  account  these  boundary  conditions  and  the  equation  of 
continuity,  it  is  possible  to  transform  equation  (1.21)  to  the  form 


+  _!!!!_  (2  4-//)  +  — -  a**  — — 

4x0%  H  *u 


~CT 


(1.23) 


where 


and  j»o 


4x 
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iHt  t  *  'v 


***.-■ 


In  this  equation,  called  the  Karman  momentum  equation 


•m ,  i 


is  the  thickness  of  displacement. 


(1.24) 


(1.25) 


is  the  thickness  of  loss  of  momentum. 

The  upper  limit  of  integration  y  -  co  corresponds  to  the  theory 
of  an  asymptotic  boundary  layer;  the  limit  y  ■  6  corresponds  to  the 
theory  of  a  layer  of  finite  thickness. 

Due  to  a  sharp  change  of  speed  in  the  interval  from  0  to  6, 

6*  and  6**  have  the  same  value  for  both  upper  limits,  i.e.,  there  are 
certain  "internal11  linear  scales,  the  one  and  the  same  within  the 
theory  of  an  asymptotic  layer,  and  in  the  theory  of  a  layer  of  finite 
thickness. 

The  analogous  integration  of  the  thermal  conductivity  equation 
leads  to  the  equation  of  the  energy  of  a  boundary  layer  (without 
consideration  of  radiation  and  internal  sources). 


rf.r  [  ir(  A  T 


fpt  j\  _  9cr 

<>  Po  “  g  w„  A  f 


(1.26) 


where  (AT)1  »(d/dx)  ( AT) ,  AT  is  the  temperature  difference, 

cpl  is  the  specific  heat  of  the  gas,  introduced 
through  a  permeable  surface, 

CpQ  is  the  specific  heat  of  the  basic  gas. 

The  quantity 


1  - 

hw< 


»)dy 


(1.27) 


is  called  the  thickness  of  loss  of  energy. 
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The  coefficient  of  friction  is 


c,  - 


(1.26) 


and  the  coefficient  of  heat  transfer  is 


«  =  Ssl- 

a  r 


The  dimensionless  forms  of  the  hee  coefficient 

and  Nusselt  numbers)  have  the  form. 


(1.29) 

'anefer  (Stanton 


St 


(  « 


,50) 


where  Iq  is  a  characteristic  linear  dimension. 

In  integral  (1,27)  $  is  a  cer~ain  dimensionless  temperature .  At 
Pr  ■  1  and  dp/dx  ■  0,  the  value  of  $  is  exactly  aqua,  to  the  ratio 

TWan-T*/TW4lii-To*,  where  T*  is  the  value  cf  the  inhibitian  temperature 

# 

at  a  given  point  and  is'  the  value  of  the  stagnat ion  temperature 

outside  the  boundary  layer. 

The  temperature  difference  is  determined  by  the  expression 


“  "  ,C7>  (1,52) 

where  Twa2.i  1®  adiabatic  temperature  of  the  wall,  i.e.,  that 
temperature,  which  is  established  on  a  given  surface,  if  the  latter 
ia  completely  heat-insulated  and  9wall  -  0, 

In  general  it  is  assumed  that 


7 J,=  r,+ 


where  r  ia  the  temperature  recovery  factor. 
In  a  turbulent  boundary  layer 


(1.35) 


l.e.,  for  gases  close  to  or  equal  to  un 


(1.54) 
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s 


l 


Fig.  2.  Diagram  of  a 
boundary  layer  on  an  axially 
symmetric  body. 


dx 


In  Fig.  2  is  given  a  diagram  of  a 
boundary  layer  on  an  axially  symmetric 
body.  Because  of  the  small  thickness 
of  the  boundary  layer  relative  to  the 
radius  of  curvature  R  ,  the  equations 

A 

of  motion  and  thermal  conductivity  do 
not  change  their  form  as  compared  to 
Planar  flow.  The  form  of  the  con¬ 
tinuity  equation  changes,  taking  the 

form: 

d(p  Wg  Rx)  _ q 

d  v  •"  . 

(1.36) 


In  connection  with  this  the  equations  of  momentum  and  energy  are 
somewhat  modified;  namely. 


Here 


tii>* 
- + 

rf  X 


(2  +  H)  f 
_ 4_  - _ 5 


(t 


4x 


<*l/l  —  fCT 

»>oto  «•#  tc^)fttWt^T 


(1.37) 


(1.38) 


* -^r  cosP)rf,':  (1-39) 

= f^O  --  co*  »)<* 

0 

*’"  =  J  n  - «)  (l  ±  cos  p)rfy.  (I*41) 
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1.4.  Characteristic  Reynolds  Numb e_r  • 
a  iBound'ary'  Layer 

The  trivial  determination  of  the  Reynolds  number  has  the  fo rmt 


where  wQ  is  characteristic  velccltyjt 

l  _  the  characteristic  linear  dimer,  si  on  of  ?  streamlined  body: 

Vq  the  kinematic  viscosity  at  a  characteristic  point  of  the 
system.  However,  in  boundary  layer  theory  sue:--  a  determination  of  e 
basic  hydrodynamic  criterion  is  ineffective.  tally,  from  the 
momentum  equation  it  is  clear  that  such  a  very  Import ®yt  integral 
characteristic  of  the  interaction  of  flow  with  6  body;  as  the  coeffi¬ 
cient  of  friction  cf,  is  connected  not  ..ith  the  length  but  with 
certain  "internal"  dimensions  of  the  boundary  layer  —  6#*  and  5*. 
Besides  the  structure  of  the  momentum  equation  shows  that  the  first 
of  these  quantities  is  the  most  important. 

The  universality  of  the  linear  characteristics  6*,  &**  and  b  ** 
becomes  especially  evident,  if  we  remember  that  ti  sir  determination 
is  not  connected  with  the  idea  of  a  layer  of  finite  thickness  and  the 
change  of  the  upper  limit  in  integrals  (1.24),  (1.25),  and  (1.27) 
from  6  to  oo  does  net  change  their  value.  In  connection  with  this 
characteristic  Reynolds  number  for  a  dynamic  boundary  layer  it  is 
expedient  to  construct  in  the  form  of  the  quantity: 


or 


JN**  = 


(1.4?) 


*4. 


(1.44) 


In  an  isothermal  turbulent  boundary  layer  far  from  the  point  of 
breakaway,  the  ratio  5*/6**  is  almost  constant  (see  formula  (?.!?)), 
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i.e.,  in  this  case  the  numbers  Re4*  and  Re*  differ  only  by  a  constant 
factor. 


For  a  thermal  boundary  layer  we  have  correspondingly 


Re, 


'*0 


(1.^5) 


The  dominating  role  of  the  quantity  6**  in  the  momentum  equation 
and  consideration  of  the  expedient  uniformity  of  the  characteristics 
of  dynamic  and  thermal  boundary  layers  force  us  to  select  as  basic 
modifications  of  the  Reynolds  the  expressions  (1.43)  and  (1.45). 
Regarding  selection  of  the  quantity  it  is  more  convenient  to 
relate  to  conditions  outside  the  boundary  layer.  Such  a  determination 
is  convenient  from  a  computation  point  of  view,  inasmuch  as  the  param¬ 
eters  of  undisturbed  flow  usually  are  known. 

The  relation  between  the  Reynolds  number 


*  '  (1.46) 

and  the  numbers  Re**,  Re^**  is  established  via  the  equations  of  momentum 
and  energy,  if  the  following  relations  are  known* 


and 


'/-/iHo 


(1.47) 


(1.48) 


These  relations  are  called,  respectively,  the  law  of  resistance 
and  law  of  heat  transfer. 

Definition  of  these  laws  is  a  basic  problem  of  boundary  layer 
theory. 


1.5.  Similarity  of  the  Velocity  and  Temperature  Fields 
The  equations  of  motion  (1.5)  and  thermal  conductivity  (1.11) 
become  identical  relative  to  the  quantities  w  and  T*  upon  fulfillment 

Jv 

of  the  conditions* 


[  —  —  ■  —  0;  Pr=\\. 


(1.43) 


But  identity  of  the  differential  equations  signifies  identity  of 
their  integrals  during  similarity  of  t  ‘  -  V  ••  •  -  ' 

Consequently,  during  preesureless  st.r-  nr'1 4 nir.g  f  a  solid  body 
by  a  polyatomic  gas  (Pr  »  1}  and  fo>*  *  •  ’.arly  fined  boundary  con¬ 

ditions  the  longitudinal  flow  velocity  and  stagnation  temperature 
fields  are  similar,  i.e.. 


Ter  -  T'  ■ 

r"_l«  '  (1.50) 

These  conditions  are  exactly  fuj.iix.Lea  during  'Ogi c\.o Inal  flow  by  an 
unbounded  stream  of  a  polyatomic,  ..ianar  gas  '  ...  a  constant  surface 
temperature . 

Similarity  the  velocity  and  tempera  ure  fields  at  Pr  =  1  means 
that  between  the  heat  flu.  q  and  the  Tangential  stresses  t  tnere 
exists  the  relationship 


(1.51) 


or 


St  = 


*L 
2  ' 


(1.52) 
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CHAPTER  II 

LAWS  OF  RESISTANCE  AND  HEAT  TRANSFER 

Definitions  of  Cyrillic  Items  in  Order  of  Appearance 

t  =  t  =  thermal,  turbulent 
ct  =  wall  =  wall 
Kp  =  cr  =  critical 

2.1.  Tangential  Stresses  in  a  Planar  Boundary  Layer 

of  a  Compressible  Gas 

Molecular  friction  is  commensurate  with  turbulent  only  in  the 
thin  near-wall  layer.  This  region  of  the  turbulent  boundary  layer  is 
called  the  viscous  underlayer. 

The  relative  thickness  of  the  viscous  underlayer 


in  a  flow  of  an  incompressible  liquid  is  1-3$,  and  in  a  flow  of  a 
compressible  gas  can  reach  10-15$. 

In  the  remaining  part  of  the  turbulent  boundary  layer,  the  so- 
called  nucleus, 

and 

3  w9  P3?  ♦  vM  v» * 

77V,? 


(2.2) 

(2.3) 


where 


—  a  coefficient,  talcing  into  account  the  influence  cf  density  pulsa¬ 
tions  on  momentum  transfer. 

Below  the  sign  of  the  average  values  of  p,  w,  and  1  will  be  Alltel 
As  can  be  seen  from  formula  (7.2),  *1  fur  :  r.mc  o 4.  ?  *  pr.'hiern  f  the 
theory  of  turbulence  is  the  connection  between  the  peerage  produo4"  of 
the  vector  components  of  the  pulsation  .1  component  of  the  flow  ve  locity 

and  the  average  velocity.  This  relation  is  expressed  very  generally 
by  the  Prandtl  formula. 


\Tv  =  it  —9  \* 
v‘v'  V  *>  )’ 


.4) 


where  l  is  a  certain  proportionality  factor,  having  the  dimension  of 
length.  Usually  the  quantity  l  is  called  the  length  cf  the  mixing  path 
Formula  (2 .h  )  follows  from  dimensional  analysis,  if  we  assume  that 
in  che  region  of  significant  changes  of  the  average  velocity,  momentum 
is  proportional  to  the  derivative 

A 

Introducing  the  value  V  ~  from  {  ’.41  in  (2.2)  we  can  write  that, 

y  \  /  \  '  7 

a  turbulent  nucleus 


-  = t. 


tv 


(2.5) 


ihis  expression  leads  to  the  form 


'!  _  P)  lt  y 

1  •  H*'  1  "• 


(2.6) 


;:neri 


wall 


is  the  distributive  law  of  tangential  stresses  along 


:he  thickness  of  the  boundary  layer. 


2.2.  Law  of  Resistance 
From  equation  (2.6)  it  follows  that 


(2.7) 


»-i  *  t 


m  * 


i 


where  y.  is  the  thickness  of  the  viscous  underlayer,  i.e . ,  the 

coordinate  of  the  lower  boundary  of  the  turbulent 
nucleus  of  flow. 

W1 

0).  =  — -  is  the  dimensionless  velocity  at  this  boundary. 

1  wu 

Equation  (2.7)  is  written  in  a  form,  not  connected  with  any  sort 
of  theory  concerning  a  boundary  layer  of  finite  thickness.  The  transi¬ 
tion  of  the  terms  of  the  theory  of  a  boundary  layer  of  finite  thickness 
is  accomplished  according  to  equation 


where 


,r _ _ il _ 

I  //ITJ  J  ir tzj  ’ 

y\  L  * 

6  is  the  thickness  of  the  dynamic  boundary  layer; 
i  is  the  dimensionless  distance  from  the  wall; 


(2.8) 


y1 

^  =  y  is  the  dimensionless  thickness  of  the  viscous  underlayer. 
The  upper  integral  of  equation  (2.7)  has  the  same  limits  in  the 
theory  of  an  asymptbtic  boundary  layer  and  in  the  theory  of  a  layer  of 
finite  thickness. 

Let  us  introduce  the  relationship 


\  c/o  !*f 


(2.9) 


where  cf  is  the  value  of  the  local  friction  coefficient  at  given  con¬ 
ditions  of  streamlining  of  a  body; 

cf0  is  the  value  of  the  friction  coefficient  at  an  isothermal, 

gradient  less  streamlining  of  an  impenetrable  wall,  i.e., 
during  isothermal  streamlining  of  a  planar  impenetrable 
plate  by  an  unbounded  flow. 

Comparison  the  friction  coefficients  is  carried  out  for  identical 


values  of  the  characteristic  Reynolds  number  of  the  boundary  layer, 
because  Integration  in  (2.8)  is  carried  out  with  respect  to  the  thick¬ 
ness  of  the  latter. 

Multiplying  both  parts  of  equation  (2.6)  by  the  quantity  ~0, 
it  is  possible  to  write: 


i 


(2.10) 


-1 


In  the  formulas  T„  Is  the  distributive  law  of  one  ouantity  — - 

0  '  \-3.n 

along  the  thiciu.efas  of  the  isctheimal  boundary  lay^i  at  an  impenetrabit 


plate,  with  streamlining  of  an  unbounded  un.'  incite  5  fl  .  (i.e.  for 


dr 

ax 


-  0). 


From  the  given  formulas  it  is  clear  that  for  establishment  of  law 
of  resistance  it  is  necessary  to  know  the  iaws,  determining  the  quanti¬ 
ties  l,  y, ,  t,  ■£— ,  and  3.  The  quantity  w.  ce;  t-  computed  cy  oalcula- 
1  p0  1 

tion  of  molecular  friction  and  molecular  thermal  conduction,  i+'  the 
q  tar.tity  y„  and  the  relation  u(T}  r, re  known. 

U- 

In  the  contemporary  semi-empirical  theories  equation  (2.7)  is 
ca_culated  on  the  assumption  that  p  =  0,  i  =  l  ~  y,  or  v/hich  is 

practically  the  very  same, 


a  n ,  :  g* 

/  Or 


)■ 


nd  the  dimensionless  thickness  . f  the  viscous  underlayer. 


maintains  the  value,  found  experimentally  for  Isothermal  flow  on  a 
plate,  if  the  viscosity  is  related  tc  a  certain  "determining"  tempera¬ 
ture,  It  is  clear  that  boundary  layer  theories  for  a  compressible 
gas  based  on  such  assumptions  cannot  lead  to  sufficiently  reliable 
results. 

However,  equation  (2,10)  obtains  special  properties  at  Re  -*■  oo, 
to  which  attention  has  not  been  paid  until  now.  We  will  now  turn  to 
an  examination  of  these  properties. 


2.3,  Values  of  the  Quant It le  3  P , 
and  at  Re  — -  co 

We  will  show  that  at  Re  -•  oo  the  quantities  p,  to,  and  approach 


".'K*rz-,‘ri  ti 


zero. 


Equation  (2.4)  assumes  the  existence  of  a  correlation  between 
the  vector  components  of  the  pulsotional  portion  of  the  velocity  of 


the  form: 


Analogously,  for  temperature  pulsation  it  is  possible  to  assume  the 


relationship: 


From  (2.2)  it  follows  that 


y  _  .  JJL 

'  *y  *y 


=Vy*  v>  V  -j-;’ 


(2.12) 


where  v 


*  =  j/" is  the  dynamic  velocity. 


With  accuracy  up  to  the  coefficients,  taking  into  account  the  in¬ 
fluence  of  density  pulsation,  the  relationship  between  the  heat  flux 
and  tangential  stress  in  a  planar  turbulent  boundary  layer  has  the  form 

— ? —  ■—  g  c  -gy..- 

*  TT-t,  • 


In  the  degree  of  similarity  temperature  and  velocity  fields  con 


dition  (1.51)  is  fulfilled,  i.e., 

v^e  _  a  7- 
VTv, 


(2.13) 


Introducing  these  relationships  in  formula  (2.5),  we  can  write 


that  in  order 


\  *o  Wo  f  f>  /  T 


(2.14) 


The  transverse  velocity  component  at  an  impenetrable  wall  in 
succession  is  equal  to  the  quantity: 


i.e,,  in  this  case 


41  i  <n**  &  tf 

gjc  tlx  0  ••• 


+  r  (2.15; 

Let  us  consider  flow  with  "disappearing  viscosity,"  i.e.,  with 
U  -*■  0.  In  this  case  Re  -►  co,  and  the  friction  coefficient  cf  -►  0. 
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Correspondingly,  the  coefficient  (3  also  approaches  zero.  At  a  permeable 

wall  the  maximum  additional  component  of  the  velocity  r.  opponent.  w. 

V 

occurs  at  y  =  0  and,  as  is  shown  in  Chapter  IV,  is  eoual  to  w.  ,,  £ 

co  a  a  .t.  x 

g  — £9-2  b  w„,  where  b  is  a  finite  o.'urMty,  Cc-se ^ently.-  at 

r*  f  1  O  Y* 


cr  O' 

pwall 


cr 


ro. 


cf  -►  0,  the  quantity  Wwg also  approaches  ze- 

The  equation  of  motion  of  a  planar  boundary  layer  in  direct 
proximity  to  a  wall  has  the  form  (since  in  this  case  v;v  ~  0). 


*SL 

dx 


d  x 

<>y 


P“'j 


d 

d  y 


(2.16) 


3w. 


At  y  <  y±  <X  6  t  =  y,  — ~  and,  in  the  first  approximation,  p  =  Pwali 

by 

M  =  Pwalli  vr  =  wwall.  Here  for  the  quantity  w  ,r  ,  is  understood  the 


waxl 


value  of  w  at  y  =  0. 

Integrating  (c.l6)  we  have 


'  - J  -tr  d> = 1  (-£• + «-  tt  j dy 


* 


y-i-Ptrtt'cr*,- 


'  O  17) 
V  *-  •  -  I  1 


Correspondingly,  the  velocity  distribution  '  n.  the  nearest  proximity 
of  the  wall  is  determined  by  the  equation: 

(2.18) 


*CT  ~~  =  +  ““  y  +  Per  0*4  - 


dj  *  dx 

By  reducing  this  equation  to  the  dimensionless  form  and  taking 


dm 

dt 


—  ^Ict 


.10), 

we 

have 

*  l- 

*/ 

i 

V 

2 

_L_  ft  a.  .gcTfcr.  \ 
**’  7  Po»o  )’ 


(2.19) 


whe  re 


’CT 


the  characteristic  Reynolds  number  of  the  boundary  layer,  referred  tr. 
the  wall  temperature. 

The  quantity 


/  = 


dm,, 


dx 


(2.20) 


-18- 


can  be  considered  as  a  measure  "of  aerodynamic  curvature"  of  a  stream¬ 
lined  surface  and  is  called  the  form  parameter. 

During  streamlining  of  an  impenetrable  surface  wwall  =  0  and 
integration  of  equation  (2.19)  gives: 

J/t 


*# 

+  _  , 

_i_/ 

►€T 

»**  l 

-LLj.) 

21**  /' 


(2.21) 


6** 

The  value  of  is  finite  for  any  Re  numbers.  The  value  of  cjd^ 
by  definition  lies  between  0  and  1.  The  coefficient  of  friction 
is  always  inversely  proportional  to  the  Re  number  in  a  degree,  less 
than  1  (see,  for  instance  formula  (3.24)). 

Taking  into  account  these  circumstances  and  reducing  equation 
(2.21)  to  the  form: 


2»**«, 
t  M  * 

'  T 


(2.22) 


we  note  that  at 

Re  —  0. 

Thus,  the  thickness  of  the  viscous  underlayer  decreases  with 
growth  of  the  Re  number  relatively  faster  than  the  thickness  of  the 
whole  turbulent  boundary  layer. 


At  f  =  0 

at  f  =  f  v  0, 
cr 

cf  =  0  and 


Hl  = 


Re"  — 

f,** 


-  +±2i 


’I  »l* 


’CT 


a  2 

w,  — 


(2.23) 


i.e.,  at  the  point  of  breakaway  of  the  boundary  layer. 


Re»p  i  — 


•in 

V«T 


(2.24) 


Since  as  the  critical  Reynolds  number  of  the  viscous  underlayer  i 
always  finite  then  at  Re  -*■  go,  oo^  0. 

An  analogous  result  is  also  obtained  for  a  permeable  plate,  i.e., 
in  general: 

(2.25) 


2.4.  Limiting  Law  cf  Resistance 

At  Re  -*•  co,  -*•0,  -*  0,  p  -*  0,  and  accoru.ng  x  ,  ormu.ia  [J.'lj.  ) 

z-YTiv-r 


f  2 .  ^ 


Let  us  expand  the  f  un.c+i^r.  7  i  ^ r'"  1 


u  degree: 

the  perturbation  factor  aria  let  nr  es  »  ate  th  sum  of  the  terms  i,rorr 

1  =  2  to  1  =  cc  by  A$.  Then  in  accordance  with  form  (2.11) 


z  =.-  z„  f  hi*  i 


;  ? . 27 ) 


*  * 

The  quantity  Zq  =  1  -  uuq,  which  follows  ‘  rom  (\  ....  ),  if  In  this 

equation  are  inserted  p  =  p^,  x  =  x^  and  H  -  1 . 

Thus,  it  is  possible  to  write  that 


Z  -  1  —  m|(i  -*  1<S>  8 


I  T 


At  Re  -*  oo  aj^  -*•  0,  cfp  — >  0.  lons'-'Ciuently,  if  the  function  Av  at 
Re  -*■  cd  remains  finite  or  approaxhes  i n x -•  riity  more  weekly  than  V 
approaches  zero,  then  there  occurs  the  conuition: 


Re  •  » 


-  ? 


(2.29) 


Thus,  during  certain  conditions  there  exist  certain  limiting  laws 
of  the  relative  influence  of  nonisothermalness,  compressibility,  and 
other  disturbing  factors  on  the  coefficient  of  friction  in  a  turbulent 


boundary  layer,  determined  by  an  integral  of  the  form: 

rf 1,1  ^ 


f 


.» I 


/  n'  — £l 
’■ 


(2.30) 


More  detailed  information  about  the  properties  of  the  limiting 

relative  laws  of  friction  in  a  turbulent  boundary  layer  can  be  obtained, 

if  one  considers  in  a  sufficiently  general  form  the  connection  between 

the  length  of  the  path  of  mixing  l  and  the  coordinate  y.  Let  us  expand 

the  function  -L(£)  in  a  series  by  degrees  of  the  coordinate  £: 
o 


-s>n- 


(2.51) 


Experimental  material  shows  that  the  quantity  w  can  be  considered 

i+l 

as  some  universal  constant,  but  the  sum  of  the  terms  is  always 

finite,  although  the  coefficients  are  also  functions  of  the  non- 
isothermalness,  pressure  gradient,  transverse  flow  of  the  substance, 
and  other  disturbing  factors. 

In  particular,  during  diffusion  flow  of  an  incompressible  liquid 
on  an  Impenetrable  surface 

0<2S*'“  ,<I- 


From  experimental  data  it  follows  (see  chap.  Ill)  that  at  Re  — ►  oo 


Introducing  these  relationships  in  (2.26)  and  expanding  the  func¬ 
tion  VT  in  a  series  by  degrees  of  £>,  we  find  that  at  Re  -+  oo 


z--^ 


—  -|n  ii 

In*#** 


Hi  *#•* 

From  the  data,  given  in  chapter  III,  it  follows  that  on  an 

d] 


(2.52) 


impenetrable  surface  at  =  0  the  relative  thickness  of  the  viscous 
underlayer  is  related  to  the  Reynolds  number  by  the  relationship: 


.=  -3! - 

‘  |  He**  * 


c  ^  ^ 

In  addition  the  values  ik  ,  - ,  ¥  always  are  finite. 

±  5 

Putting  this  value  at  in  formula  (2.32),  we  find  that  in  the 
considered  case  at  Re  oo: 


Z-l- 


r  - 

In  In  **”  -  In  — 1 — j  1 
_ il _ 

In**” 


1. 


In  Chapter  V  it  is  shown  that  during  diffusion  flow  (-^  >0)  of  an 
isothermal  boundary  layer  on  an  impenetrable  surface  at  the  point  of 


breakaway  of  the  boundary  layer: 


roust 


1  Rf** {it> 

Putting  the  value  ^  in  (2.32),  we  find  that  at  **  —  m  7  — ►  0.65. 

Thus,  in  the  general  care  the  integral  v ^ ;  . ..  qua!  l  c  oertal 
value  Z  ,  not  being  a  function  ct  F  -yc'  1  in  number.  V"  -  c-s  the  pa^am- 

eter  is  exactly  equal  to  unit  or  cl  se  to  it. 


2.5.  Approximation  of  r.be  Jiangc-ntta ^ 

Stresses  rroi fro* 

From  determination  of  a  dynamic  boundary  layer  1  finite  thickness 

we  have  the  condition: 


'6  =  0.  *r  -  tfT;  | 

.6”  ».  -~e. 


) 


From  the  condition  of  smoothness  of  the  *  unction  1  (£)  at  the  point  £  =■• 

*  1  it  follows  that 

1  •  0 

i  *«”!;  «  (2.34) 


In  the  region  ^  *  0  with  accuracy  up  to  small  quantities  of  the  second 
order  equation  (2,17)  should  be  satisfied. 

Conditions  (2.17),  ^2 .33;,  and  (2.34;  art  outisi i-u  by  th-  cubic 
parabola, 

t  =  1-3S*-f  2V-h[-\l-r  ~l)\  (2.35) 


where 


?Pcr  P'ct 
tf  hi  «•<* 


(2.36; 


(2.37) 


The  first  of  these  quantities  is  a  certain  modification  of  the 
form  parameter.  The  second  quantity  characterizes  the  influence  of  the 
feeding  or  removal  of  a  substance  through  the  surface  of  a  streamline a 
body.  We  will  call  it  the  parameter  of  wall  permeablity. 

Subsequently  we  will  deal  again  with  several  modifications  of  the 
form  parameter  and  permeability  parameter. 


At  gradiant  less  streamlining  of  an  impenetrable  wall  according 


to  the  given  approximation: 

;  =  1  —  3t*  +  2;* . 


(2.58) 


Correspondingly 

—I—=  I -r  (Ac  +  £,••)/($).  (2.39) 

where 

A 

The  value  of  this  function  changes  from  1  at  £  =  0  to  at  £  =  1. 
To  formula  (2.39)  corresponds  the  relation 

=  *+  (A.S-1 (2.40) 

H 


where 


A#  = 


2> 


/; 


for*«T 


(2.41) 

(2.42) 


2.6 .  Approximation  of  the  Temperature  Profiles 

Inasmuch  as  pressure  across  boundary  layer  does  not  change,  then 
in  accordance  with  (1.7)  the  density  in  equation  (2.2)  at  a  uniform 
boundary  layer  can  be  expressed  by  the  density  of  an  undisturbed  flow 
by  formula 

»=*-*-•  (2.45) 

Thus,  for  solution  of  equation  (2.30)  it  is  necessary  to  know  the 
connection  between  the  temperature  and  velocity  fields.  This  relation 
was  known  previously  only  for  the  case  of  similarity  of  the  velocity 
and  inhibition  temperature  fields  examined  in  section  1.5. 

For  gases  the  Pr  numbers  are  equal  to  or  differ  little  from  unity, 
i.e.,  one  of  the  basic  conditions  of  the  existence  of  such  similarity 
is  always  fulfilled.  Therefore,  it  is  possible  to  put  dependency 


(1.50)  in  the  base  of  the  unknown  relation.  However,  lx 
to  give  it  a  form,  taking  Into  account  the  disturbance  t 

due  to  the  independent  effect  of  the  inMbitl  .0  •>*  ha  oy 

f  1,  which  is  seen  from  the  structure  of  equation  L. .  1 1 ) . 
Let  us  take  the  re  tat  >  .hi.- 


ic  necessary 

similarity, 

<— 

nv; 

nt  Pr  y 

o 


frr  T 
Tn  - 


(2.44) 


'  1. 


T*  —  T  {  r  (:)  - 


(2.45) 


The  form  of  the  functions  e(£)  and  r(y),  in  general,  depends  on  the 
pressure  gradient,  non! sothe rmai ne s 3 ,  and  mas.  o answer. 

At  e  =  r  =  1,  T  =  T  and  equation  n,  transforms  to  (1.50), 
i.e.,  there  occurs  exact  similarity  :.f  toe  velocity  and  stagnation 
teirpe rat ure  fields. 

Let  us  require,  that  on  the  1  undai  i.us  ji  trie  thermal  layer  the 


u entity 


70  ■  J> 

r<  C 


(2.46) 


satisfied  the  same  conditions  as  the  quantity  co  on  the  boundaries  of 
a  dynamic  layer.  Then  at  y  =  6^3-  =  1,  to  which  correspond  the  values 
=  T  and  r(£)  =  r.  In  a  turbulent  boundary  layer  the  quantity  r 
has  the  order  i.e,,  for  gases  close  to  unity.  Therefore,  with¬ 

out  apprecialbe  error  it  is  possible  to  assume  in  all  sections  of  a 
thermal  boundary  layer  r(£)  «  r  and 

r*r':''UFT'  (2-47) 

The  problem  of  the  function  e(£)  is  considered  in  section  2.7. 

The  value  of  the  dynamic  temperature  can  be  connected  with  the 
dimensionless  velocity 


w. 


(2.48) 


by  the  relationship 


where 


2tf,T„ 


iff**.* 


(2.49) 

(2.50) 


is  the  Mach  —  Mayevskiy  number,  referred  to  parameters  outside  the 

c 

boundary  layer;  k  =  —E-  is  the  index  of  Poisson's  adiabatic  line, 

cv 

From  (2.44)  and  (2.47)  it  follows  that 


or 

=  1>  (*•  !)<•»-, 

»• 

(2.51) 

— J—  =  ^- —  ($•  —  1)  t«s. 

(2.52) 

The  quantity 

(2.5?) 

is 

called  the  temperature  factor; 

The  quantity 

=  ~r-  »  1  -f-  r  1  M- 

M  * 

(2.54) 

is 

the  kinetic  temperature  factor,  determining  the  degree 

of  aero- 

dynamic  heating  of  a  body; 

(2.55) 

is  the  heat-transfer  factor,  since  at  =  0  adiabatic  streamlining 
of  a  body  occurs; 

At  Aip  >  0  a  body  gives  up  heat  to  a  flow,  at  AiK  0  a  body  receives 
heat  from  a  flow. 


2.7.  Coefficient  of  Nonslmllaritv  of  the 
Temperature  and  Velocity  fields 

The  velocity  and  temperature  fields  in  the  nucleus  of  a  turbulent 
boundary  layer  are  well  approximated  by  exponential  formulas  of  the 


form: 


In  flow  regions,  not  close  to  the  point  of  breakaway  of  the 

boundary  layer  (see  chaps.  IV  and  V),  both  exponeiv  „  i . ,  1.  t> 

the  profiles  are  very  complete. 

From  (2.56)  It  f ol low 3  that 


where 


(  .00 

( a.  58) 


Far  from  the  point  of  breakaway  n  n,  ,  1.  ,e„.  trio  re  exists  ? 
relative  similarity  of  the  temperature  ar- d  ve_  tv  fields,  expressed 
by  tne  formula 

0«eo>.  (1.59) 


'T'he  ratio  in  a  uniform  hour. lore  layer  i  -'ete mined  by  formula 

Po 

•  ;,v)  and  equation  (2.51 ) .  In  addition,  two  cases  ar.  s.otinguisbed : 
when  the  thickness  of  the  dynamic  boundary  layer  is  less  than  the 


’hi  ' - s  of 

the  thermal  and 

who n  the r  .  is  ar.  invei 

Tn  the 

first  case 

6  < 

ar.  d  l i •  e  roc e  s  s  0 f  h 

iv-  * h  whole 

thi.ck.nes" 

•P  +  0 

d"T'<anr!  c  layer.  Cor>v*e 

r>n  tri  OT1 

C  <  y  <  >: 

fcrmul 

a  '2 . 52 )  ' s  valid , 

in  the  second  case  6  >  5^,  i.e.,  in  the  region  &t  <  y  <  f,  heat 
transfer  .1  s  absent  and  the  temperature  is  determined  by  the  condition 
TT  =  ^onst.  Correspondingly,  in  the  region  0  <  y  <  5^  formula  (2.52) 
Is  valid,  and  in  the  region  5^  <  y  <  6  —  formula 


(2.60) 


2.8.  Limiting  Law  of  Heat  Transfer 
In  the  nucleus  of  a  planar  turbulent  boundary  layer  the  heat  ilux 
along  the  normal  to  the  streamlined  surface  is  determined  by  formula 

f ~f.=  -  c,rv^(l  -3t)r  (2.61) 

where  is  a  coefficient,  taking  Into  account  the  influence  of  density 
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pulsations  on  turbulent  heat  transfer. 
Analogously  to  (2.4)  we  can  write  that 

Vp>  =  t±ZL  .JZ 

'  i  y  dy 


(2.62) 


In  the  general  case  /  i.  From  (2.57)  -c  follows  that 


dm 


•  *  —  *? 

*  "r  ft  *»  . 


(2.63) 


At  c  =  const  equation  (2.61)  Is  reduced  to  the  form 

St  —  *  ?  >  1  -  ftr  ji  d*  i# 


tp  T«  »o 


djr  dy  * 


(2.64) 


where  q  =  ■?- 


^wall 

In  a  region,  not  close  to  the  point  of  breakaway  n  «  n^,  and  from 
(2.64)  it  follows  that: 


where 


H, 


%•  \ 

‘  fit,;-’ 


(2.65) 


(2.66) 


z,=  V  st. 


JVt? 

*IT 


dy 

Is 


(2.67) 


Transition  to  a  thermal  boundary  layer  of  finite  thickness  is 
carried  out  according  to  equation 


(2.68) 


where 

t  6t 


In  these  equations  StQ  is  the  value  of  the  Stanton  number  for  the 
conditions  P  “  P0'  lx  =  °'  wwall  =  0. 

The  quantities  Bt0,  and  possess  the  properties  of  the 
quantities  co^,  c^0,  p  and  Z,  i.e.,  at  Re  -*■  oo,  St0  -+  0,  0,  p^  -►  0, 
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Correspondingly,  at  Re  -*■  so  equation  (2.6b)  mu.  •  ^  -.miming  soiutioi 
of  the  form: 

*<  ( i  V  f 

'  V  '  1  -? 

This  equation  is  ana]  ego  us  :  equu  :■  . ;  r )  T‘or  vm  i  smiting  law  o.f 

resistance . 


2.9.  Approximation  of  tm  U  t  yrofil, 

From  determination  of  the  thermal  l  un fo.ry  _.„yer  of  finite  thici- 

we  have  the  conditions : 


»» 


-0, 


! 


( o  7  0  '1 
\  —  •  1  / 


she  condition  of  smoothness  of  Ur.  -.unction  r(i*tU  at  the  point 
i  it  follows  that : 


ii-j  =(!. 

if  ^  1 


(2.71 


The  thermal  conductivity  equation  near  the  wall  (v;  -  0,  T  -  T+) 

X 


v’^  form: 


iq  d  Tr 

=  «>Te»Vrt - 


f o  70' 

dy  . .  <*v 

Integrating  (2.72),  we  find  that  in  the  neighborhood  of  the  wall 

9-  Vct  h  0-Tcttti,  f*.  (2.T7) 


These  conditions  are  satisfied  by  the  cubic  parabola: 

J-l-  35;  +  25?-t*ltO(l  ;T)«. 


Correspondingly 


where 


J#=l  -  35t  +  2tf  ; 


Per  *ct  . 
5»poVo  ’ 


(2.74) 


(2.75'' 

(2.76) 

(2.77) 
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vr-tA-Miw "  -a 


b  =  —  *rT<>CT  . 
’  * 


/O.)  -•  (»,  + 1) 


-i 


(2.78) 

(2.79) 


Cp  /  const  the  quantities  b^t  and  b^  should  be  multiplied  by 


c 

the  ratio  -P  a11. 

Cp0 

At  Pr  =  1  and-|5-  =  0  8  =  8,,  8  =  1,  and 
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CHAPTER  HI 

LONGITUDINAL  STREAMLINING  OF  .  HI  HIT"  FABLE  PI  ATE 

Definitions  of  Cyrillic  Items  in  Order  ’ ;  pea ranee 
ct  =  vail  =  wall 
Kp  =  cr  =  cr’.t-  'a) 
t  =  t  =  thermal ,  .  rsnlent 

1.1.  iso  Thermal  roan. :ary  Lay  er 
trie  chapter  the  stre&m] j n:’ ng  of  flat  plate-  by  an  unbounded 
.  of  gas  for  tne  conditions 

•  zz  'f  and  7t  const 

dx 

'  d  e 1  ‘o  d  . 

•  ooundary  layer  for  the  variables  and  T  is  considered  in 

:  r  er  V. 

A  large  quantity  of  experiments  confirms  for  the  region  £  <  0.4, 
m-  logarithmic  distribution  of  velocities  in  a  well-developed  tur- 
. t.  isothermal  boundary  layer  is  expressed  by  the  Prandtl-Kikuradze 
ula: 

<p=*5,5+2,5  Inrj.  (5«l) 

Fere 


The  last  quantity  is  called  the  dynamic  velocity. 
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,  r  > 


14  i  ’ 


'  » 


Formula  (3.1)  is  not  applicable  at  a  large  distance  from  the  wall 
and  because  in  it  an  unlimited  increase  in  y  leads  to  unlimited  growth 
of  wx,  whereas  the  latter  quantity  at  y  -*•  co  is  equal  to  wQ,  Thus, 
this  formula  should  be  used  within  the  framework  of  the  theory  of  a 
boundary  layer  of  finite  thickness. 

Since  in  the  region  £  >0.4  in  a  turbulent  boundary  layer  on 
an  impenetrable  plate  o>  >  0.9,  then  the  indicated  circumstances  do 
not  introduce  appreciable  errors  in  the  calculations  of  friction  and 
heat  transfer. 

In  Fig.  3  is  given  a  graph  of  <p(q),  also  including  data  for  the 
region  of  the  viscous  underlayer.  Distribution  of  velocities  in  it 
is  determined  by  the  expression 

f-n.  (5.2) 

which  corresponds  to  formula  (2.21)  at  f  =  0. 
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Fig.  3.  Universal  velocity  profile 
on  a  flat  plate. 


Intersection  of  the  lines,  calculated  by  the  formulas  (3*1)  and 
(5.2),  gives  the  calculated  thickness  of  the  viscous  underlayer  in  a 
two-layered  diagram  of  isothermal  turbulent  flow  t|^q  =  11.6.  Dimen¬ 
sionless  velocities  on  this  boundary  are  equal  to: 

fit- 1 1,6;  (5.5) 

Such  a  diagram,  which  nominally  divides  the  stream  into  a  viscous 
underlayer,  in  which  p.  »  p.^  and  a  turbulent  nucleus,  in  which  p.  «  p.^. 
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turns  out  to  be  fully  acceptable  for  calculations  of  f^i^tion.  .It 
is  also  applicable  for  calculations  of  hea‘  tranc  for  Pr  ■=  .1 


..  V* , 


‘O’ 


/I 


— -  -  5,5-  2,5  !i:  /? 


wnere 


.  <-■ 


He. 


Further,  we  have 


ft  -  IT, 


V, 


.*,0 

v 


‘1  r.  ■  n 


«*.» 


** 


Law  of  resistance: 


/■ 


e/0 


=  5,5  2,5  In 


;a. 


onvenient  formula  of  T.  Karman: 


(2.5  In  *<••*  4-  3,8)* 


)  It  follows  that 

;  - 

V 

(5.9) 

!.5|/  -**- 

(5.6) 

net-  ,  '  '  5  ,  we  will 

obta  Ln 

fie*' 

(5.7) 

>2-5  1  /  '%■ 

in  (s.5  - 

)  chs.T'P;cs 

*  # 

of  Re  and  in  practice 

uie  tc  use  4 ne  signif 

'leant  ly 

tv  ' 

(3.3) 

From  (3.1)  it  follows  that 

(JkL)L.a,,,„p|, 


•  32 

7" 


(5.9) 


P 


(5.10) 


'me  logarithmic  velocity  profile  is  an  envelope  of  a  family  of  expo  • 
nential  profiles 

?=**•.  (5.11) 

For  many  calculations  the  use  of  such  an  approximation  of  the  velocity 
profile  is  very  convenient. 
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I 


r  t 


By  putting  in  (5.11)  the  values  of  w 

that 

—  =  (*  M)^. 

'/• 


wQ  and  y  =  6,  we  find 


(3.12) 


For  the  thickness  of  displacement  and  loss  of  momentum  we  obtain: 


i*  _  <i 
t  l  +  M  * 

>»*  _  n 

I  <1  +  m)(1+2m)  ’ 

«  =  -£-  =  >+*»• 


(5-15) 


From  this  the  exponential  law  of  resistance  follows: 

(3.1*0 

where 


Q  —  2A  “  »+«  |  0  +  "H1  t  **)  |t  +  n  . 


m  = 


2m 


l  +  M 


(3.15) 


For  the  thickness  of  the  boundary  layer  and  the  relative  thickness  of 
the  viscous  underlayer  we  have: 


(~)  *  M-)*  ; 

(3.16) 

\  *  /o  \  */*  1 

»  .  .  » 

5,,=  11,6*"  . 

(3-17 ) 

The  momentum  equation  takes  the  form: 

fix  2 

(3.18) 

or 

*****  _  *! 

(3.19) 

flUe,  2  * 

Putting  in  (5.19)  the  value  of  cf  from  (5.14)  and  considering 

boundary  layer,  developing  turbulent ly  from  a  certain  cross  section 

x„^,  we  find  that 
cr' 

ft'"' +  ’ - /?<,'  *  "  =  -L jZ-Btfte,- Be,.,  (3.20) 

If  xcr  =  0,  i.e.,  the  turbulent  layer  in  practice  starts  at  the 
leading  edge  of  the  plate. 
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Re"  =(~4 ,  {Z.?D 

to  which  corresponds  the  value: 

(-/„  -  H ,  /?f,‘  m  ,  <  .  •  ■  •  , 

W  M'O 

_  l_  .  ■■  | 

1  i  n  1  '  !  { -2  ■>  v 

j  * 

•Vi  I 

r;»  -  *  .  _ 

1  »  -  /*!  / 

The  values  of  the  coefficients  1'  exponent  ia  ■*  laws  of  resis¬ 
tance  and  the  velocity  distributions  L«/e-  Table  3 .  i .  In  the 

s  :i  ?no  place  are  also  given  certain  other  p.ant  :  •  •  ensuing  from  these 
laws  . 


Table  3.1 
f icients 

(3-14),  ( 

^  /v. 

• 

1  n  Foi 
5*22) 

°  s  >' 

Si''  1  S  ' 

’  *  n 

)  > 

«  1 

iif 

1  1  ' 

Lm. 

* 

8,71 

r»  1 1 
j*  ,  «  . 

lf',6 

!  11,5 

J  ♦* 

i 

0,0975 

i.'/nso 

j  0  0  2.  H 

f  .07*7 

H 

».'» 

1 

1  1.22 

1.20 

m 

Q.'.'.iG  j 

0,222 

0  ?oo 

0,182 

H 

0.0252  1 

0,0206 

1  0.0  90 

0,"|4H 

0.2..0  l 

0.182 

:  0.167 

i  0.154 

0.U5/6  j 

0,1.1') 

1  .1,0362  i 

0,0i"a 

*  *  j.| 

In  practice  in  the  region  of  Re  numbers  <  10  ,  it  is  possible 
u.>  the  formulas  for  the  distributive  law  of  velocities  according 

the  exponent  n  =  1/7. 

4  **  6 

In  the  region  10  <  Re  <  10  the  Faulkner  formula  gives  gooa 

■  suits: 


r/0  =  0,0131  Re*<  "  . 


(5.24) 


For  a  turbulent  layer,  developing  from  a  cross 
formula 

i 

cfQ  -  0,0263  Re 7  7 


section  x  =  0,  the 


(5.25) 


corresponds  to  relationship  (3.24).  In  Fig.  4  is  given  a  comparison 
of  the  given  formulas  with  experimental  data. 
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i 


Cf  to* 


lining  of  a  flat  plates  1)  Blasius'  law  (laminar  conditions);  2) 
the  Prandtl  law;  3)  the  Prandtl-Schlichting  law;  3a)  transition 
region;  4)  the  Schultz -Grunov  law.  Measurements  by  Wlselsberger, 
Hebers,  Frude,  Kempf,  Schenger. 


3.2.  Heat  Transfer  Coefficient  for  f  —  1 
For  polyatomic  gases  Pr  *  1  and  in  the  case  considered  the 
dependency  (1.52)  is  fulfilled. 

Putting  in  it  the  value  of  c^  from  (3.14)  we  have: 

(5.26) 

For  a  plate,  on  which  dynamic  and  thermal  boundary  layers  are 
developing  from  a  cross  section  x  =  0,  the  dependences 

»„  =,  Rt;  "•  (5.27) 

or 

(5.28) 

correspond  to  the  formula  (3.26). 

For  a  region,  in  which  we  will  apply  distributive  law  of  veloc¬ 
ities  by  the  degree  n  =  1/7,  B1  =  0.0576,  and  =  0.20.  Correspon¬ 
dingly,  at  Pr  =  1 

AT*,  =  0,0288  Re°A  (5.29) 
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In  Fig.  5  are  given  the  experimental  data  of  B . 


F.  Petukhov, 


A.  A.  Detlaf,  and  V.  V.  Kirillov  on  the  local  values  of  the  Nusselt 


'  Hupr'0,4 


Local  values  of  Nu 

A 

.  ig  subsonic  s  ti  eamlining 
plate  by  air  (according 
.  u lie  experiments  of  E.  S. 

fi'nov):  O  —  Wq  =  26  0  m/sec; 

1*  —  v  -  24b  m/sec;  0  —  w)  = 

m/sec;  x  —  wQ  =  lo8  m/sec; 
,,  =  162  m/sec ;  +4-  —  Wq  = 

1  ; 1  in/ tec ;  ©  —  v/g  =  62.5 
,  Sec;  ©  —  wQ  =  42.5  m/sec J 
=  54.1  m/sec.. 


number  during  subsonic  streamlinin 
,>f  a  .. .  -  -e  .  i  . 

As  can  b*-  s '  ■  ;  ,  introducti  on 

(  >  l 

o j  factor  Pr  *  '  puts  these 
data  on  the  lJne  cr‘  formula  (5.29) 
'Ini  1  r cum st ante  is  ^ on firmed  by 
the  experiments  of  Emas,  Frank, 
et  al. 

There:;  u  •* ,  for  gases  it  is 
puss  ;\  ’  th  great  accuracy  to 
postulate  that 

Nu,-  -&  P/' Rc\ -  p.;  o) 

to  which  corresponds  the  dependenc 

ReM  (5.51) 

4p 

At  (dp/dx )  -  0  and  Twall  - 
=  const  the  energy  equation  takes 

the  form: 

(3-32) 

dK*x 

Putting  in  (5.52)  the  value  of  StQ 


i rom  (5. 51),  we  find  that  for  the  conditions  examined: 


Bx  Pr~  0,5 

V(l  -  m,( 


(5.55) 


■ y  replacing,  according  to  this  formula,  Bex  in  (5.51)  by  Ret  ,  we 


find  that 


0.<  (I 


(5.54) 


Inasmuch,  as  the  experimental  value  of  the  exponent  for  the 
Prandtl  number  is  determined  for  m  =  0.25,  and  the  number  itself  for 


gases  is  close  to  unity,  it  is  possible  to  consider  with  a  sufficient 
degree  of  accuracy  that 


Stu 


Pr~  0,73  •  Re* 


From  formulas  (5.35)  and  (3.14)  it  follows  that 

*=-*-*-*a(— 

Comparing  formula  (3.21)  and  (5.55)  we  find  that 


(3.55) 


(3.56) 


»r 


(3.57) 


Consequently,  according  to  formula  (2.58)  the  nonsimilarity  coef¬ 
ficient  for  the  velocity  and  temperature  fields  for  a  plate,  com¬ 
pletely  covered  by  dynamic  and  thermal  boundary  layers,  is  equal  to: 

i^Pr™'.  (5.38) 

By  this  formula  for  air  e  0.97#  i.e.,  it  differs  little  from  unity. 

5.5.  Law  of  Resistance  for  a  Nonlsothermal  Boundary  Layer 
Let  us  integrate  equation  (2.10)  for  the  case  of  nonisothermal 
streamlining  of  a  flat,  impenetrable  plate  by  an  unbounded  flow  of 
gas.  Assuming  ~  and  determining  the  ratio  p/pQ  by  formulas 

(2.45),  (2.52),  and  (2.59)#  i.e.,  for  conditions  of  relative  simil¬ 
arity  of  the  temperature  and  velocity  fields,  we  will  calculate  the 
right  integral  of  the  equation. 


For  the  case  6  <  5^  we  obtain: 


<+•-!)**  [ 


arc  sin 


2<+*  —  11 4  + 


V  <  !+•-  1)<++ A+)  +  (.Ai)* 

-«*  .m  -  - 1 

V  «(♦*-■)  (♦•  +  .*)  +  (•**)•  J 


(5-59) 


In  the  case  5  >  6^ 

j iVt'-f- 


1 

f 


4m 


+•--<+*  l)«» 


(3.40) 


_  **"7 


where  o>^  is  the  dimensionless  velocity  at  the  point  y  =  6^. 

For  the  relative  similarity  o>t  =  e  ,  because  at  the  point  y  = 

-  6t  *  -  i- 


Carrying  out  the  integration,  we  find  th* t  at  f  >  5 


t 


i 


<*•  -  I'l 


arc  sin 


2(f  l)« 


-  arc  sin 


I  V  4  <+♦  IK**  -;  i+)  +  <»A+)-. 

2(4*-  -h 


+  arc  sin 


V  4<r  -  »><*•  +  Ai>  <.*{,)• 

arc  sin  «  •>*  -  *  \ 

V  **  V  f.L*  J 


(3.41) 


It  is  necessary  to  consider  that  in  formula  (3.39)  e  <  1,  and 
in  formula  (3.41)  e  >  1.  These  formulas  can  be  used  not  only  for 
theoretical  analysis,  but  also  for  the  experimental  determination  of 
the  functions  Z  and  o^. 


3.4 .  Limiting  Law  of  Resistance 

Assuming  in  formulas  (3-39)  and  (3.41)  =  0  and  Z  -  1,  we  find 

that  at  He  -*  oo: 


a)  at  e  <  1 


*•-11 


trc  sin 


2lv*  -  h  +  «A* 


—  tit  sin 


i  4^*- 1)(**  +  A» +  (•!»• 
_ 


I  4<**  -  !)<*•  +  +  <•**)• 


V 


b)  at  e  >  1 


(3.42) 


_! — fare  sin  —  ^j^±^==  - 
l  /  4(*« -!)(**  +  Ai) +  (•*•>, • 

-  avc  tie  + 


(3.43) 


+•*«•»  y -  «*  *•»  V^r 

Formulas  (3.42)  and  (3.43)  express  the  limiting  laws  of  friction 
for  a  nonisothermal  turbulent  boundary  layer  on  an  impenetrable  plate. 
They  do  not  contain  empirical  "constants  of  turbulence"  and  are  not 
connected  with  any  special  type  of  semlemperical  theories  of  turbulence. 
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The  quantity  cfQ  in  the  limiting  law  is  defined  independent  of 
the  method  of  its  derivation,  i.e.,  can  he  defined  both  from  special 
theoretical  considerations  (for  instance,  proceeding  from  any  sort  of 
semiempirical  theory  of  turbulence  in  isothermal  flow),  and  directly 
from  experimental  data. 

During  adiabatic  flow  the  limiting  ratio  c^/c^q  does  not  depend 

on  the  degree  of  nonsimilarity  of  the  temperature  and  velocity  fields. 

/ 

This  follows  from  formulas  (2.52),  (5.42),  and  (3.45),  which  at  tap  =  0 
give  one  and  the  same  expression. 

( - V^T)  (3.444- 

V-l 

* 

For  f  -*■  1,  i.e.,  for  subsonic  flows,  we  obtain: 


a)  for  e  <  1 


V  =»[ 

2 

T. 

(5.45) 

l 

]• 

t>) 

for  e  >  1 

<r =r. 

a  , 

•  IV. 

(3.46) 

c) 

i 

for  e  =  1  and  V'  m 

«(i'  *  +t) 

l 

■  ]• 

i*  —  /  a 

(5.^7) 

\  M  +> 

) 

The  formulas  obtained  show 

that  the  quantity  e  affects  the 

relative 

change 

of  the  friction  coefficient  due 

to  nonisothermalness 

the  most 

noticeably  during  subsonic 

velocities. 

The  degree  of  this 

influence 

is  seen  from  Table  3.2. 

It  is  interesting  to  note  that  the  influence  of  the  nonsimilarity 
of  the  velocity  and  temperature  fields  on  the  quantity  T  during 
cooling  and  heating  of  a  gas  is  opposite  and  small. 

In  the  majority  of  practically  important  cases  the  ratio  6/6^ 
is  located  in  the  limits  from  0.5  to  2,  which  allows  one  to  assume 
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In  the  formulas  for  cf/cfQ  the  quantity  e  =  1.  This  circumstance 
essentially  simplifies  calculation  of  boundary  layers. 

Table  3.2.  Values  of  (cf/cfo)Re** 


During  Subsonic  Velocities  from  the 
Limit ing  Formulas  (3.45)  and  (3*46) 


1 

•t 

a1 

4  -  o,ar, 

1  ,"3 

1 

:< 

s 

• 

0 

4.00 

9.00 

1 

0,50 

0,33 

0,25 

0.25 

0.26 

0.71 

2.65 

2.35 

1 

0,62 

0.45 

0,35 

0.29 

0.6 

0.90 

1,68 

1.45 

1 

0,65 

0.50 

0,4*. 

0.35 

o.« 

0.97 

1.01 

1.41 

1 

0,67 

0,52 

0,43 

0.37 

1 

1.00 

1.78 

1.38 

1 

0.69 

0,54 

0.45 

0.38 

2 

1.10 

1.69 

1.33 

1 

0,71 

0,58 

0.47 

0,43 

5 

1.26 

1.59 

1.29 

l 

0.75 

0,63 

0,54 

0.49 

10 

1.39 

1.54 

1,16 

1 

0.77 

0,65 

0.58 

0.54 

<30 

3U 

1 

1 

1 

1 

1 

1 

1 

•••The  quantity  e  is  defined  by 
formula  (2.58)  at  n  =  1/7. 


During  supersonic  flow  the  influence  of  e  decreases  with  the 

# 

growth  of  V'  j  i.e.,  the  number  N.  In  general,  the  quantity  e  does 

not  affect  the  limiting  friction  law  during  adiabatic  flow. 

* 

In  Fig.  6  is  shown  the  dependence  of  ¥  on  if/  and  ,  calculated 
by  the  formula: 


f  =  — - — f 


arc  sin 


2(4*  -D  +  A4 


—  arc  sin 


|  4(4“  -  l)(f+  Ai)-  <*«• 

X  * 


uy  I 


(3.48) 


i  4<i-  -  im  ;*  ;  **)+ 
inis  formula  is  from  (3.42)  and  (3.43)  at  e  ■  1. 

As  is  seen,  the  function  Y(^  }  ty)  has  a  complicated  character, 

whereby  the  degree  of  Influence  of  the  heat  transfer  factor  Atf' 

# 

decreases  with  growth  of  the  M  number.  The  upper  limit  for  V  is 
determined  in  this  case  by  the  fact  that  at  M  *  8-10  a  noticeable 
dissociation  of  the  gas  begins  due  to  high  temperatures,  developed 
during  inhibition  of  flow  in  the  boundary  layer. 

Cooling  of  a  gas  (tap  <  0)  causes  an  increase  of  the  resistance 
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Fig.  6.  Dependence  of  Y  on 
if*  and  ^  from  the  limiting 
formula  3.48  (r  =  0.9). 


of  friction  of  flow  against  the  wall; 
heating  of  a  gas  (luff  >  0)  lowers  this 
resistance. 

During  subsonic  flow  f*  *»  1  and, 
according  to  formula  (3.47),  the  ratio 
cf/cf0  a<t  ®  depends  only  on  the 
temperature  factor  i/.  fhe  solution  is 
asymmetric  in  relation  to  heating  and 
cooling.  During  cooling  of  a  gas  the 
temperature  factor  ^  cannot  be  less 
than  zero,  to  which  corresponds  at 
e  ■  1  the  limiting  value  Y  *=  4.  In  the 
region  of  heating  the  growth  of  i/  is 


unlimited  (if  the  problem  of  dissociation  is  excluded)  and  the 
quantity  ¥  can  vary  from  1  to  0. 

In  the  neighborhood  of  ^  -  1  expansion  of  the  right  side  of 
formula  (3.47)  gives 


T  =  *-m  = 


t  +  t 


(5.^9) 


The  quantity 


^  T  I  _  Tft  ***  Tn 


is  the  dimensionless  arithmethic  mean  temperature  of  the  boundary 
layer . 

The  dependences  (3.49)  satisfactorily  approximate  the  exact 
solution  of  (3.47)  in  the  practically  important  range  of  values  of  ? 
from  0.5  to  3. 


3.5.  Limiting  Law  of  Heat  Transfer 
If  5  <  6^,  then  in  the  region  0  <  y  <  0  distribution  of  the 
temperatures  is  determined  by  formula  (2.31).  In  the  region  y  >  5 
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co  =  1,  i.e.,  at  6  <  y  <  6, 


Besides  the  integral 


ii). 


>-)■ 


(5.50) 


(3-51) 


following  from  (2.69)  at  <1  ~  Qq  ariG  p  =  pQ^ri  breaks  up  into  two: 


»t»  (7-)' 


(3-52) 


m\  »> 


Putting  in  (3.51)  the  corresponding  expressions  for  T/Tq  and 
integrating,  we  obtain  an  expression  for  the  limiting  relationships 
of  the  Stanton  numbers: 
a)  for  6  <  &t(e  <  1) 

f  *  * 

*  --  I  arc  sin - — - * - - - 

1  v"  * 1  |/  4— — (L'  -  AC)  «  (iL)» 

-  arc  sin  Bl  - 1 4-  (5-55 


(5.55) 


+  [/  I  +  A+(l  -  ■)  -  »'  I  4  A*  ])! 


b)  for  6>6<t(e>l) 


r* = y~T*H  = 
L  I  4  — {y  *  ** + 

-  arc  sla  -  —  - — - K 

y  1  (V  4  **)  4  (Aiy  j 


(3.5*-: 


During  subsonic  flows,  when  ^  »  1,  the  limiting  relationship  of 

the  Stanton  numbers,  in  general,  does  not  depend  on 
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r'~(  /T+<  )  ’  (5-55) 

The  given  analysis  shows  that  during  the  majority  of  ratios  of 
the  thicknesses  of  the  thermal  and  dynamic  boundary  layers  met  in 
practice  in  a  region  not  close  to  the  point  of  breakaway,  it  is  pos¬ 
sible  to  assume: 

<*.  =  *•  (5.56) 


5.6 .  Comparison  of  the  Limiting  Law  of  Resistance  with 
Experimental  Data  for  Supersonic  Flow 

In  Fig.  7  is  given  a  comparison  of  calculations  by  the  limiting 
formula  (3.48)  with  experimental  data,  obtained  at  fairly  large  values 

of  Re  ,  and  mainly  at  very  high  velo¬ 
cities  and  intense  heat  transfer. 

Not  only  the  qualitative,  but  also 
full  satisfactory  quantitative  agreement 
of  theory  and  experiment  is  clearly 
revealed. 

Thus,  the  limiting  law  of  the 
relative  change  of  resistance  of  friction 
with  nonisothermalness1  sufficiently 
well  describes  real  flows  with  finite 
Re  numbers.  This  result  agrees  well 
with  the  known  experimental  fact  of  the  weak  influence  of  the  Is  number 
on  the  ratio  cf/cf0.  In  this  case  it  has  an  exclusively  important 
value . 2 

Indeed,  the  limiting  formula  (3.^8)  is  obtained  from  the  general 
solution  of  (5.59)  by  means  of  conversion  to  limiting  values  of  the 


Fig.  7.  Comparison  of  cal¬ 
culations  by  the  limiting 
formula  (3.48)  with  experi¬ 
mental  data  for  large  fac¬ 
tors  of  heat  transfer: 

O  —  experiments  of  Lobb, 
Winckler,  and  Persh;  ©  — 
experiments  of  Sommer  and 
Short . 


1  Subsequently,  for  brevity,  the  term  "limiting  law  of  resistance" 
will  be  used. 

Additional  confirmation  was  obtained  in  the  last  work  of  Matting, 
Chapman,  and  Nyholm. 
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quantities  Z  and  o)^,  which  depend  on  the  Reynolds  number.  But  close 
agreement  of  the  limiting  ratios  c^/c^q  with  the  corresponding  ratios 
for  fully  finite  Re  numbers  means  that  the  joint  influence  of  the 
functions  Z  and  is  small.  Therefore ,  it  is  fully  permissible,  for 
calculation  of  the  influence  of  the  Re  number  on  the  relative  change 
of  the  coefficient  of  friction  due  to  the  nonisothermalness  of  flow, 
to  introduce  in  (3.39)  the  second  known  limiting  value  of  these  func¬ 
tions,  corresponding  to  isothermal  flow: 


(5.57) 


For  a  plate  910  *  11.6. 

With  such  a  substitution  we  obtain  the  formula  (at  e  <  1): 


1*  « 


l 


arc  Sin 


2  ft*  -  1)4 


<4*  _  |)  (l  -8,2 1  r/o  )*  I  I  f  Vr*  U  <**  +  **>  <•*»« 

-  .rc  .1. - T. 

V  I )  (i*  +  &i>)  +  (*  J 


(5.58) 


For  adiabatic  flow  ■  0)  we  have 


/ 

1  /  «C  tin  J/ 

(5.59) 

~l\ 

1-8,21'  V  • 

In  Fig.  8  is  given  a  comparison  of  a  large  number  of  experimental 

data  with  calculations  by  formula  (3.58).  Characteristics  of  the 

experiments  are  given  in  Table  3.3.  Experimental  data,  obtained  in 

the  presence  of  heat  transfer  are  reduced  to  adiabatic  conditions  by 

## 

recalculation  by  formula  (3.58).  As  can  be  seen,  change  of  the  Re 
number  from  its  smallest  value,  at  which  a  turbulent  boundary  layer 
can  still  exist  to  infinity  leads  to  a  change  in  the  quantity  f  less 
than  two  times  for  N  *  10.  Essentially,  the  zone  outlined  by  the 
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O  1,0  2,0  3,0  4,0  3,0  0,0  2,0  0,0  0.0  *.i 


* 

Fig.  8.  Comparison  of  experimental  data  with 
formula  (3.59)  (designation  see  Table  3.5). 

Experimental  data  with  heat  transfer  are  re¬ 
duced  to  heat -insulated  conditions  by  equa¬ 
tion  (3.58). 

##  7  ## 

theoretical  dependences  for  Re  *  10^  and  Re  ■  00  embraces  a  zone 

of  scattering  of  experimental  points. 

A  fairly  distinct  tendency  toward  subdividing  the  experimental 

** 

data  according  to  Re  numbers  can  be  noted,  which  corresponds  to 
theory. 

Appearemce  of  certain  groups  of  points  10-20%  from  the  theoret¬ 
ical  value  hardly  can  be  considered  significant,  bearing  in  mind 
the  whole  complexity  of  carrying  out  the  experiments  in  supersonic 
streams. 

On  the  graph  are  placed  also  data  on  heat  transfer.  They  are 
disposed  together  with  data  on  aerodynamic  drag.  Thus,  there  is 
direct  experimental  confirmation  of  the  validity  of  formula  (3*56) 
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for  a  boundary  layer  of  a  gas  in  the  range  of  values  e  close  to  1 

Table  5.3.  Parameters  of  Experiments,  from 
Which  the  Graph  Pig.  8  is  Constructed 

a)  Turbulent  friction  on  a  plate 


Authors 


Coles 


Hill 


Korksgi 


#J  CD 

S* 

£E 


o 

6 


Ke** 


‘wn  11 


iTwnl  1 


I 


2.6 

2.6 

3,7 

4.5 

4.5 

4.5 

4.5 


6600 

10200 

4100 

4560 

2000 

3470 

5M0 

6590 


l 


1.0 

1.® 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 


Lobb, 

Wlnckler, 

Ptrah 


9 

e 

9 

9 

9 

9 

9 


8,99 

9.0* 

9.07 

9.10 

8.22 

8.35 

8.27 

8.29 

8.29 


❖ 

<► 

♦ 

♦ 


9 

G 

G 

9 


5.7*7] 

5.77 

5.71 

5,8051 


1245 

1607 

1900 

2287 

2>«l 

2498 

2885 

3202 

3151 


0.M8 

0,460 

0.47| 

0.495 

0,403 

0.497 

0,500 

0,500 

0,50"? 


2177  i  1.0 

I 

2780  !  1,0 


3429 

4040 


4.« 

53S0 

0.369 

5.01 

6410 

0.713 

5.03 

7960 

0,575 

5,06 

7370 

0,535 

0.33 

8530 

0.613 

6.33 

12610 

0.500 

0.73 

1400 

0.513 

0.7* 

7960 

0.457 

7.67 

1410 

0.465 

1.0 

1.0 


Conditi ons 
of  experi¬ 
ment 


Ff.at- 

insul&ted 

plate 


Cooled  wall 
of  a  coiii- 

cal  nozzle 


Heat- 

insulated 

plate 


Cooled  wall 
of  a  plana 
nozzle 


Method  of  deter¬ 
mination  of  C„ 


Direot  measure¬ 
ment  *<ith  the 
aid  of  a  float 


By  the  velocity 
gradients  or. 
the  wall 


Direct  measure¬ 
ment  with  the 
aid  of  a  float 


by  velocity 
gradients  or 
the  wall  and 
via  ths 
Reynolds ’  * 

analogy  by 
measurements 
of  heat 
fluxes  in  the 
inveetl  gated 
Motion 
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Table  3.3a  Continued 


Conditions 
of  experi- 
msrrt 


Heat- 

insulated 

plat* 


Streamlin¬ 
ing  of  a 
heat- 

insulated 
oylinder 
In  the 
longitudi¬ 
nal  direc- 
tion 


Shot  by  a 
hollow 
oylinder  in 

a  wind 
tunnel 
toward  the 
flow 


Cylindrical 

■urfaoe 


Method  of 
determination 
of  C 


By  changes  of 
velocity  pro¬ 
files  and 
momentum 
equation 


Direct  measure¬ 
ment  of  mean 
values 


Direot  measure¬ 
ment  of  local 
values 


B|y  measurements 
of  the  velocity 
of  flight  of  a 
model 


By  measurements  of 
velocity  profiles 


By  veloolty 
profiles 


e  on  plate 


Pappas 


2,516 
It  1*0 
3,410 


1,123 

2.210 


Streamlining 
of  a  cone 


Plate  with 
uniform  pre¬ 
heating 


Method  of  deter¬ 
mination  of  St 


Measurement  of  heat 
transfer  by  a  non¬ 
steady  state  method 


By  eonsiaeption  of 
electric  power  and 
wall  teaperature 


Table  3.5b  Continued 


Authors 

Point*] 
(Fig-  8  ] 

M 

Conditions  of 
experiment 

Method  o*'  u.  t  w- 
rri nation  of  St 

■e 

2.00 

Plats  with 

By  consumption  of 

Shulberg 

2,So 

309 

uniform  pre¬ 
heat 

heat  transfer  t.y 

•».  r.or.s4  3*'ite 

i.,ethu(i 

Slack, 

Ml  it 

~  i 

fata  not  ci:*<;  1 

3.7.  Plate  with  Initial  Adiabatic  Section 
A  diagram  of  the  problem  is  shown  in  Fig.  Q.  A  dynamic  turbulent 
boundary  layer  is  developed  from  the  leading  edge  of  plate.  On  the 
section  with  length  Xq  heat  transfer  is  absent.  From  section  x  =  xQ 
heat  transfer  begins  between  the  plate  and  the  gas.  The  forming 
thermal  boundary  layer  is  submerged  in  the  dynamic  layer,  i.e.,  5^  < 

<  5. 


Fig.  9.  Plate  with  initial  heat- 
insulated  section  (Pr  =  O.72). 


Assuming  Tr  =  1,  p  =  p Q,'Pt  =0,  t^=  l ,  q  «  qQ,  and  5  >  6t,  we 
bring  equation  (2.64)  to  the  form: 


#»  *• 


(5.60) 


#jr  dy 

On  the  other  hand,  for  these  conditions  it  follows  from  (2,6)  that 

'-Sr-  V'-t;*  ■  (5.61) 

Combining  these  equations  and  integrating,  we  obtain  the  dependence 
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IM  *  *f 


*J  V  * 


*1- 


!»• 


(5.62) 


From  (3.61)  we  have 


•  -f  '* 


(3.63) 


Consequently, 


where 


2  * 


!•»  —  •,«)=!  —  I, 


IM 


(-H- 


(3.64) 


(3.65) 


The  dimensionless  temperature  difference  In  the  viscous  under¬ 
layer  : 


3-  >  1  ■  -  —  St ^  11 1_  "  • 


l.e. , 


2Mt 


I..  =  ”  Pr* 


(3.66) 

Putting  this  value  of  in  (3.64),  we  find  that  at  Pr  *  li 

(3.67) 

In  a  first  approximation,  for  n  =*  1/7  and  PT  =  1, 


2  st  _  _|_ 
«>  •» 


_**r_ 

**•* 


0,0129 

0.0IW 

He* 


(3.68) 


Let  us  integrate  the  first  of  these  equations,  assuming  that  at 
x  »  xQ  6t  =  0.  Let  us  integrate  the  second  equation,  assuming  that 
at  x  =  0  6  b  0,  As  a  result  we  obtain  relationship  for  the  thicknesses 
of  the  boundary  layers: 


_*l_  =  i  x~x"  \M 

h  \  X  f 


(3.69) 
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4-fcfc.  I  • 


MM  1 


Correspondingly , 

**  =  i  *  r**  (5.70) 

As  can  be  seen  from  Fig.  10,  this  formula  is  well  confirmed  by  cxpcr  - 
ment. 


Fig.  10.  Heat  transfer  to  a 
plate  with  initial  heat- 
insulated  section:  -  cal¬ 

culation  by  equati  on  (5.70); 

O  —  calculation  by  the  equa¬ 
tion,  describing  the  experi¬ 
ments  of  W.  Reynolds,  W.  Kays, 
and  S.  Kline. 


_  ## 

Since  in  this  case 


(Va)*«  * ,  then 


or 


St_  _  /  X  -  Xu  V0  "** 
St,  \  x  } 


(J.71) 


(3-72) 


where 

st,  mOflm  if#;*  ** . 

A  second  approximation  slightly  refines  this  result. 

The  correction  for  isothermalness  can  be  calculated  by  the  for¬ 
mulas  given  earlier  upon  substitution  in  them  of  the  value: 

•=K(-rhr)T 
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IM  tt 


*  .. 


t* 


Table  3.4.  The  Ratio  0t/8to  by 
Formula  (3.72) 


X. 

1 

u* 

Ul 

i.e 

in 

* 

r 

m 

A 

0 

0.S1 

O.srj  0.91 

0.91 

0.9I| 

|  9,95 

1 

For  quaall8othermal  conditions  it  is  possible  to  assume: 

_  i.OI»  I  X-Xm  Y~ 

*>"*;*''  \  *  / 

Hence,  for  a  plate  with  a  completely  turbulent  boundary  layer: 


(3.74) 


Correspondingly, 


or 


*< 


_  0.9156 


st = o.ons  tr  j 


(3.75) 

•H- ) 

(3.76) 

.1  *  \ 

(3.77) 

r\  *  /’ 

Here 


«• 


J; 


(3.78) 


=  -f(^-  I 

The  function  qp(x/xQ)  is  represented  graphically  in  Fig.  11. 

At  (x/x0)  >1.50  with  an  error  not  exceeding  8#,  it  is  possible 
to  assume: 


*•4  =  0,0288  Pi*'**'?, 


(3.79) 


i.e.,  to  calculate  heat  transfer  by  the  usual  formula  with  substitution 
in  it  of  the  true  length  of  the  heated  section: 

x,  =  jr  —  x„.  (3.80) 
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Formula  (3.79)  was  offered  by  M.  A.  Mikheev  on  the  basis  of 
treatment  of  the  experimental  data  (Fig.  12). 


Fig.  12.  Comparison  of  experimental  data 
for  heat  transfer  to  a  plate  with  formula 
(5*79) •  Dotted  lines  —  calculation  of  the 
initial  laminar  section. 
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3.8.  Solution  of  Equations  of  Momentum  and  Energy  for 
a  fronisothermal  Boundary  tayer  for 

Twall  '  con8t 

The  analysis  conducted  showed  that  the  difference  in  the  thick¬ 
nesses  of  the  dynamic  and  thermal  boundary  layers  slightly  affects 
the  relative  change  of  the  coefficient  of  friction  in  connection  with 
the  nonlsothermalness  of  the  flow.  This  Important  result  allows  in  a 
majority  of  the  practically  interesting  cases  to  assume  in  the  expres¬ 
sion  for  T  the  quantity  e  »  1. 

Just  as  important  is  the  fact  of  the  weak  influence  on  the  T 

Reynolds  number.  The  latter  allows  introduction  in  the  equations  of 

## 

momentum  and  energy  the  quantity  f ,  referred  to  the  Xt  number,  aver¬ 
age  in  the  considered  section,  or  simply  the  limiting  law,  strictly 
self -simulating  relative  to  the  lie  number. 

Both  these  circumstances  are  the  result  of  the  small  value  of 
the  exponent  in  formula  (3.11),  i.e.,  high  population  of  the  velocity 
profiles  in  turbulent  flow. 

Let  us  consider  nonisothermal  turbulent  layers,  developing  from 
the  leading  edge  of  a  plate  at  T  -  const.  We  have 


Introducing  in  these  expressions  the  values  of  cfQ  and  §tQ 


out  integration  for  the 

conditions 

f  - 

const,  x  ■  0,  Rt 

we  obtain 

(-*-)  =1 

f  i 

i 

\  =**!♦« 

l  </.  m*M  1 

/**, 

(3.81) 

(3.82) 

and  carrying 

_  **  . 

1,0 

(3.85) 


Thus,  the  ratio  of  the  coefficients  of  friction  and  heat  transfer 
for  identical  values  of  Rax  numbers  more  weakly  depends  on  the  non- 
isothermalness  than  for  identical  values  of  M  and  . 


-53- 


¥ 


t 


s 

4- 


f 


i 


In  Pig.  15  1b  given  a  comparison  of  calculations  by  the  theoret¬ 
ical  formulas  (3.58)-(5.85)  with  calculations  according  to  a  method, 
based  on  the  semiempirical  Prandtl-Karman  theory.  At  the  same  plac- 
are  placed  experimental  points.  The  known  fact  that,  at  large  M 
numbers  Bind  Intense  heat  transfer,  the  contemporary  semiempirlcal 
theory  does  not  give  satisfactory  results,  is  clearly  confirmed. 


Fig.  15.  Comparison  of  experimental 
data  with  calculations  by  various 
methods:  1)  Van  Dr test’s  method;  2) 
by  formula  (5.85);  •  —  experiments 
of  Sommer  and  Short. 
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CHAPTER  IV 

LONGITUDINAL  STREAMLINING  OF  A  PERMEABLE  PLATE 

Definitions  of  Cyrillic  Items  In 
Order  of  Appearance 

ct  =  wall  -  wall 
up  *  cr  =  critical 
t  «  t  ■  thermal,  turbulent 
was  =  rad  ■*  radiated,  radiation 
oth  «  rel  *  ralation,  relative 

4.1.  Formulation  of  the  Problem 
The  problem  of  calculation  of  a  boundary  layer  on  a  surface, 
penetrated  by  a  flow  of  substance,  has  an  extraordinarily  Important 
significance.  Such  processes  arise  during  protection  of  parts  of 
machines  from  the  influence  of  high-temperature  gas  streams  (so-called 
"porous  cooling"  of  blades  of  gas  turbint-i,  the  enclosing  surfaces  of 
combustion  chambers,  and  so  forth)  during  evaporation  and  condensation, 
the  presence  of  chemical  reactionr  1  *  the  flow  and  at  the  wall,  hard¬ 
ening  of  a  liquid,  and  fusing  of  solicit* . 

Below  is  considered  a  turbulent  boundary  layer  of  a  gas  on  a 
surface  permeable  at  all  points.  If  gas  is  injected  into  a  boundary 
layer  or  is  pulled  from  it,  then  the  number  of  openings  in  the  wall  is 
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great,  and  their  dimensions  are  small.  Therefore,  the  temperature  cf 
the  penetrating  gas  should  he  equal  to  the  wall  temperature.  A 
diagram  of  the  flow  problem  is  shown  in  Fig.  14. 


Mass  velocity  of  a  strt  t* 

of  gas  intersecting  a  wall, 

(4.1) 


Fig.  14.  Diagram  of  a  Boundary  Layer  Inside  the  boundary  lay^r 

on  a  Permeable  Plate. 

this  flow  causes  an  additlona 

component  of  the  velocity  vector  w  and  gradually  it  is  dispersed  so 

e 

that  .t  th.  wall  and  in  the  region  y  *  6  transverse 

flow  of  the  injected  substance  is  equal  to  zero. 

Transverse  flow  of  the  substance  causes  the  turbulent  viscous 
underlayer.  This  circumstance  favors  application  of  the  method  of 
conversion  to  the  limiting  laws  of  friction  and  heat  transfer,  corre¬ 
sponding  to  Be  -*  oo. 

In  connection  with  the  above  further  research  is  built  on  the 
basis  of  equation  (2.30). 

Distribution  of  tangential  stresses  is  determined  by  formula 
(2.40)  at  A  «  0. 

Ve  have 


(4.8) 


This  integral  has  a  finite  value,  i.e.,  growth  of  the  integrand 
should  be  limited.  Since  the  quantities  a>  and  f(£)  lie  in  the  interval 
from  0  to  1,  then  changes  of  the  integrand  in  this  case  are  connected 
with  the  change  of  the  permeability  factor  b.  During  injection  of  a 
gas  through  the  wall  into  the  boundary  layer  the  quantity  b  is  positive. 
Consequently,  there  should  exist  some  limiting  value  of  the  permeability 
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factor  b  ■  bQ  ,  at  which  the  value  of  f  turns  out  to  be  equal  to  zero. 
This  phenomenon  can  be  identified  with  separation  of  the  boundary  layer 
from  the  streamlined  surface. 1 

The  critical  value  of  the  permeability  factor  will  be  determined 
by  equation 


i 


(4.3) 


which  ensues  from  (4.2)  at  T  -  0. 

Subsequently,  we  will  be  limited  by  solutions  of  equations  (4.2) 
and  (4.3)  in  an  approximation,  corresponding  to  the  conditions  f(()  - 
«  1  and  e  *  1. 

Theoretical  calculations,  and  the  principal  comparlsion  with 
experimental  data,  show  the  full  acceptability  of  these  conditions  for 
the  majority  of  practically  important  cases. 

Putting  in  (2.69)  the  value  of  V,  •  from  formula  (2.76),  we 

f* 

obtain  an  expression  for  the  limiting  law  of  heat  transfer  on  a 
permeable  surface: 


|r,  + 

Since  f(4t)  *  f($)#  then  solution  of  equations  (4.2)  and  (4.4) 
at  e  *  1  have  the  same  form,  l.e., 

**,(*»)  =  *<*)•  (4.5) 

The  gas  introduced  through  the  wall  can  differ  from  the  substance 
of  the  main  flow;  therefore,  we  will  distinguish  uniform  and  nonuniform 


1For  greater  detail  see  Chap.  V. 
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boundary  layers*  A  nonuniform  boundary  layer  appears  during  Injection 
of  a  gas,  differing  from  the  gas  of  the  main  flow.  In  this  case  the 
boundary  layer  consists  of  a  mixture  of  gases,  whereby  at  the  wan  the 
concentration  of  Injected  gas  has  the  largest  value,  but  at  the  external 
boundary  of  the  layer  Is  equal  to  zero. 


4.2.  Law  of  Resistance  for  a  Uniform. 

Isothermal  Boundary  Layer 

Assuming  In  equations  (4.2)  and  (4.3)  p  =  p0  and  f(£)  -  1,  we  find 

that 


▼HI- <US*)*; 


(4.6) 


(4.7) 


Correspondingly,  formula  (4.6)  can  be  written  In  the  form: 


(4.8) 


As  will  be  shown  below,  this  simple  formula  possesses  great  universal¬ 
ity. 

The  moment  In  equation  for  gradientless  streamlining  of  a  permeable 
plate  has  the  form: 


fcfjfe?  _  c; 
tLx  2 


(4.9) 


or 


Let 

boundary 

For 


=  (*+%•) 


rf*#,  ' '  2 

us  Integrate  equation  (4.10),  assuming  that 
layer  Is  developed  from  the  leading  edge  of 
the  case  of  b  ■  const  we  obtain: 


(4.10) 

the  turbulent 
plate. 


=  ▼<*  -r  1*) 


*  *i  . 


(4.11) 


1 
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VW»H»i  wnwin 


rm’.+Ht* 


-  wli  yxijm—m/w  CMT'W* 


i» 


When  the  quantity  ?  is  determined  by  (4.6), 


(4.12) 


Here  m  is  the  exponent  in  formula  (3.14). 

At  i  >  0,25  and  b  =4,  from  (4.12)  it  follows  that 

Ci 

=  <•  -  0.25*)’(l  i  (4.13) 


For  the  case 
awkward ,  but  at  m 


Pwallwwall 

Powo 


=  const  the  solution  obtained  is  very 


-  0.25  Is  well  approximated  by  formula 


MM-  =  <l-0,25*)*(l+0.25*:-« 
\  tfn  R*. 


(4.14) 


Formulas  (4.13)  and  (4.14)  agree  with  an  error,  not  exceeding  2#. 

Cf* 

In  Fig.  15  is  given  a  comparison  of  the  values  of  compared 

cf0 

for  identical  values  of  Re**  and  R#x.  The  corresponding  curves  are 
close  to  each  other. 


In  Fig.  16  is  shown  a  comparison  of  the 
theoretical  formula  with  experimental  data. 

It  is  possible  to  ascertain  full  satisfaction 
of  theory. 

Due  to  the  large  slope  of  the  relation 
Y(b)  in  the  pre -breakaway  region,  experi¬ 


Fig.  15.  Dependence 
of  the  coefficient  of 
friction  on  the 
permeability  factor: 

1)  -=-  for  Re**  =  idem; 
cf0 

Cf* 

2)  7T=-  for  Re„.  =  idem. 

CfO  x 


mental  determination  of  the  point  of 
breakaway  is  inaccurate. 

In  practice  in  e3q>eriments  the  condition 
of  breakaway  of  the  boundary  layer  for  the 
case  considered  can  be  fixed  already  at 
values  of  b  of  the  order  of  3. 


-59- 


S> - 


i.iwtrft. ’ t  d  k  1 


Fig.  16.  Comparison  of  calculations  by  formula 
(4.14)  with  the  experiments  of  Hacker,  Mickley, 
Pappas,  and  Okuno, 

i 

O  —  experiments  of  Hacker  -O —  experiments  of 


Mickley 


(ao) 


©  —  experiments  of  Mickley 


C4)1 


—  experiments  of 
Pappas  and  Okuno. 


Since  the  quadratic  term  in  formula  (4.8)  begins  noticeably  to 
have  an  effect  only  in  the  range  of  values  of  the  permeability 
parameter,  close  to  the  critical,  then  for  b,  not  very  close  to 


b  , 
cr 


v^t-  2 


(*.15) 


The  linear  relation  between  the  quantities  ¥  and  b  was  noteu  by  several 
experimenters  and  is  well  explained  by  the  theory  expounded.  At 

bcr  ■  4  “d  b  «  bcr 


▼•1-0,5  6. 


(4.16) 
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4 . 3.  Law  of  Heat  Transfer  for  a  Uniform, 
Quaslisothermal  Boundary  Layer 


For  small  temperature  differences,  when  f  *>  1,  the  physical 
properties  of  a  gas  in  a  boundary  layer  can  be  considered  constant. 

By  formulas  (4.5)  and  (4.8)  the  limiting  law  of  heat  transfer  for 
a  quasiisothermal  boundary  layer  at  a  permeable  wall  has  the  form: 

(4.17) 

In  addition: 

*'*  =  *•*  (4.18) 


In  the  presence  of  relative  similarity  of  the  velocity  and 
temperature  fields  it  is  possible  to  apply  formulas  (3.36)  and  (3.37). 
Then 


b, 


bPt 


(4.19) 


During  development  of  both  boundary  layers  from  the  leading  edge 
of  the  plate  condition  (3*38)  is  fulfilled,  and 

t>,*bPr'\  (4.20) 


4.4.  Calculation  of  Cooling 

Usually  in  the  calculation  of  a  boundary  layer  on  a  cooled  wall 
distribution  of  the  temperature  Twan(x)  and  the  initial  temperature 
of  the  liquid  coolant  are  assigned.  It  is  necessary  to  determine 
the  flow  rate  of  this  liquid  through  each  cross  section  of  the 
protected  surface. 

The  quantity  of  heat,  transmitted  by  the  main  flow  to  the  wall, 
in  a  given  section  is 

(ft, =  1  ( 7,i  T„)  -  fiuj. 

where  qrad  is  the  heat  flow  due  to  radiation. 


(4.21) 


Let  us  assume  that  part  of  this  heat  q^  is  transmitted  through 
the  wal3*  not  in  connection  with  introduction  through  it,  of  the  cooling 
medium  (for  instance,  by  means  of  thermal  conduction  through 
structural  metal).  Then  the  quantity  of  heat,  which  the  cooling 
medium  should  receive  is: 

~  Vn  —  Vi  -•  S  cpiJ t  ( 7*ct  r,)/  / 1.  no\ 


where  J  ^  -  Piwl  ^ow  of  the  maBS  of  cooling  medium  through  the 

wall.  Combining  these  equations,  we  find  that 


v.= 


_  *r 


rt 


'CT 


(4.23) 


where 


j . .  i 


*HT4  - 


♦CT 


On  the  other  hand,  for  ?  we  have  the  dependence  (4.17).  Combining 

s 

equations  (4,17)  and  (4.23)  and  solving  them  relative  to  the  value  of 
the  thermal  permeability  factor  of  the  wall,  we  obtain  the  value 


t-'  »  -m-mI- 


(4.24) 


where 


Au 


/»  =  0  when  /f>0  *f<V 

4.5.  Laws  of  Resistance  and  Heat  Transfer  for  a 
Uniform,  Nonisothermal,  Subsonic  Boundary  Layer 

In  a  uniform  subsonic  boundary  layer  of  a  gas  e  =  1  according 

to  formula  (2,52) 
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P  <t 


(4*25) 


Po 


Putting  this  value  of  in  equation  (4.2)  and  solving  it  at 
f(£)  =  i,  we  find  that: 
a)  at  if/  <  1 


=  4  fin  ('-»><■  +  ».>  +» P, 

(l”*  ».  [  /l-*  +,/  J 


(4.26) 


b)  at  f  >  1 


^  |/^ ,T-7  *i,vr,fClt/ &]' 


(4.27) 


At  ^  =  1  we  have  formula  (4.6). 

Correspondingly,  the  critical  value  of  the  parameter  of 
permeability  of  the  walls,  calculated  by  equation  (4.5)  upon  substi¬ 
tution  in  it  of  the  value  of  ~  from  (4.25)  and  f(£)  «  1,  are  determined 
by  the  formulas: 
a)  at  <»  1 


(4.28) 


b)  at  >  1 


'•-Trl 


arc  cot 


IzjL 


•)*• 


(4.29) 


Table  4.1.  Values  of  b„  by  formulas 
(4.28)  and  (4.29).  or 


♦ 

04S 

MO 

0,75 

t 

ia> 

4.00 

*<5P 

9.4 

6.2 

4.S 

4 

2.5 

1.4 

/JJSuu) 

\  /6r 

2.3 

3.1 

3.6 

4 

5,0 

5.7 
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The  relation  b  ($/)  is  shown  graphically  in  Fig.  17. 

C  i 


Fig.  17.  Influence  of  heat  transfer  on  the 
parameter  of  breakaway  b  for  a  uniform 

boundary  layer.  . — -  Calculation  by 

equation  (4.28)  for  the  case  if  <  1  and  by 
equation  (4.29)  for  the  case  ip  >  1, 


It  is  interesting  to  note  that  in  the  region  ip  >  1  the  critical 

relationship  of  velocities  changes  significantly  weaker  than  the 

critical  relationship  of  mass  flow  rates.  In  the  region  ip  <  1  both 

parameters  change  almost  equally,  although  there  is  observed  a 

/wwall\ 

somewhat  greater  conservation  of  the  ratio  ( - - 1  . 

\w0  /cr 

For  engineering  practice  the  case  f  <  1  has  the  greatest 
significance  since  in  this  case  introduction  of  a  substance  through 

the  wall  into  the  boundary  layer  protects  the  streamlined  body  from 
the  thermal  influence  of  the  main  flow  of  gas. 

In  Fig.  18  is  given  a  comparison  of  calculations  by  formulas 
(4.26)  and  (4.27)  with  calculations  by  formula 


( 


/T  +  1 


(4.50) 
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As  can  be  seen,  this  simple  combination  of  formulas  (3.47)  and 
(4.8)  well  approximates  the  exact  solution. 


|  Calculation  by  equation 

1 

(«•) 

W.W) 

*  0,7  0,4 

0.6 

0.8 

2 

4 

8 

Designation  0  3 

0 

0 

0 

In  Fig.  19  is  given  a  comparison  of  calculations  by  formula  (4.30) 
for  heat  transfer  on  a  plate  and  in  the  inlet  section  of  a  pipe  with 
experimental  data.  In  spite  of  a  significant  scattering  of  the 
experimental  points,  they  are  grouped  in  a  mass  around  a  theoretical 
line. 
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o  at  0.4  as  at  v> 


Fig.  19.  Influence  of  injection  of  a  gas  on 

convective  heat  transfer.  - calculation 

by  equation  (4.30),  A  —  experiments  of  Mickley 
(plate);  experiments  of  Friedman  (pipe). 
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4.6.  Law  of  Resistance  for  a  Nonuniform, 

Isothermal  Boundary  Layer 

During  supply  through  a  wall  into  a  boundary  layer  of  a  foreign 
gas  the  process  of  diffusion  arises.  The  partial  density  of  the 
injected  gas  is  changed  from  the  value  p  ^  on  a  streamlined  surface 
to  zero  in  the  region  y  =  5^,  where  6^  is  the  thickness  of  the  diffusion 
.boundary  layer.  The  field  of  concentrations  can  be  connected  with  the 
velocity  field  by  a  formula  of  the  type  (2.44).  For  the  flow,  not 
close  to  the  point  of  breakaway,  it  is  possible  to  assume  the  relation1 


fa-*' 


Pc* 


(4.31) 


1In  the  general  case,  for  instance,  during  evaporation  of  water  in 
humid  air,  Pq  ^  0  and  in  the  denominator  of  (4.31)  one  should  write  the 
difference  pwall  -  pQ. 
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where 


(4.32) 


Since  in  gases  the  diffusion  Prandtl  number  is  close  to  unity, 
then  all  the  conclusions,  made  in  chapters  II  and  III  about  the  weak 
influence  of  the  degree  of  nonsimilarity  of  the  temperature  and  velocity 
fields  on  the  relative  changes  of  the  coefficients  of  friction  and 
heat  transfer,  are  also  valid  for  the  process  of  diffusion.  Therefore, 
for  a  plate  it  is  possible  to  consider  the  limiting  law  of  friction, 
assuming  «  e  »  1. 

The  gas  constant  of  the  mixture  is  connected  with  the  density  of 
the  mixture  and  partial  density  of  the  injected  gas  p'  by  a  known 
relationships 


(4.33) 


Here  R,  R»,  RQ  are,  respectively,  the  gas  constants  of  the  mixture, 
injected  gas,  and  basic  gas. 

On  the  other  hand,  the  density  of  the  mixture  p  is  connected  with 
the  density  of  basic  gas  outside  the  boundary  layer  pQ  by  the  equations1 


* 


(4.34) 


i.e.,  under  isothermal  conditions 

f _ _  /?« 

*  "IT-  (4.35) 

Combining  formulas  (4.31) >  (4.33) >  and  (4.35) ,  we  find  that  at  e  =  1 
and  T  =  const 

_JL=e_fa_+(|_^.)„.  (4.36) 


1  Since  in  the  boundary  lever 
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Introducing  this  value  of  in  equation  (4.2),  we  find  that: 

a)  at  R*  >  Rq  (is  injected  with  a  small  molecular  weight  (j. *  <  |i0) 

.  ln/M.  +/ 

»  "m*»-  >)  /»T+/  +1-1 

b)  at  R*  <  Rq  (gas  is  injected  with  a  large  molecular  weight 
H*  >  M-q) 


i  (4.37) 


-  1  -  ;  t' 

*1  l  V  (I -♦.)*, 


x[„„  V/iJ-J' 


(4.38) 


where 


♦«  =  -£-•  (4.39) 

is  a  quantity  to  a  certain  degree  analogous  to  the  temperature  factor 

In  a  uniform  flow  of  gas,  obeying  the  Clapeyron-Mendeleev 

equation  ^ 

Definition  of  the  parameter  depends  on  the  method  of  supply  of 
the  foreign  gas  to  the  boundary  layer. 

If  a  gas  is  introduced  through  a  porous  wall,  then  the  flow  of 
substance  is 


where  D 


D-br), 

ry  is  the  diffusion  coefficient. 


From  (4.31)  it  follows  that: 


fc* 

(in) 

\  *y  L  0 

•n 

v  *>  /. 

(4.40) 


(4.41) 
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On  the  other  hand. 


Hence, 


where  PrD  «  Is  the  diffusion  Prandtl  number. 

From  these  equations  it  follows  that  at  Prjj  *»  eD  1 

■*CT  _ 

fct  I  + 


Taking  into  account  (4.33) »  we  determine  the  quantity 


(4.42) 


(4.44) 


(4.45) 


In  Fig.  20  is  given  a  comparison  of  calculations  by  formulas 
(4.37) ,  (4.38) ,  (4.45)  with  experimental  data.  The  theory  well 
describes  the  results  of  the  experiments  with  such  heterogeneous  pairs 
of  media,  as  helium-air  and  freon-12-air. 

At  b  =  bcr  separation  of  the  boundary  layer  occurs  and  the  surface 
of  the  wall  is  covered  by  a  film  of  injected  gas.  Thus,  in  this  case 
Pwall  =  Pwall  and  ^or  condi't^on  T  =  const 


(4.46) 


Putting  this  expression  in  equation  (4.3)  and  carrying  out 
integration  at  f(4)  =1,  we  find  that, 
a)  at  R'  >  Rq 


(4.47) 
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b)  at  R*  <  RQ 


Fig.  20.  Influence  of  injection  of  a  foreign 
gas  on  the  coefficient  of  friction  on  a  plate. 

—  —  calculations  by  equations  (4,6)  —  (air-air), 

(4.37)  —  ( helium -air ) ,  and  (4.38)  —  (Freon-12-air) , 

O  —  experiments  of  Hacker  (air-air),  8  —  experiments 
of  Mickley  (air-air).  Experiments  of  Pappas  and 
Okunoj  O  —  air-air,  9  —  helium-air,  9  —  Freon-12-air. 


In  Fig.  21  is  shown  the  dependence  of  the  critical  value  of  the 
parameter  of  wall  permeability  on  the  relation  of  the  molecular 
weights  of  the  injected  gas  and  basic  flow. 

In  Fig.  22  is  given  a  comparison  of  the  experimental  data  with 
calculation  by  formula  (4.12)  during  determination  of  the  quantity  b 

v  X 

by  formulas  (4.47)  and  (4.48). 

In  that  case,  when  the  streamlined  surface  is  covered  by  a  film  of 
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Fig.  21.  Dependence  of  the 

critical  parameter  of  injection 

b  on  the  ratio  of  the  molecular 
cr 

weights  (Y  =  1):  a)  l)r—  by 
equation  (4.47),  2)—  by  the 
equation 

*■9  =  I  +  3—  ; 

b)— by  equation  (4.48).  In  the 
range  of  from  1  to  10  is  the 

H0 

linear  approximation  is  valid: 


liquid  and  evaporation  occurs,  the 
value  of  the  density  of  the  mix¬ 
ture  at  the  wall  is  determined 
by  the  total  pressure  and  partial 
pressure  of  the  volatile  liquid. 
The  latter  with  a  great  degree  of 
accuracy  is  equal  to  the  saturated 
pressure  at  the  evaporation 
temperature.  In  addition: 

♦,=- - *- 


-(-r)4 


*  (4.49) 


where  pQ  is  the  pressure  outside 
the  boundary  layer. 

For  quasiisothermal  condi¬ 
tions,  when  ip  **  1, 

i 


♦.= 


(4.50) 


-  1.47  -  2,53*-, 


It  is  also  possible  to  use 
these  formulas  for  calculation  of 
heat  transfer  during  quasiisother¬ 


mal  conditions  with  replacement  of  the  factor  b  by  the  factor 


t»\i\ 


(*.51) 
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£kTu* 


o  a*  v  as  Q4  as  as  v  as  as  to** 

Pig,  22.  Comparison  of  exper¬ 
imental  data  with  formula  (4.12), 
O  —  experiments  of  Hacker,  $  — 
experiments  of  Mickley,  Exper¬ 
iments  of  Pappas  and  Okuno;  (D  — 
air-air,  0  —  helium-air,  Q  — 
Freon-12-air, 


4.7*  Law  of  Resistance  for 
a  Nonuniform,  Nonisothermal, 
Subsonic  Boundary  Layer 

During  nonisothermal  flow 
the  intensity  of  the  process  of 
diffusion  is  more  correctly 
determined  by  the  difference  of 
partial  pressures.  Let  us  assume 
that  the  following  condition  is 
fulfilled1 


Pcr-P' 


(4.52) 


which  at  T  =  const  coincides  with  (4.31). 
Taking  into  account  (4.33) ,  we  obtain: 


-N'-'-f  (•-  -£)» 


(4.53) 


At  the  wall 


“  T- [•- 4- 


(4.54) 


Hence, 


(4.55) 


Determining  the  temperature  by  formula  (2.52)  at  =  1  and 


^■In  the  more  general  case  in  the  denominator  one  should  write  the 
difference  P^all  -  Pq. 
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r  » 


fc&  i# 


e  *  1,  we  have 


~  ±±±y 


♦  +  d-*)- 


(4.56) 


In  these  formulas  the  factor  ip^  is  defined  by  (4.59).  Introducing 
this  value  of  -P—  in  (4.2)  and  (4.5)  and  assuming  f(£)  =  1,  we  obtain 

Po 

the  equation: 


r  *  1 

+.  + 

l'~  4,  )* 

1-  Ml 

+)  HO  +  *i'-> 

_*<>  \ 

«■  +l' 

*  r 

1+  +  (i 

-  m|  III 

d  in 


(4.57) 


(4.58) 


During  derivation  of  equation  (4.58)  it  is  taken  into  account 


R * 

that  at  the  point  of  breakaway  Pwall  =  p^all  and  • 

In  Fig.  25  are  given  the  results  of  the  numerical  solution  of 

equation  (4.58)  in  the  range  of  values  of  from  0.1  to  1.0  and 

K  0 

values  of  ip  from  0,2  to  2.0.  Along  the  axis  of  ordinates  on  the  graph 

R  * 

is  placed  the  ratio  of  the  quantity  b_„  at  given  values  of  ip  and  — — 

cr  rq 

p  i 

to  the  quantity  bcr  at  the  same  value  of  ip  and  =  1  (i.e.  is 

calculated  by  formula  (4.28)  at  ip  <  1  and  by  formula  (4.29)  at  ip  >  1). 
These  results  are  shown  the  graph  by  the  points.  The  line  given  by 
these  points  is  described  by  the  formulas: 

at  u’  >  M«0 

^,,(0,37  +  0,67-4-);  (4<59) 


at  u1  <  M-o 


(o,2S  +  0,75  )  • , 


(4.60) 


'I  J 


In  these  and  the  other  formulas 

M.*  Is  the  molecular  weight  of 

injected  gas  and  u0  the  molecular 

weight  of  the  gas  of  the  basic  flow, 

For  the  case  u'  «  the  integral 

(4.58)  is  substantially  simplified 

and  for  b _ we  have  the  formula: 

cr 


Fig.  4.25.  Dependence  of  the 
critical  parameter  of  injection 
the  ratio  of  the  molecular 
weights  and  the  temperature 
factor. 
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Thus,  the  temperature  factor  affects 
bcr  in  approximately  the  same  manner 
as  it  does  the  coefficient  of  friction  during  streamlining  of  an 
impenetrable  plate. 

In  Fig.  24  are  given  the  results  of  the  numerical  solution  of 
equation  (4.57).  Along  the  axis  of  ordinates  is  placed  the  ratio  of 
y  (for  given  values  of  ^  and  )  to  ¥,  calculated  by  formula  (5.55). 
Along  the  abscissa  is  placed  the  ratio  of  b  to  b  ,  calculated  by 
equation  (4.58).  Curve  1  corresponds  to  formula  (4.8). 

As  can  be  seen,  the  exact  calculation  for  the  case  of  injection 
of  freon-12  in  a  flow  of  air  is  near  the  approximate  formula  (4.50). 

In  any  case  the  divergences  are  of  the  same  order  as  with  the 
experimental  data  in  Fig.  22.  Exact  calculation  for  the  case  of  helium 
injection  in  a  flow  of  air  confirms  satisfactoriness  of  the  approxi¬ 
mation  of  the  influence  of  the  temperature  factor  ty,  but  deviates 
more  significantly  from  approximation  of  the  influence  of  the  factor 

FIT 


cr 


From  the  data  of  Fig.  22  on  the  injection  helium  into  air,  it 
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Fig.  24.  Solution  of  equation  (4.57).  1) 

calculation  by  equation  (4.8);  2,  3)  by 
experimental  points. 

is  clear  that  this  divergence  decreases  upon  transition  to  a 
calculation  by  Rex. 

Considering  this  circumstance,  and  also  the  absence  of  experi¬ 
mental  data  for  a  nonisothermal  nonuniform  boundary  layer,  in  a 
first  approximation  it  is  possible  to  assume  the  approximate  relation 
(4.30). 
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CHAPTER  V 


STREAMLINING  OF  A  CURVED  SURFACE 


Definitions  of  Cyrillic  Items  In  Order  of  Appearance 
ct  =  wall  =  wall 
Kp  =  cr  =  critical 
T  =  t  =  thermal,  turbulent 


5.1  Limiting  Parameters  of  Breakaway  of  an  Isothermal 
boundary  Layer  on  an  impenetrable  Surface 

During  streamlining  of  a  curved  surface  the  flow  velocity  on  the 


external  boundary  of  a  dynamic  boundary  layer  changes  along  the  stream 
lined  contour  and,  consequently,  ^  0. 

Nozzle  flow  is  distinguished,  when  ^  <  0,  and  divergent,  when 


0. 


During  nozzle  flow,  the  stream  is  accelerated,  the  direction  of 
motion  of  the  liquid  coincides  with  the  direction  of  the  pressure  force, 
and  the  boundary  layer  is  always  stable  in  the  sense  that  it  is  not 
detached  from  the  streamlined  surface. 

During  diffusion  flow  the  stream  is  slowed  down,  pressure  increases, 
and  its  action  is  directed  toward  the  motion  of  the  liquid. 

The  pressure  gradient  in  the  boundary  layer  is  determined  by 
formula  (1.10),  i.e.,  a  change  of  pressure  occurs  in  strict  conformity 
with  a  change  of  veolcity  wQ. 
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Inside  the  dynamic  boundary  layer  flow  is  inhibited  by  friction 
and  wx  <  Wq.  Therefore,  the  reserve  of  the  kinetic  energy  of  flow  inside 
the  boundary  layer  is  insufficient  for  full  surmounting  of  the  action 
of  the  field  of  pressures  directed  toward  it.  As  a  result  a  positive 
pressure  gradient  evokes  inside  the  boundary  layer  inhibition,  and  then 
stopping  and  reverse  of  the  current  of  liquid  near  a  streamlined  body. 
This  phenomenon  is  called  separation  of  the  boundary  layer. 

Beyond  the  point  of  breakaway  a  vortex  motion  of  the  liquid  appears, 
accompanied  by  sharp  growth  of  the  resistance  of  pressure  and  still  not 

yielding  to  theoretical  calculation. 

Schematically  the  sequence  of 
deformations  of  the  velocity  profile 
in  the  region  of  divergent  flow  is 
shown  in  Fig.  25. 

Inasmuch  as  separation  of  the 
boundary  layer  is  characterized  by 


Fig.  25.  Diagram  of  deforma¬ 
tions  of  the  velocity  profile 
during  flow  of  a  boundary 
layer  with  a  pressure  gradient. 


conversion  of  the  flow  in  direct  proximity  to  the  wall,  i.e.,  there, 
where  the  flow  is  the  most  inhibited,  then  the  point  of  breakaway  is 
determined  by  the  condition 


(5.1) 


Correspondingly,  at  the  point  of  breakaway  cf  =  0. 

Actually,  separation  of  the  boundary  layer  does  not  occur  at  any 
exactly  fixed  point,  but  embraces  a  certain  finite  domain.  In  this 
region  the  velocity  profiles  are  deformed  from  the  mimimum  stable  form 
to  a  form,  corresponding  to  condition  (5.1).  In  addition,  pulsations 
are  possible,  accompanied  by  oscillations  of  the  point  of  breakaway 
within  the  limits  of  a  certain  region  Ax. 

Measurements,  conducted  by  many  researchers,  show  that  during 
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nozzle  flow  the  velocity  profiles  are  more  populated,  and  during  dif¬ 
fusion  flow  are  less  populated  than  the  velocity  profile  for  the  case 


2r  = 

dx 


0. 


For  the  point  of  breakaway  of  an  isothermal  boundary  layer  at  an 

impenetrable  wall  we  have  the  condition: 

{*/  ^0;  p  -  p.w  )  -  ft;  j,  -  o; .  ( 5 . 2 ) 

Putting  in  equation  (2.5)  the  value  of  the  considered  quantities,  we 
obtain  the  equation,  determining  the  velocity  profiles  at  the  point  of 
breakaway  of  an  isothermal  turbulent  boundary  layer  at  an  impenetrable 
wall ; 


(5.3) 


The  index  "cr"  indicates  that  the  corresponding  quantities  are 
referred  to  a  cross  section,  in  which  separation  of  the  boundary  layer 
appears . 

Distribution  of  the  tangential  stresses  across  the  boundary  layer 
at  an  impenetrable  wall  in  a  sufficiently  general  form  it  is  possible 
to  write: 

*  =  —-  =  *•(«)  +  A«?i  0).  (5.4) 


During  approximation  by  a  cubic  parabola  in  accoradance  with 
(2.35) 

fid  =  0-y.  I  (5.5) 


and,  taking  into  account  that  the  form  parameter  A  is  determined  by 
formula  (2.56)  and  at  the  point  of  breakaway  t  =  0,  we  have 


=  “V  *  /-M 

>*0  P*o  **  \  V/Uf 


(5.6) 


dw. 


C  VA¥1q 

where  f  =  —  '  is  the  form  parameter,  not  related  to  the  quantity 
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(5.7) 


Thus,  if  the  integrand  and  lower  limits  of  integration  in  equa¬ 
tions  (5.7)  and  (5.3)  are  known,  then  the  breakaway  velocity  profile 
in  the  turbulent  nucleus  of  a  boundary  layer  and  the  critical  values 

5** 

of  the  parameters  f,  — ,  and  H,  are  calculated,  i.e.,  all  the  quanti¬ 
ties,  characterizing  the  phenomenon  of  breakaway. 

Equations  (5.7)  and  (5.8)  can  be  solved,  if  we  accept  the  assump¬ 
tion  of  the  conservation  of  the  dependence  of  the  length  of  the  path 
of  mixing  l  on  the  form  parameter  f  and  we  approximate  tangential 
stresses  profile  by  a  cubic  parabola. 

Both  of  these  assumptions  always  can  be  considered  as  a  first 
approximation  of  the  real  distributive  laws  of  the  quantities  l  and  t 
in  the  critical  cro  s  section,  and  in  any  case  they  do  not  emerge  within 
the  framework  of  assumptions,  accepted  in  the  contemporary  semiempirical 
theories  of  a  turbulent  boundary  layer. 

Any  further  more  precise  definitions  of  these  dependencies  can  be 
easily  considered  by  the  corresponding  change  of  integrands  in  equations 
(5.7)  and  (5.8).  » 

Regarding  the  lower  limits  of  integration,  i.e.,  cr  and  , 

then  for  their  determination  in  the  general  case  knowledge  of  the  law 
of  stability  of  the  viscous  underlayer  is  necessary. 
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At  Re  -*  oo  -*  0,  od1  -*■  0,  and  the  critical  parameters  of  the  break¬ 
away  of  the  isothermal  turbulent  boundary  layer  at  an  impenetrable  wall 
are  determined  by  equations; 


f(-H 


(5.9) 


f(-H/ 


Taking  as  primary  experimental  fact  logarithmic  velocity  profile 
at  ^  =  0,  we  find  that 


(5.10) 


According  to  formula  (5.1)  h  =  0.4. 

According  to  the  other  experimental  data  the  value  of  this  constant 
lies  in  the  range  0.58  to  0.41,  i.e.,  very  close  to  the  quantity, 
defined  by  Nikuradze. 

Putting  relations  (5.5)  and  (5.10)  in  equation  (5.9) 1  we  obtain 
the  limiting  velocity  profile  in  the  region  of  the  breakaway  of  the 
turbulent  layer: 


Accordingly, 


_  In  (2/  4€«  +  2S  +45+0 
In  (2  / T  +  5) 


= 2H" ,n  (j  r' 6  +  s)F  ■ 


(5.11) 


(5.12) 


From  formula  (5.11)  it  follows  that  for  the  accepted  assumptions 
the  limiting  velocity  profile  in  the  cross  section  of  the  breakaway 
of  the  boundary  layer  does  not  depend  on  the  empirical  constants  of 
turbulence.  The  limiting  critical  value  of  the  form  parameter  f 
depends  on  the  constant  h. 
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At  H  =  0.4,  for  the  considered  consitionc  (T  =  const,  =  0, 

Re  oo) ,  we  have1 : 

(-/-£-)„ =ft062:  Pr  )„ - a30; 

(-^-)  =0,16;  ».,=  l,87;  A,  =  0,010. 

The  velocity  profile  (5.11)  is  approximated  by  the  exponential 
relation 

-  =  P.«  (5.14) 


(5.13) 


This  value  of  the  exponent  for  the  cross  section,  in  which  break¬ 
away  of  the  turbulent  layer  from  the  impenetrable  wall  occurs,  is  very 
close  to  the  value  n  „  =  0.5,  found  earlier  from  other  considerations 
in  the  works  of  oss,  Stratford,  and  Townsend.  However,  it  seems  to 
us  that  in  these  investigations  less  general  assumptions  are  accepted 
and  not  all  of  the  complex  of  parameters,  characterizing  the  region  of 
breakaway  are  obtained.  In  particular,  the  critical  values  of  the  form 
parameter  f  were  not  calculated. 

From  the  given  calculations  it  is  clear  that  just  the  form  param¬ 
eter  f  preserves  a  finite  critical  value,  even  at  Re  -*■  co,  whereas  the 
critical  value  of  the  Buri-Loytsyanskiy  form  parameter 


r_=-2C  (5.15) 

<h 

approaches  infinity  (since  at  Re  -+  oo,  c  -*•  0).  In  the  region  of  not 
very  large  values  of  Re  number  the  quantity  Tcr  has  fully  finite  values. 

Qualitatively,  this  conclusion  was  earlier  substantiated  in  the 
works  of  L.  E.  Kalikhman  and  G.  M.  Bam-Zelikovich. 

In  Fig.  26  is  given  a  comparison  of  the  velocity  profile,  calcu¬ 
lated  by  formula  (5.14),  with  results  of  the  measurements  of 


1Values  of  6*  and  6**  are  found  by  formulas  (1.24)  and  (1.25)  upon 
substitution  in  them  of  o>  from  (5.11)  and  p  =  pQ. 
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Fig.  26.  Comparison  of 
calculated  velocity  pro¬ 
file  at  the  point  of 
breakaway  of  the  turbu¬ 
lent  boundary  layer  wich 
experimental  data:  1) 
calculation  by  formula 
(5.11)l  2)  calculation 
according  to  the  1/7 
power  law;  O  —  experi¬ 
ments  of  I.  Nikuradzej 
•  —  experiments  of  A.  I. 
Leont’ev,  A.  N.  Oblivin, 
P.  N.  Romanenko. 


I.  Nikuradze  and  A.  I.  Leont’ev,  A.  N. 
Oblivin,  and  P.  N.  Romanenko. 

It  is  clear  that  the  theoretical 
limiting  velocity  distributions  in  the 
prebreakaway  region  of  flow  is  satisfac¬ 
torily  confirmed  by  the  experimental  data 
for  finite  Re  numbers. 

The  critical  value  of  the  form  param¬ 
eter  H  from  the  experiments  of  I.  Niku¬ 
radze  are  equal  to  1.8,  from  the  experi¬ 
ments  of  E.  Grushwitz  —  1.9,  i.e.,  very 
close  to  the  found  theoretical  value. 

According  to  some  other  experiments 
the  quantity  H  attains  2-2.5.  It  is 

L  X 

possible  that  this  increase  of  the  value 


Hcr  is  connected  with  the  roughness  of  the  streamlined  surface  or  some 
other  sort  of  factors,  not  considered  by  the  accepted  approximations 


of  the  functions  i(y)  and  T(y). 


5.2.  Condition  of  Stability  of  the  Viscous  Underlayer 
During  streamlining  of  an  impenetrable  plate  by  an  isothermal  flow 
with  ^  =  0,  the  condition  of  stability  of  the  viscous  underlayer  can 
be  written  in  the  form 


£#1  —  — ~  =  tJ0  =  const .  ( 5  •  i6 ) 

V 

In  the  two-layered  diagram  of  a  turbulent  boundary  layer  »  1,6 
and  Re^  =  154  [-  ?  -  original  illegible].  During  gradient  flow  the 
condition  of  stability  is  changed  due  to  a  change  of  the  configuration 
of  the  velocity  profile  in  the  viscous  underlayer.  The  character  of 
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this  change  can  be  clarified  by  proceeding  from  the  following  consid¬ 
erations  . 

Let  us  assume,  that  formula  (5.14)  approximately  describes  the 
velocity  profile  in  the  cross  section  of  breakaway  and  at  finite  Re 
numbers.  Such  an  assumption  is  more  valid  the  larger  the  Re  number. 
Furthermore,  it  is  also  confirmed  by  the  experimental  data  (see 
Fig.  26). 


In  the  viscous  underlayer  the  velocity  distribution  is  determined 
by  formula  (2.21),  which  at  cf  =  0  takes  the  form 

i 


Q\  —  — 


Re 


a* 


(5.17) 


Intersection  of  profiles  (5.14)  and  (5.17)  at  the  values  of  the 

critical  parameters  from  (5.13)  will  give  the  values  of  5  and  ou 

l  c  r  1 


cr 


to  a  first  approximation. 
We  have 


e  -  2.84 

"■ • 

t,S7 

— s:..  • 


(5.13) 


Hence, 


•* Re  *-  ’"i'.irt 81  “ 


28. 


In  addition  the  sum  of  the  exponents  for  Re**  in  formulas  (5.18) 
is  rounded  off  to  1  with  accuracy  up  to  0.1,  which  is  fully  permissible, 
taking  into  account  the  approximate  character  of  these  dependencies. 

Thus,  the  Reynolds  number  for  the  viscous  underlayer  essentially 

decreases  with  an  increase  of  divergence. 

Obviously,  such  a  decrease  of  stability  of  the  viscous  underlayer 
is  connected  with  a  corresponding  distortion  of  the  velocity  profile 
with  an  increase  of  the  pressure  gradient.  In  some  measure  this  cir¬ 
cumstance  is  taken  into  account  in  the  Reynolds  number,  constructed 
from  the  velocity  derivative  on  the  external  boundary  of  the  viscous 
underlayer: 


■I 


I 


(5.19) 


•  Q 

At  (dp/dx)  =  0  Re1  is  identically  equal  to  R«1  or  t|1q.  Assuming, 
as  in  all  the  preceding  calculations,  the  quantity  t)1q  =  11.6,  we 
find  that  at  f  =  0  Re.^  =  134. 


At  the  point  of  breakaway 


(5.20) 


Introducing  here  the  value  of  the  critical  parameters  from  (5.15) 
and  (5.18)  and  rounding  off  to  two  the  exponent  for  Rft  number  in 
the  cube  of  the  quantity  ,  we  find  that 

0.161  “5/- 

Thus,  Re1  actually  is  stabler  than  Re1. 

The  condition 

to,-  const  (5.21) 

is  identical  to  the  Stseblevskiy  condition  introduced  earlier 

(Hr/ 7),.,.  (5.22) 


5.3. 

Above  was  shown  the  method  of  calculation  of  the  values  of  the 
form  parameters  at  the  point  of  breakaway.  The  value  of  the  form 
parameters  at  (dp/dx)  =  0  is  calculated  from  the  experimental  univer¬ 
sal  velocity  profile  <p(*n ) . 

According  to  the  degree  of  the  population  of  the  velocity  pro¬ 
file  that  occurs  with  an  increase  of  the  form  parameter  f,  the  magni¬ 
tude  of  the  exponent  n  decreases.  Assuming  that  at  f  -*•  co  n  -*•  ao 
from  (5.15)j  we  find  that  in  this  case  H-+1. 

Thus,  we  have  the  conditions: 

/=/- ri 

/  —  0;  //=  H0\  (5.25) 

/=  H=  1. 


troximation  of  the  Relation  Between  the  Form 
Parameters  H  and  f 
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In  addition  2  >  H  >  1,  i.e.,  changes  minutely  in  the  segment 


from  f  to  f  =  co. 
cr 

All  these  conditions  are  well  satisfied 
polation  formula: 

h-\  T 

M*  —  I  ~~\  —  l  /  * 


by  the  simple  inter- 

(5.24) 


where  7  =  (f/fcr). 

To  the  logarithmic  velocity  profile  corresponds  the  value 


Hu  *=  1  r 


\  2fj, 


(5.25) 


'-V  iff,, 

Consequently, 

//=//</;  (5.26) 

In  Pig.  27  is  given  a  comparison  of  calculations  by  formula  (5.24) 
with  the  experimental  data  of  I.  Nikuradze,  plotted  in  the  coordinates 
(H;  f),  where  H  =  (H/Hcr). 


Fig.  27.  Dependence  of  the  pa¬ 
rameter  H  on  the  parameter  f: 

1)  by  formula  (5.24);  2)  cal¬ 
culation  for  a  laminar  boundary 
layer;  points  from  the  experi¬ 
ments  of  Nikuradze. 


We  can  see  the  satisfactory  agreement  of  calculation  with  experi¬ 


ment  . 
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5 . 4 .  Law  of  Resistance  for  an  Isothermal  Boundary  Layer 
on  an  Impenetrable  Surface  at 
-  (dp7dxTV  0 - 

Fiom  equations  (2.6)  and  (5.4)  it  follows  that  the  velocity  pro¬ 
file  in  the  turbulent  nucleus  of  an  isothermal  boundary  layer  on  an 
impenetrable  curved  surface  is  determined  by  the  equation: 


=••+] 


r 


T'» 


(5) 


» ri  (*)  d 5. 


(5.27) 


where,  as  before,  ¥  =  (cf/cfO^R#*# * 

Distribution  of  velocities  in  the  viscous  underlayer  is  described 
by  equation  (2.21). 

The  condition  of  stability  we  will  write  by  means  of  the  quantity 


Re^. 

k  = 


Taking  into  account  equation  (2.19)  for  the  case  T  =  const  and 
0,  we  have 


y 

2 


(5.28) 


At  the  point  of  breakaway  ¥  =  0  and  from  (5.28),  (5.19)  follows. 

The  system  of  equations  (5.27),  (2.21),  and  (5.28)  determines 
the  quantity  ¥,  if  the  function  £■(£);  <Pq(£);  ^(S);  ft«^(f)  and  cf.0(Re**) 
are  known. 

Actually,  by  assuming  in  (5.27)  o>  =  1  and  correspondingly  the 
upper  limit  of  the  right  integral  is  also  equal  to  1,  we  have  the 
system  of  equations: 


i 

[ 


r 

V  * 


-f-  *.<« 


I** 


n?.<5>  </;; 


(r-f-  )(-£-  *•")'  =  *».• 


(5.29) 


When  the  indicated  functions  are  known,  then  the  system  of  equa¬ 
tions  (5.29)  is  solved  in  a  first  approximation,  if  we  introduce  in 
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value  5  /6,  defined  via  the  quantity  H  from  the  formulas  (3.13),  and 
the  quantity  H,  in  turn,  is  determined  by  approximation  (5.24). 

Further,  by  introducing  in  equation  (5.27)  the  thus  found  values 
of  and  we  will  calculate  the  quantity  5  /5  in  a  3econd 

approximation  and  again  solve  the  system  of  equations  (5.29).  As  a 
result  in  the  second  approximation  are  calculated  not  only  f;  and 
cd1,  but  also  6  /5  and  H. 

The  authors  together  with  N.  N.  Kirillova  and  G.  P.  Zykin  solved 
the  system  of  equations  (5.28)  in  a  first  approximation.  In  addition 
the  quantity  1/6  was  determined  by  formula  (5*10)  at  h  =  0.4;  the 
functions  9q(£)  and  q>^(£)  —  by  formulas  (5.5)j  the  quantity  c -  by 
formula  (3*7) j  the  quantity  H  —  by  formula  (5.24);  the  quantity  H 

v  X 

was  taken  equal  to  I.87. 

Thus,  in  the  semiempirical  theroy  of  a  turbulent  boundary  layer 
developed  here,  as  a  primary  experimental  fact  is  taken  a  definite 
universal  velocity  profile  <p( rj ) ,  and  not  some  special  hypothesis  on 
the  length  of  the  path  of  mixing.1 

The  full  system  of  equations  of  the  examined  first  approximation 
takes  the  form: 


(5.30) 


xIn  the  majority  of  works  Z  =  >ty  (Prandtl  formula)  or  l  ■ 

2  2 

=  n[(dwx/dy)/(d  wx/dy  )]  (Karman  formula)  is  assumed.  "Constants  of 
turbulence"  h  and  rj^0  are  determined  from  the  same  experimental  velo¬ 


city  profile.  Experiment  shows  that  these  dependences  are  sufficient!,/ 
accurate  only  in  a  narrow  prewall  region.  Usually  t  =  t  ^  *  const 

is  assumed,  i.e.,  a  condition  valid  also  in  direct  proximity  to  the 
wall. 


The  hypothesis  on  the  conservation  of  the  function  l( y),  naturally, 
is  also  necessary  in  these  theories. 
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The  critical  parameters  are  determined  preliminarily  from  this 
system  of  equations  at  ¥  =  0  and  f  =  1. 

Calculations  were  carried  out  on  the  electronic  computer  of  the 
Calculating  Center  of  the  Siberian  Section  of  the  Academy  of  Sciences 
of  the  USSR.  The  results  are  given  in  Table  5.1  and  in  Figs,  28,  29, 
and  30. 


Table  5.1.  Values  of  the  Parameters  of  an 
Isothermal  Boundary  Layer  at  the  Point  of 
Breakaway  from  a  First  Approximation _ 


MU' 

i*to* 

s-w 

1  •!•*'• 

5.  JO' 

i.|» 

•»or 

0,133 

0,077 

0,0416 

0,0354 

0.0294 

0,0183 

*ior 

jo.0103 

0,00625 

0,00215 

0,00135 

0,00046 

0,00029 

(t)„ 

0,16 

O.lti 

0,16 

0,16 

0.16 

0.16 

""/or 

0.010 

0,010 

0,010 

0,0100 

0,0010 

0.0010 

—/or 

4.7 

6.7 

9,0 

10,0 

12,8 

14,1 

##or 

m  * 

1.87 

1,87 

1,87 

1,87 

1,87 

1,87 

As  can  be  seen,  the  thickness  of  the  viscous  underlayer  increases 
in  the  region  of  diffusion  flow,  and  the  velocity  on  the  external 
boundary  of  the  viscous  underlayer  decreases.  Besides  the  larger  the 
Re  number,  then  the  relatively  stronger  are  these  changes. 
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IThe  relative  change  of  the  coefficient  of  friction  in  the  dif¬ 
fusion  region  is  very  significant  and  depends  both  on  the  form  param¬ 
eter  and  Reynolds  number,  i.e., 

*=**(/;  Re**).  (5.3 

In  addition,  the  greater  the  value  of  Re  ,  the  more  rapidly  the 

coefficient  of  friction  decreases  with  the  growth  of  the  parameter  f. 
Calculations  were  performed  for  two  conditions  of  stability: 


V/Re,  -  T^0  *  11.6  and  Vie ^  *=  11.6  -  4. Of.  The  latter  is  a  linear 
interpolation  between  the  value  of  the  quantity  iJL-t/j-]  , 

determined  experimentally  for  the  case  f  =  0,  and  the  value  of  this 
quantity  at  f  =  1,  calculated  from  the  limiting  velocity  profile. 

Calculations  by  both  conditions  of  stability  practically  coincide, 
i.e.,  possible  inaccuracy  in  approximation  of  the  function  Re1(f)  is 
immaterial  for  calculation  of  the  parameters  of  a  turbulent  boundary 
layer. 


5.5.  Solution  of  Equation  of  Momentum  for  an  Isothermal 
Boundary  Layer  on  an  Impenetrable  Surface 

It  is  convenient  to  write  the  momentum  equation  in  the  form: 

-!-/„«»*  (I  H.,  (5.J2) 

_  ##  ww 

Here  Re  =  (wQ6*  /v)  is  the  flow  value  of  the  Reynolds 

number,  constructed  according  to 
the  thickness  of  the  momentum  loss; 

ReL  =  (w0I/v)  is  the  Reynolds  number,  constructed 
h  according  to  the  local  value  of  the 

velocity  wQ  and  the  characteristic 

linear  dimension  of  a  streamlined 
surface  L; 

X  =  (x/L)  is  the  relative  distance  from  the 
inlet  edge; 

f  ;  H  are  the  values  of  the  form  param¬ 
eters  f  and  H  at  the  point  of 
breakaway  of  the  boundary  layer; 
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f  =  (f/f  );  H  =  (H/H  )  are  relative  values  of  the  form 

parameters; 

is  relative  change  of  coefficient 
of  friction  under  the  influence  of 
a  pressure  gradient. 

H  and  are  complex  functions  of  7  and  R«*#.  Therefore,  an  exact 
solution  of  the  momentum  equation  can  be  obtained  only  by  numerical 
methods  for  a  given  velocity  distribution  outside  the  boundary  layer. 
The  solution  is  radically  simplified  by  linearization  of  equation 


(5.31). 


We  will  represent  this  equation  in  the  form: 


2 

€/m  k»L 


***** 


=  F{f). 


where 


At  ?  =  0  we  have: 


at  f  =  1  we  have: 


Fif)  --**,-(1 4 


F{f)  |; 


(5.33) 


(5.3*1 ) 


*V=0;  h-u  /=■{/»  =  i  iit/gr.. 

Linear  Interpolation  gives  the  formula: 

r</>=i  (i  +  H.jfrv  (5.35) 

Calculations  by  the  approximate  formula  (5.35)  coincide  very 
closely  with  calculations  by  the  exact  formula  (5.3*0  upon  substitution 
of  the  values  of  H  and  7,  which  are  determined  by  formula  (5.25)  and 
the  graph  in  Fig.  28,  at  |fcr|  *  0.01. 

Introducing  in  (5.33)  the  value  of  F(f)  from  formula  (5.35)  and 
cf0  from  formula  (5.14),  we  obtain  the  equation: 

— +(!+«.,)  —  -37  -  —  '  -g£m~  (5.36) 


Integrating  this  equation,  we  find  that: 


X 

He**  =  u'o  ''  p  "  0 Re0  f  Wo  dX  + 


(5.37) 
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where 


*i  —  1  +  Htp\  *j  —  I  -|-  (1  4-  m)  (l  -f  »„  — 


—  _•» 


Re  u  = 


IPoi  l* 


wQ1  is  the  velocity  of  the  current  incident  on  the  leading  edge  of  the 
streamlined  body;  is  a  dimensionless  coordinate  of  the  beginning 
of  the  turbulent  boundary  layer. 

By  equation  (5.37)  the  Re  is  distributed  along  the  contour 
of  the  body,  covered  by  the  turbulent  boundary  layer.  In  addition, 
the  flow  velocity  outside  the  boundary  layer  wQ(X)  should  be  given. 

The  local  value  of  the  form  parameter  is  determined  by  formula 


t-' 


to**  dw„ 


**»«.  dX 
*# 


(5.38) 


For  given  values  of  f  and  Re  the  value  of  is  given  by  the 

graph  in  Fig.  28,  and  by  formula 

*/  =  *'/•.  (5.39) 

the  local  value  of  coefficient  of  friction  is  calculated. 

The  values  of  the  coefficients  B  and  m  during  solution  of  the 

momentum  equation  are  selected  for  a  given  interval  of  values  of 

•ft# 

Re  numbers  in  accordance  with  Table  3.1. 

The  quantity  cf0  in  formula  (5.39)  can  be  calculated  by  formula 
(5.8). 

For  not  very  large  values  of  the  form  parameter  f,  it  is  possible 
to  put  H  =  1.3  and  if  =  1  in  the  momentum  equation.  In  this  case 


Re* 


1  +  m 


B  Re, 


•P"  1 

x, 


3  m 


dX  r 


(5.40) 


{Re**  w 


i,Jj*  ji  -i 


#*  4 

For  Re  numbers  <  10  it  is  possible  to  assume  m  =  0.25  and 


B  ■  0.0258.  Then 


^•  =  WM[»M)16l«.p«dX  (5.41) 
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In  the  works  of  Gruschwitz,  Dengot,  and  Tetervin  it  is  shown 
that  the  form  parameter  H  may  be  a  good  criterion  for  determination 
of  the  point  of  breakaway  of  the  turbulent  boundary  layer.  Besides, 
as  was  noted  above,  the  critical  value  of  the  form  parameter  H,  found 
experimentally  by  Gruschwitz,  practically  coincides  with  those  cal¬ 
culated  by  us.  In  the  opinion  of  the  above  researchers,  distribution 
of  the  form  parameter  H  along  the  contour  of  a  streamlined  body 
depends  not  only  on  the  local  value  of  the  Re  number,  but  alBo  on 
the  distributive  law  of  pressure,  i.e.,  on  function  wQ[H].  In  this 
connection  empirical  relations  of  the  type, 

4**  d  H  ~ 

-r-  •  -7JT  = 

are  introduced. 

The  degree  of  reliability  of  these  relations  is  vague,  inasmuch 
as  the  corresponding  experimental  data  are  not  published  in  sufficient 
volume . 

The  presence  of  some  sort  of  influence  of  the  character  of  the 
function  wq(H)  on  the  local  parameters  of  the  boundary  layer  does  not 
contradict  the  general  considerations  about  similarity  and,  within 
the  framework  of  our  theory,  should  be  taken  into  consideration  by 
introduction  of  the  corresponding  correction  in  the  approximation  of 
the  tangential  stress  profile. 

Regarding  the  quantities  f  ,  Hcr,  and  f(f;  Re  ),  the  influence 
on  them  of  the  indicated  factor  has  not  yet  been  observed. 

5 .6 .  Law  of  Heat  Transfer  in  the  Diffusion  Region 
of  a  Quasiisothermal  Boundary 
layer  on  an  Impenetrable 
- Waff - 

The  presence  of  a  pressure  gradient  essentially  disturbs  the 
similarity  between  the  processes  of  friction  and  heat  transfer  in  the 
boundary  layer.  In  addition  the  process  of  heat  transfer  possesses 
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a  significant  degree  of  conservation  which  is  revealed  already  upon 
comparison  of  the  distributive  laws  of  the  tangential  stresses  and 
heat  flux  along  the  cross  section  of  the  boundary  layer. 

By  formulas  (2.35)  and  (2.74)  for  an  impenetrable  surface  (b  =  0): 


*  =  l  +  A  5  (2A  ~  3) -  (A  -j-  2) 

f  =  1  -Sfc’  +  Kj. 


(5.42) 


Within  the  limits  of  the  given  approximation  the  distribution  of 
the  heat  flux,  in  general,  does  not  depend  on  the  pressure  gradient, 
whereas  the  distribution  of  tangential  pressures  changes  substantially 
with  a  change  of  the  quantity  (dp/dx). 

In  Fig.  31  is  shown  the  distribution  of  tangential  stresses  and 
heat  flux  along  the  cross  section  of  the  boundary  layer  at  the  point 

of  breakaway,  calculated  by  formulas 
(5.42).  In  the  wall  region  the  path  of 
the  curves  q(£)  and  ^t/PqWq^  ( £ )  ls  Qul^e 
different. 

Let  us  carry  out  a  calculation  of 
the  intensity  of  heat  transfer  at  the 
point  of  breakaway  at  Pr  *  1  and  6^.  <  6, 
In  this  case  *  l  and  equation  (2,64) 
for  quasiisothermal  conditions  will  take 
the  form 

’M-H'-tt  --3T  (5-4?) 


\ 

j  M  *  l\ 

i  1  1  \ 

v] 

m 

mmm 

SSI 

o  Q.I  a«  Q*  ? 

Fig.  31.  Distribution  of 
tangential  stresses  1)  and 
heat  flux  2)  at  the  point 
of  breakaway  of  the  boun¬ 
dary  layer  by  formulas 
(5.42). 


or,  taking  into  account  (5.10)> 


SI  as  0,16;* 


d  u 


dh 


d\  di 

Taking  the  velocity  distribution  from  formula  (5.14),  we  have 

1 JL,. 

J  M  ' 


(5.44) 


(5.45) 


«  «  0.0685  — 


u. 

J  t«.« 


=  0.0295  r, 


0.43 


-K-r 
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Here 


».  si  Pr  Re"  ~  (, . 


Disregarding  the  quantity  [ (6/6t )5± ]0>^5  as  compared  to  1  and  assuming 
Pr  =  1,  we  obtain  the  value: 


Stuf^ 


0,005 


I  +  0,0395 


(5.^6) 


Putting  in  (5.46)  the  value  of  (5  /6)  =  0.16  and  from  formula 

(5.18),  we  have 

AOtf 

(5.47) 


Stlf^ 


*r**°r(l  +0,71  1 


Assuming  St  =  0.0129  Re**“0,2^’,  we  find  that  in  the  considered 


conditions 


st>, 

St, 


3.5 


’  (5.48) 

"R 

By  formula  (5.48),  in  the  region  of  Re  numbers  *  (5-10) *10 
the  ratio  fltCI/flt0  is  equal,  on  the  average,  to  1. 

With  growth  of  the  Re  number  the  critical  value  of  the  St  number 


becomes  less  than  St. 

Thus,  a  theoretical  appraisal  shows  that  for  the  practically 
most  frequently  met  values  of  Re  ,  the  law  of  heat  transfer  almost 
does  not  change  with  a  change  of  the  pressure  gradient.  This 
important  conclusion  is  well  confirmed  by  the  experimental  data  shown 
in  Figs.  52,  55,  and  54.  It  is  distinctly  clear  that  upon  a  signifi¬ 
cant  decrease  of  the  coefficient  of  friction  and  an  abrupt  deformation 
of  the  velocity  profile  with  an  increase  of  the  diffusion  nature  of 
flow,  the  values  of  the  criterion  St  and  the  temperature  profile 

almost  do  not  change. 

Nonetheless,  the  theory  shows  that  at  Re  oo  the  value  of  St, 
a’ though  slowly,  approaches  zero.  This  tendency  also  is  revealed 

in  experiments. 
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-0.0282;  O  -  f  Re**0,25 
-0.0740;  (D  —  f  Re**0*25 
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5.7.  Solution  of  the  Energy  Equation  for  a  Quasiisothermal 
Boundary  Layer  on  an  Impenetrable  Surface 


The  equation  of  the  energy  of  a  boundary  layer  (1.26)  during 
constancy  of  the  physical  properties  of  the  flow  and  an  impenetrable 


wall  can  be  written  in  the  •‘orm 


4  Rdf  i 

dX  'AT 


dX  T 
dX 


=  St  •  ReL 


Here  X  =  (x/L)  is  the  relative  length; 


Nd 

XPr 


(5.49) 


Re,  =  (w0L/v)  is  the  Reynolds  number,  calculated  from  the  local 
L  value  of  the  velocity  wQ  and  from  the  character¬ 

istic  linear  dimension  of  the  body  —  L; 


St  =  (o./cyv!{  )  is  the  local  v.viue  of  the  Stanton  number; 

Nu  =  (ax/A)  is  the  local  value  of  the  Nusselt  number. 

Due  to  the  conservation  character  of  the  law  of  heat  transfer 


relative  to  the  pressure  gradient, in  the  majority  of  practical  cal- 

**  il 

culations  (at  Re  <  10  ),  it  is  possible  to  assume  the  value  of  the 
St  number  equal  to  its  magnitude  at  (dp/dx)  =  0.  Consequently,  for 
quasiisothermal  flow  formula  (5.55)  can  be  used.  Putting  this  value 
df  St  in  (5.49),  we  obtain  a  linear  equation  relative  to  the  quantity 
Re**l-Hn.  carrying  out  the  calculation,  we  find  that 


»».'■  -  |-~5r  |'  n'  ~“X  ■: 

X, 

-i  r)i* “l57*. 


(5.50) 


where  Ren  =  (wni  I/v)  is  the  Reynolds  number,  calculated  on  the  velocity 
0  U1  of  the  approaching  stream  Wq^  and  the  charac¬ 

teristic  length  of  the  body  L; 

w0  =  (Wfv/wQ.)  is  the  relative  velocity  outside  the  boundary 
0  layer. 

Putting  the  value  of  1#^*,  determined  from  formula  (5.50),  in 
formula  (5.55),  we  find  the  local  value  of  the  Stanton  number  and 
correspondingly  the  heat  transfer  coefficient.  In  addition,  the 


functions  Wq(x)  and  AT(x)  should  be  given  beforehand. 


5.8.  Influence  of  Nonlsothermalness  on  the  Parameters  of 
Breakaway  of  the  Boundary  Layer  from 
an  Impenetrable  Surface 


The  velocity  distribution  in  the  region  of  the  breakaway  of  a 


nonisothermal  layer  on  an  impenetrable  wall  is  determined  by  equation 


(5.51) 


The  solution  of  the  right  integral  is  expressed  by  formula  (5.8). 
Prom  here  we  obtain,  that  in  the  limiting  case,  when 

st  -*  o,  *»,  -*■  o,  ?  -*>  o, 

t 


(5.52) 


At  p  -  pQ  we  obtain  formula  (5.9)  for  the  limiting  case  (“i  =  0) 
Thus,  the  relationship  of  the  limiting  critical  values  of  the 
form  parameter  for  nonisothermal  and  Isothermal  flows  Is  determined 
by  integral 


(5.53) 


The  value  of  this  integral  for  a  constant  coefficient  of  non¬ 
similarity  of  the  temperature  and  velocity  fields  was  calculated  in 
Chapter  III.  For  gradient  flow  in  the  general  case  e  =  e(£).  Thus, 
during  isothermal  flow  at  the  point  of  breakaway  the  velocity  profile 

is  determined  approximately  by  formula  (5.14).  At  the  same  time  In 

4 

the  region  of  Re  numbers,  at  least  up  to  10  ,  the  law  of  heat  trans 
fer  almost  does  not  change  with  a  change  of  the  form  parameter,  i.e., 
Ilf.  m  n^Q  w  (1/7)  .  For  these  conditions  from  formula  (2.57) s 
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.  =  (5.54) 

where  Eq  Is  the  value  of  the  coefficient  of  nonsimilarity  of  the 
temperature  and  velocity  fields  at  (dp/dx)  -  0. 

But  at  Re  -*>  oo  St  -*■  0,  i.e.,  the  processes  of  friction  and  heat 
transfer  in  the  diffusion  region  become  similar.  Therefore,  inasmuch 
as  integral  (5.53)  expresses  the  relationship  of  the  limiting  critical 
values  of  the  form  parameter,  its  solution  coincides  with  formula 
(5.48). 

This  result  will  also  apply  for  finite  Rt  numbers.  Actually, 
even  for  moderate  values  of  Re  ,  from  formula  (5.54)  the  quantity 
e  >  0.5£q  already  at  £  =0.1j  at  £=0.5  e=  0.82e0j  at  £  ■  1  e  =  eQ. 
At  the  same  time  the  same  quantity  e,  as  was  shown  in  Chapter  III, 
weakly  affects  the  considered  integral. 

From  the  above  it  follows  that: 

a)  during  moderate  flow  velocities  of  a  gas 


(5.55) 


b)  during  significant  flow  velocities  of  a  gas1 

*  — — - — -fire  ... - +  ^ 

i!r  —II 


1  * 

f\ 

% 

/_L_ 

f\ 

(  j 

/wjn 

!)<*•  + A*) +  (A*)* 


—  arc  sin 


V. 

I  W  »)<**  + A*) +  <*+)•  J 


(5.56) 


Taking  into  account  the  value  of  ncr,  obtained  in  (5.14),  we  find 
that  the  limiting  velocity  field  at  the  point  of  breakaway  is  deter¬ 
mined  by  equation 

(5.57) 


jl/>d"-UjV  -Td- 


1The  formulas  (5.55)  and  (5.56)  were  obtained  previously  by 
L.  E.  Kalikhman,  but  not  as  limiting. 
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\ 


For  subsonic  flows,  when  it  *  lj 


-(i  *  -Oft41!* 


(5.58) 


For  flows  with  high  velocities,  when  it  is  noticeably  larger 
than  unity, 


arc  tin- 


n;*  -iv±  a; 


-  ire  >ln  -7-)  - 


(5.59) 


where 


£=1  41^*- iH^-f  a^)  +  (a;)»  . 

In  Fig.  35  is  shown  the  influence  of  the  temperature  factor  on 
the  limiting  velocity  profile  at  the  point  of  breakaway  of  a  boundary 

layer  of  gas  during  subsonic  flow.  In  the 


Fig.  35.  Influence 
of  heat  transfer  on 
the  limiting  velo¬ 
city  profile  at  the 
point  of  breakaway: 

1)  it*  m  1,  f  «=  0.25; 

2)  it*  *  1,  it  =  1.0; 

5)  /  -  1,  V  -  2.0; 

4)  »  6,  ii  «  6. 


same  place  is  given  a  graph  for  the  adiabatic 
supersonic  flow  of  gas  at  it  *  6  which  corres¬ 
ponds  to  M  =  5  at  r  =  1  and  M  «  5.?  at  r  «  0.9, 
The  temperature  factor  comparatively 
weakly  deforms  the  velocity  profile  at  the 
point  of  breakaway  of  a  boundary  layer  of  gas 
from  an  impenetrable  wall. 

In  Figs.  36,  37 >  and  38  are  given  the 
values  of  the  limiting  critical  parameters 

f  ;  Hcr;  (B/6**)cr  in  dependence  on  the 

* 

temperature  factor  it  at  subsonic  flows  {it  m  1 ) . 

The  quantity  H  almost  linearly  changes  with  the  growth  of  the 

v  A 

temperature  factor^.  Cooling  the  surface  {it  <  i)  increases  the  sta¬ 
bility  of  diffusion  flow. 

During  heating  of  a  streamlined  surface  {it  >  1)  the  stability 
of  diffusion  flow  decreases.  In  Fig.  39  is  shown  the  change  of  the 
quantity  of  the  critical  value  of  the  form  parameter  fcr  during 
supersonic  streamlining  of  a  body. 
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t 


l 


As  can  be  seen,  during  high  flow  velocities  the  region  of 
existence  of  a  stable  turbulent  boundary  layer  at  (dp/dx)  >  0  is 
extremely  limited. 


on  V'  and  on  ^  and 


In  Fig.  40  is  given  a  graph  for  the  value  of  Hcr  during  super¬ 
sonic  flow  of  a  gas.1 


5.9.  Solution  of  Equations  of  Energy  and  Momentum  for  a 
Nonlsothermal  Boundary  Layer  on  an 
Impenetrable  Surface  burlng 
Subsonic  iFlow  of  a  Gas 


The  Integral  energy  equation  for  the  region  of  finite  values  of 


** 


Bet  numbers.  In  which  one  may  assume  that  the  St  number  little  depends 
on  the  form  parameter  f,  will  take  the  form 

«/  Re'  Re**  d*T 


dX 


1  T 


d  X 


=  ReL  4  tS/„ 


(5.60) 


where 


_ 2 _ 

+  1.; 


(5.61) 


is  the  correction  for  nonisothermalness  taking  into  account  the 
finiteness  of  the  Ra  number. 


*A  more  conservative  parameter  (f*H)  which  is  explained  by  mutual 

V  X 

compensation  in  the  change  of  and  Hcr. 
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Formula  (5.61)  is  obtained  from  formula  (3.58)  at  V'  m  t  m  1. 

The  limiting  value  of  the  function  during  subsonic  flow  of  a 
gas  is  determined  by  formula  (3.55),  which  coincides  with  the  corres¬ 
ponding  limiting  formula  (3.47)  for  the  ratio  of  the  coefficients  of 
friction. 

In  practice  in  calculations  it  is  possible  to  be  limited  by  the 
limiting  formulas  for  or  to  introduce  this  function,  referred  to 

the  mean  value  of  c^q  for  a  given  streamlining.  Then  *»  Y(V')*  i.e., 

*# 

it  can  be  considered  independent  of  Re  .  The  solution  of  equation 
(5.60)  under  these  conditions  has  the  same  form  as  equation  (5.50)* 
but  under  the  integral  sign  stands  the  product  T^WqAT  .  The  value 
of  the  local  Stanton  number  is  determined  for  the  obtained  value  of 
the  Re^  number  by  formula 

St  =  *1'T  S/„,  (5.62) 

where  st0  is  determined  by  formula  (3.35). 

During  linear  interpolation  of  the  function  F(f)  its  value  during 
nonisothermal  flow  will  be  determined  by  formula: 

F(/l  af,-1 1  (5.6  J) 

For  the  coefficient  of  friction  we  will  obtain  the  expression 

(5-64) 

where  is  the  correction  for  the  Influence  of  the  pressure  gradient 
during  isothermal  flow. 

As  was  explained  in  Chapter  III,  the  function  for  the  coef¬ 
ficient  of  friction  differs  from  the  analogous  function  for  the 
coefficient  of  heat  transfer  by  the  form  of  calculation  of  the  coef¬ 
ficient  of  nonsimilarity  of.  the  velocity  and  temperature  fields  e.  It 
was  also  shown  there,  that  for  the  most  frequently  met  ration  0.5  < 

<  (6  /6^  )  <  2,  in  practical  calculations  the  influence  of  the  coef¬ 

ficient  e  on  can  be  ignored  and  e  =  1  can  be  assumed. 
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Solution  of  the  momentum  equation  takes  the  form 


where 


Trr)[-f  Re-  K 

xmpm^+oF^ 

j=$  (t+«.r)-^— . 


(5.65) 


(5.66) 


Here  Hcr  is  a  function  of  the  temperature  factor  V'. 

During  a  constant  wall  temperature,  the  integral  (5.65)  is 
simplified  and  can  be  written  in  the  form 


where 


«  =  I  .  H. 


up- 


(5.67) 


In  that  case,  when  it  is  necessary  to  determine  the  distribution  of 
the  wall  temperature  at  a  given  law  of  feed  (discharge)  of  heat 
qwaii(x)i  the  energy  equation  is  written  in  the  form 

_  9"ix)  (5.68) 


Hence 


R*7  =  — ft— :  ■  r  Re,:Xt. 
T 


(5.69) 


Taking  Into  account  formulas  (5.62)  and  (3.55).  we  find  that  at 
Jta.  Y  m  0  the  local  difference  of  temperatures  of  the  wall  and  the 

flow 


V 

frcr  L 


t 


I  *  m 


(5.70) 
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Putting  in  (5.70)  B  =  0.0258,  m  ■  0.25,  and  from  the  limiting 
formula  (3.55),  we  have 


Here,  as  before. 


H1"  (f^rr 


1 >=-y-.  i.  = 


(3.71) 


5.10.  Resolution  of  the  Equations  of  Energy  and  Momentum 

During  Supersonic  Flow  of  das 

* 

At  V'  >1  the  energy  equation  preserves  the  form  of  equation 


(5.60),  if  we  assume: 


A  T=  7-;t  -  rtt; 


/?*;/  = 


>*,*•* 


h,  W„L 


(5.72) 


!•<*> 


where  m-qq  lfl  the  viscosity  of  the  gas  at  the  temperature  of  stagnation 
outside  the  boundary  layer. 

The  integral  of  this  equation  has  the  forms 


*'■  =  “a  f  [-W*-  »*»-  f  i-T-fx 

a. 

—  '  1_i_ 

xt/ii  -  £/-)  a  r1  "<**  -(/}#;' a  r 


(5.7?) 


Here 


Re»>  - 


.  U  = 


m, 


•  1  ~~/rrr*T'' 


The  function  Tt  is  calculated  from  the  respective  formulae  of 
Chapter  III  for  supersonic  flow.  The  limiting  value  of  f  at  e  *  1 
is  determined  by  formula  (3.58). 

The  Stanton  number  will  be  determined  by  the  formula 

St  —  _ iLli—  / _ _*• _ \m 
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•It 


«* 


wh ete  the  value  of  Rt^  is  found  from  equation  (5.73). 

Taking  into  account  the  compressibility,  the  integral  of  the 
momentum  equation  will  take  the  form 


*# 


•• 


=  “*i  ■  T^r.ilf  Rt-) l^r)" x 


»  i 


where 


X  U\  1  —  U\  exp(/)</ A 


j*".. 


(5.74) 


At  Twall  "  const 

*,<• u  .[-Li-. u *«„ {  v, (Ju-f  >■ 

•^4 

—  j_JL_ 

X  ( 1  —  £/*<*""*  U'*"'m'‘dXTlRe**U')lXl "l'  •• 

•“  I +  //>„. 

The  value  of  the  form  parameter  is  determined  by  formula 

_ fr** 


(5.75) 


where 


*•  nVO  U*) 


I 

*-  »* 


±  r 

d.X 


(5.76) 


The  critical  value  of  the  form  parameter  is  determined  from  the 
graphs  of  Fig.  39. 

The  drag  coefficient  is  calculated  by  formula 

=  (5.77) 

## 

where  cfQ  is  referred  to  the  value  of  Jte  ,  determined  by  equation 
(5.74)  or  (5.75). 

In  practice^  equation  (5.75)  can  be  used  and  in  the  case  of  not 

very  strong  changes  of  the  wall  temperature  along  the  length  of  the 

contour.  Connected  with  this  is  the  fact  that  errors  in  the  deter- 
* 

** 

mination  of  the  value  of  Rt  weakly  affect  the  value  of  c^q  due  to 
the  small  value  of  the  exponent  m  in  the  law  of  resistance  for  a 
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iu.  *+'  mmt. 


turbulent  boundary  layer. 


5.11.  Solution  of  the  Equations  of  Energy  and  Momentum 
for  Ax  1  symmetric  Flow  of  a  3a s 

For  an  axisymmetric  boundary  layer  the  equation  of  energy  and 
momentum  can  be  written  in  the  form: 

4**"  -  i  dir  i  art 


it.  \ 


*”147 


I 

MX  ft 

,  He1"  Jft 
ft  ‘  MX 


itX 


)  Re,  Sl  . 


Re,  F\f\ 


(5.78) 

(5.79) 


4X  . 

where  R  is  the  current  radius. 

The  law  of  heat  transfer  remains  practically  the  same,  as  for  a 
flat  boundary  layer. 

The  integral  of  equation  (5.78)  differs  from  the  integral  of  the 
energy  equation  of  a  flat  boundary  layer  only  by  the  fact  that  the 
product  RAT  enters  into  it  instead  of  the  value  of  AT. 

The  integral  of  equation  (5.79)  at  Twall  -  const  has  the  same 
form  as  formula  (5.75),  but  instead  of  the  value  of  UK  in  it  one  should 
put  the  product  RUx.  During  subsonic  flow  in  formula  (5.67)  the 
quantity  «0  18  replaced  by  Rw^. 

During  flow  of  a  gas  in  a  cylindrical  nozzle  (internal  problem) 
the  integral  of  the  energy  equation  is  conveniently  written  in  the 
form 


ReT  =  [ 


«  _  r  0  +  «)**»,>,  (  d»,  y 4  m(±z. L\  (  -  \  h 

2  A  T*  +  "  '  D  J  1  *  -  I  /  [  ft  +  I  ' 

.<  J  r,  a  t*  *  *  j'  ‘ 


(5.80) 


Besides,  it  is  assumed  xhat  the  turbulent  layer  starts  from  the  inlet 
section  of  the  nozzle. 

Equation  (5.80)  is  obtained  from  the  basic  solution  of  the  energy 
equation,  if  we  bear  in  mind  that  during  flow  of  a  gas  in  a  nozzle, 
with  a  good  degree  of  accuracy,  it  is  possible  to  assume 
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(5.81) 


I 

U{\-U*)rT 


1 

*+ 1 


2 

*  +  I 


•eby  in  this  case  (ftCI/ft)  ■  (DCI/I>)2, 

*e  ft  is  the  cross-sectional  area.  In  these  formulas  is  the 

cr 

leter  of  the  critical  cross  section  of  the  nozzle. 

The  value  of  the  local  Nusselt  number  will  be  determined  by 
mla 

*•-=  f  ***■*- «v*r  (2=ff  (5.82) 

:e 

"-“-XT' 


characteristic  linear  dimension  L  should  be  one  and  the  same  in 

VU  and  Re  numbers. 

* 

The  quantity  f  during  flow  of  a  gas  in  a  nozzle  is  a  single- 
led  function  of  the  dimensionless  cross-sectional  area  of  the  noz- 
that  essentially  simplifies  calculation  of  the  function  ^(X)  and 
sulation  of  the  thermal  boundary  layer. 


fig.  41.  Comparison  of  calculations  of  heat  transfer  in  a  super¬ 
ionic  nozzle  by  formula  (5. 82)  with  the  experimental  data  of 
l.  I.  Leont'ev. 


In  Fig.  41  is  given  a  comparison  of  the  calculations  of  heat 
isfer  in  a  supersonic  nozzle  by  the  proposed  method  with  the  results 
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of  an  experimental  investigation.  It  is  possible  to  ascertain  a  very 
good  agreement  of  theory  and  experiment. 


5.12.  Solution  of  the  Equations  of  Energy  and  Momentum  During 
Subsonic  frlow  of ~a  das  on  a  Permeable  slightly 

TEurved  Surface 

The  laws  obtained  in  Chapter  IV  of  friction  and  heat  exchange 
for  a  permeable  plate  can  be  used  for  solution  of  the  equations  of 
nomentum  and  energy  during  streamlining  by  a  subsonic,  nonisothermal, 
turbulent  boundary  layer  of  a  permeable  slightly  curved  surface. 


The  energy  equation  for  the  considered  case  is  written  in  the 


rorm 


4X 


rtS  T 


(5.83) 


(here  the  function  ?  is  determined  by  formula  (4.30)»  i.e., 

)’  (5.84) 

rhere  is  chosen  from  formula  (3.47). 

In  the  general  case,  the  quantities  b^,  and  b^  are  func- 
ions  of  the  coordinate  X.  Then  the  integral  of  equation  (5.85)  has 


he  form 


=  Tr  *  h  ('  -  ~£r)' +  H * 

X, 

'  1 

XAt  ’  mdX  ~  |*»r4  7)1.*  *)  *  "• 


(5.85) 


or  the  given  functions  wQ(x),  AT(x),  and  bt(x)  the  value  of  X«t  is 
alculated  from  equation  (5.85)  as  a  function  of  the  coordinate  X. 
Distribution  of  the  feed  of  the  cooling  gas  along  the  length  of 


he  circuit  is  calculated  by  formula 


_2l_  =  _5ei_s*#*T. 

t9 , 

The  local  Stanton  number  is  determined  by  formula 


SI  =1,(1- 


1 *«.. 


(5.86) 


(5.87) 
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The  local  heat  flux 


#« = i**H  *#  A  TSi.  (5.88) 

From  the  known  consumption  of  the  cooling  gas  Its  initial  temper¬ 
ature  is  determined t  , 


r, 


*w«  ’ 


(5.89) 


During  evaporation  and  sublimation  the  surface  temperature 
usually  is  given.  The  quantity  of  gas  (vapor)  discharged  from  this 
surface  is  subject  to  determination. 

In  the  considered  case 

=  /'!/,.  (5.90) 

where  r  is  the  latent  heat  of  vaporization. 

Hence 


r 


Si 

SI.. 


-  (5.91) 


where  K  »  (r/cplAT)  is  the  criterion  of  phase  transition. 

On  the  other  hand,  the  value  of  F  is  determined  by  expression 
(5.84).  Combining  formulas  (5.84)  and  (5.91)*  we  find  that 

-')•  (5.92) 


Introducing  the  value  of  T  from  (5.91)  and  bt  from  (5 
(5.83),  we  find 


92)  in  equation 


(5.95) 


From  the  found  value  of  the  value  of  fltQ  is  determined; 
the  value  of  b^  from  formula  (5*92);  and  from  formula  (5*86)  the 
flow  rate  of  the  cooling  gas. 
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For  relatively  small  intensities  of  evaporation  or  sublimation 

(which  is  usually  the  case)  it  is  possible  to  be  limited  during  cal- 

## 

culation  of  to  the  first  approximation,  i.e.,  to  conduct  a  cal¬ 
culation  of  this  quantity  by  formula  (5.40)  taking  into  account  the 
nonisothermalness  of  the  flow  from  Section  5.6. 

We  will  consider  the  case,  when  the  parameters  are  given  of  the 
external  flow,  the  wall  temperature,  and  the  initial  temperature  of 
the  cooling  gas,  passing  into  the  boundary  layer  through  a  porous 
surface.  It  is  necessary  to  determine  the  needed  distribution  of 
the  cooling  gas  along  the  circuit  of  the  warmed  surface.  In  this 
case 


r=b, 


AT, 

IT 


(5.94) 


where  AT^  *  Twall  “  T1  and  T1  is  the  temperature  of  the  gas  at  the 
entrance  to  a  porous  surface. 

The  value  of  b^  will  be  determined  by  formula  (5.92)  upon  replace¬ 
ment  in  it  of  the  value  of  the  criterion  K  by  the  ratio  (AT^/AT).  The 

## 

value  of  Re^  will  be  determined  by  formula  (5*95)  upon  the  same 
replacement  of  K  by  (AT^/AT).  Further  calculations  are  also  con¬ 
ducted  just  as  in  the  preceding  problem. 

For  an  axisymmetric  boundary  layer,  solution  of  the  energy  equa¬ 
tion  has  the  same  form  as  in  the  case  of  a  flat  boundary  layer,  but 
instead  of  AT  the  product  RAT  enters,  where  R  is  the  radius  of  the 
streamlined  body.  The  whole  remaining  course  of  calculation  does  not 
change. 

The  momentum  equation  for  a  flat  boundary  layer  on  a  permeable 
wall  of  small  curvature  has  the  form 


war* 

rfjr 


4X 


<1 


(5.95) 


=  <**  + *. 
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wh  re 


■-M" --tr- 


(5.96) 


We  have  the  dependences  *t(x),  b(x),  bcr(x).  Assuming  H  -  HQ 
1.5,  we  obtain 


^[iY(i --£)%•*] x 

I 

xrfx.HH»-Se*)Vr“} 


(5.97) 


For  an  axisymmetrlc  boundary  layer  the  solution  remains  the 
same,  but  instead  of  the  quantity  the  product  appears. 

The  local  value  of  the  coefficient  of  friction  will  be  determined 
by  formula 

-£-)V  (5.98) 

When  the  distribution  of  the  flow  rate  of  injected  gas  j^(x) 
along  the  contour  of  a  streamlined  body  is  known,  the  problem  is 

solved  by  the  method  of  successive  approximations.  The  first,  and 

** 

fully  sufficient  approximation  for  the  determination  of  Re^  is 
calculated  by  the  formulas  for  an  impenetrable  wall  of  small  curvature. 

The  methods  presented  allow  calculation  of  heat  transfer  and 
friction  on  a  porous  surface  of  the  frontal  part  of  bodies  (flame 
regulators,  the  frontal  parts  of  a  sphere  and  a  transversely  stream¬ 
lined  cylinder,  the  head  part  of  axisymmetrlc  bodies). 
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CHAPTER  VI 
FLOW  IN  A  PIPE 

Definitions  of  Cyrillic  Items  in  Order  of  Appearance 
ct  ■  wall  *  wall 
t  -  t  «  Thermal,  Turbulent 

6.1.  Distribution  of  Velocities.  Friction,  and  Heat 
Transfer  During  Quasi isothermal.  Stabilized  Flow 

During  flow  in  a  pipe  a  flow  stabilization  section  and  a  section 
with  stabilized  flow  are  distinguished.  Stabilized  flow  approaches 
after  merging  of  the  boundary  layers,  which  arise  in  the  inlet  section 
of  the  pipe,  and  for  isothermal  conditions  is  characterized  by  con¬ 
stancy  of  all  the  flow  parameters.  In  the  inlet  section  (stabiliza¬ 
tion  section)  the  flow  parameters  are  changed  in  connection  with  the 
build-up  of  the  boundary  layer  on  the  wall  of  the  pipe.  A  diagram 
of  the  flow  in  the  inlet  section  of  the  pipe  is  depicted  in  Fig.  42. 

During  stabilized  isothermal,  tur¬ 
bulent  flow  the  velocity  distribution  In 
a  smooth  pipe  is  well  described  by  formula 

Fig.  42.  Diagram  of  the 

build-up  of  the  boundary  (5.1).  This  is  explained  by  the  fact  that 
layer  in  the  Inlet  sec¬ 
tion  of  a  cylindrical  flow  in  a  pipe  is  nozzle  flow  with  com- 

plpe. 

paratively  small  values  of  the  form  para¬ 
meter.  The  pressure  drop  in  a  pipe  is  determined  by  formula 
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where  C  Is  the  drag  coefficient; 

D  is  the  internal  diameter  of  the  pipe; 
w  is  the  average  consumption  flow  velocity; 

pQ  is  the  density  of  the  flowing  medium,  referred  to  the  average 
u  consumption  temperature  of  the  stream. 

During  isothermal  flow  p  *  pQ  *  const. 

Prom  the  condition  of  equilibrium  (Fig.  43)  it  follows  that 
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(6.2) 


(6.5) 


Fig.  43.  Diagram 
of  the  action  of 
forces  on  an  ele¬ 
mentary  cylinder 
of  fluid  in  a  pipe. 


The  average  consumption  velocity 

<6-4) 

ii 


Putting  in  (6.4)  the  velocity  profile  (3.1)  and  assuming  p  - 


*  const,  we  find  that 


where 


From  here  the  law  of  resistance  follows: 


(6.5) 


— =  0,8*  In  (Rt0  \  )  —  0.9. 

I  C 


(6.6) 


In  the  region  5*10^  <  ftp  <  1*10^  the  Blasius  formula  gives  good  results 


c=  mi 


(6.7) 
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where 


A*o  = 

The  velocity  distribution  to  the  n  ■  1/7  power  corresponds  to 
formula  (6.7). 

The  form  parameter  f  can  be  calculated  by  formula  (2.36), 


/=- 


4x 


v* 

l 


IL 

2 


(6.8) 


In  the  region  of  stabilized  flow,  when  the  boundary  layers  merged, 

2t 

6  «  Rq.  According  to  formula  (6,2)  —  ^ 


Hence  during  Isothermal  stabilized  flow  in  a  pipe 


At  n  -  1/7 


4 

i.e.,  at  HeD  >  10 
At  iteD  >  105 


f — — 


/<0,003; 


/<0,002. 


(6.9) 


(6.10) 


(6.11) 


(6.12) 


The  law  of  heat  transfer  in  the  region  of  stabilized  flow  of  a  gas  in 
a  pipe  is  described  by  formula. 


where 


Correspondingly 


N*o  =  0,023  Pr**  R*)?. 


Nmn  — 


51  = 


o.sa 

to”*# 


(6.13) 


(6.14) 


-115- 


6.2.  Nonisothermal,  Stabilized  Flow  of  a  Gas 

Stabilized  flow  In  a  pipe  of  constant  cross  section  can  exist 
only  at  M  <  1. 

Inasmuch  as  in  this  case  the  boundary  layers  are  merged,  then 
condition  la,  *  idem  is  equivalent  to  the  condition  Re**  »  idem.  There¬ 
fore,  the  limiting  formula  (2.47)  for  given  case  takes  the  form 

(6.15) 

In  Fig.  44  is  given  a  comparison  of  calculations  by  formula 
(6.15)  with  experimental  data.  The  limiting  formula  in  this  case 
also  correctly  depicts  the  real  character  of  the  investigated  relation. 
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Fig.  44.  Influence  of  temperature  factor 
on  friction  and  heat  transfer  during 
stabilized  flow  of  gas  in  a  pipe:  1  —  cal¬ 
culation  by  formula  (6.15)#  ®  —  Jordan, 

O  —  Oreber,  •  —  II*  in,  A  —  Ivashchenko, 
pr  —  NASA  (averaged  by  series).  Experiments 
on  the  resistance  of  friction:  O  —  II* in. 


If  we  consider  the  finiteness  of  It  the  numbers,  then  we  should 
use  formula  (2.58),  which  at  f*  M  1  takes  the  form 


_ 9 _ 

+i 


(6.16) 


For  stabilized  flow  in  a  pipe  this  formula  should  be  written 


with  consideration  of  the  ratio 


„Ji,  inasmuch  as  in  the  theory  of  a 
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In  the  inlet  section  of  the  pipe  a  boundary  layer  is  developed 
just  as  during  external  streamlining  of  a  body  as  long  as  the  opposite 
points  of  the  external  boundary  of  the  layer  do  not  touch. 

The  over-all  flow  rate  of  a  gas  along  the  whole  pipe  is  constant, 
and  the  flow  velocity  in  the  undisturbed  nucleus  of  flow  changes  due 
to  the  build-up  of  the  boundary  layer.  Taking  into  account  this 
circumstance  all  the  equations  of  a  turbulent  boundary  layer  derived 
above  can  be  used. 
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We  will  assume  that  upon  entry  Into  the  pipe,  the  velocity  dis¬ 
tribution  Is  uniform,  and  a  turbulent  boundary  layer  Immediately 
arises.  The  latter  condition  Is  always  fulfilled  because  of  the 
turbulence  causing  action  of  the  finite  thickness  of  the  leading  edge. 

We  will  be  limited  by  the  case,  when  the  pipe  is  warmed,  starting 
from  the  leading  edge.  The  build-up  of  dynamic  and  thermal  boundary 
layers  occurs  simultaneously  along  the  whole  length  of  the  inlet 
section  of  the  pipe.  Correspondingly,  the  parameters  of  the  gas  In 
the  undisturbed  nucleus  of  the  flow  also  change. 


The  equation  of  continuity  of  the  flow  will  be  written  in  the 


form 


fa  wt,  Ft%  -  2 1  ?wdR  =  const, 


(6.18) 


where  pn<1  and  wn<t  are  the  density  of  the  gas  and  flow  velocity  in  the 
1  1  Inlet  section  of  the  pipe. 

For  a  cylindrical  pipe  according  to  formula  (1.59) 


‘•■K1- SI!' --£-)*■ 


(6.19) 


Taking  into  account  this  expression  of  the  thickness  of  displace¬ 
ment,  it  is  possible  to  bring  the  equation  of  continuity  (6,18)  to 


the  form 


—  Por«  (l  *"  2  ^ 


(6.20) 


where  p0  and  wQ  are  the  density  and  flow  velocity  outside  the  boundary 
u  u  layer  in  the  cross  section  x. 

The  density  of  the  gas  pQ  is  raalted  to  the  density  p00  for  the 

parameters  of  stagnation  (p00;  To)  and  dimensionless  velocity  U 

by  known  the  expression 


=  (1  -  U'Y 


(6.21) 
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ij*  a  iwmumtmMmmm* 


* 


"TUI***- 


which  has  already  been  used  In  Chapter  V  during  derivation  of  formula 


(5.73). 

Combining  equations  (6.20)  and  (6.21),  we  find  that 


where 


(6.22) 


(6.23) 


Introducing  the  parameter  H,  we  obtain  the  value 


i  l 

*.*•=-^[4/(1  -  t/,(l  ~0?)*-,l> 


(6.24) 


where,  as  before 


/?#••  =  . z=  -iaSaL*, 

Mi  l*w 


The  momentum  equations  will  take  the  form 


i 

—  r  —*  •  -^r  (I  +H)  =  Ui\-  {/»)*-'  X*T-^- 


4X 

where 


4X 


fc.il*- ')  ‘  (6.25) 


x  = 


We  will  be  limited  by  flows  for  not  very  large  values  of  the  M 
number.  Then  supersonic  flow  of  a  gas  in  a  cylindrical  pipe,  which 
is  always  diffusion  flow,  flows  at  values  of  the  form  parameter  f, 
not  close  to  the  critical  value.  In  this  case  the  function  I  can  be 
equated  to  the  function  i.e.,  will  depend  only  on  the  parameters 
if*  and  i/.  On  these  parameters  also  depends  the  quantity  H. 

As  L.  E.  Kalikhman  showed,  the  quantity 


(6.26) 


depends  on  the  parameters  of  nonisothermalness  significantly  weaker 
than  the  quantity  H. 

Between  the  quantities  H  and  H'  there  exists  the  dependence 
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•»  * 


Ml  i 


| 


where 


H=  - ! - (w- :  /»'  -  H\  , 

i-w  \  r%  ) 


(6.27) 


(6.28) 


Upon  solution  of  equation  (6,23)  it  is  possible  in  practice  to  assume 


*e*# 


Re 


jL_  *  Hq  during  isothermal  flow,  i.e.,  1.3, 


Then 


He**  = 


» 

(1 

-£/»)*  “‘l/  (1- £/})*"' (/Ji 

(!  V*) 

4 

• 

(/»  +  1.3  | 

» 

C  -  r„  \ 

] 

K  ) 

The  momentum  equation  will  be  rewritten  thus; 


(6.29) 


(6. jo) 


HZ -  +  -gL- [,  v ,,3(,  -  .<r« - \]  JS. B 

<v  t/(i-r>)  l  \  t;  j] 


1 

f  c/d-c/y7'*,,. 


(6.31) 


*» 


Putting  in  this  equation  Re  from  equation  (6,30),  we  obtain 


[ 


i  _ 


mi 

4--^) 


’M-V1) 


l/»]  •  |(1  -  U8)77'  -  ~j-  U*\\ -  + 

+  24/0,  (l -*/?)”']  "  |C1  — 
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f  ■  I  1  J  I  • 


(6.32) 


For  a  given  law  of  change  of  the  wall  temperature  and  known 
velocities  and  stagnation  parameters  at  the  inlet  in  the  channel, 
equation  (6.32)  allows  us  to  determine  the  law  of  the  velocity  change 
in  the  undisturbed  nucleus  of  the  flow  along  the  length  of  the  pipe 
X. 

## 

From  equation  (6.30)  the  change  of  Bi  along  the  length  of  the 
pipe  is  calculated  and  by  formula  (5*77)  the  local  values  of  the 
coefficients  of  friction  are  determined. 

The  values  of  the  Nusselt  number  are  determined  by  formula 


= -*■  (6-3?) 

It  is  necessary  to  note  that  the  presented  method  of  calculation 
can  be  applied  only  if  Impactless  entrance  of  a  gas  into  a  cylindrical 
pipe  takes  place.  For  the  region  of  subsonic  gas  velocities  it  is 

p 

possible  to  consider  U  «  1.  In  this  case  we  have  T  ^  ■  const 


where 


(6.34) 


«V 
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Assuming  m  ■  0.25#  B  -  0.0258#  and  taking  the  limiting  value  of 
the  function  for  subsonic  flow,  we  obtain 


X- 

X 


,T(^  ir;. 

few  -  2  (ew  *  •  )**  1 1 


2  #rctf  X 


I  2  fc-  if* 

i-fe-.r 


( ,  T  *  ’ 


(6.55) 


where 

ft. 


Fig.  46.  Dependence  of  the 


•0  PR 

parameter  on  wo 

1 

and  it  by  equation  (6.55)*. 
1)  ^  -  0.4;  2)  it  *  0.6; 

3)  *  -  0.8;  4W  -  1.0; 

5)  ^  *  1.4;  6W  «  1.8; 

7}*  •  2.2;  8}  f  =  2.6; 

9)  V'  «  3.0;  10)  it  -  4.0. 


eu«  D 

• 

In  Fig.  46  Is  shown  the  depen- 

-0  28 

dence  of  the  parameter  X-*V  on 
Wq  and  it,  calculated  by  equation 
(6.55). 

Using  this  graph,  we  can  deter¬ 
mine  the  velocity  change  In  the  nu¬ 
cleus  of  flow  along  the  length  of  the 
pipe  for  given  parameters  of  the  gas 
at  the  pipe  Inlet  and  for  a  given 
wall  temperature. 

For  subsonic  velocities  from 
equation  (6.50)  we  have 

(6.36) 

*.*♦ 


The  local  coefficients  of  friction  are  determined  form  the  law  of 
resistance: 


o.osn 

(/T  +  i),jr***M5  ’ 


(6.37) 
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Tne  local  coefficients  of  heat  transfer  are  calculated  by  formula 
(6.33).  Using  equation  (6.35)  from  equation  (5.40),  generalized  for 
nonlsothermal  flow,  we  obtain  the  value 


* 

D«"  _  ,  |,2+ 1*4) _ 

*.*♦  I'  ^T’  (S*-!)1* 

[■ 


a | * («••.- 1)**  -  r»  wo,  > 

i  -  («,  -  ir" 

st  (fc-ir-.T{i-ir»+r 


1M 


(6.38) 


Thus,  during  flow  of  a  gas  (P*  *  1)  In  the  Initial  section  of 
the  pipe  we  have 

I.2S+I.C24 
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ST  (5,-ir-»^T(5.-i>**+i 
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(6.59) 


R» 


** 


Fig.  47.  Dependence  of 


6** 


on 


X*HtI0,2^  and  by  formula  (6.39): 
U1 

1)  V'  *  0.4;  2)  V'  =  0.6j  3)  V'  ■  0,8; 

*)  *  -  |)  *  -  6)  *  "  1-8J 

7)  V'  -  2.2;  8)  V'  =  2.6;  9)  V'  “  3.0; 

10)  ^  =  4.0. 


In  Fig.  47,  Is  shown  the 

on  X»ReJ^#2^  and  f.  As 
Re  Dl 

can  be  seen  from  the  graph, 
during  subsonic  gas  velocities 
in  the  initial  section  of  a 
cylindrical  pipe,  simultaneous 
build-up  of  thermal  and  dynamic 

layers  1>  <  1  and  R»*  »  li#*. 
Taking  into  account  this  cir¬ 
cumstance  it  is  possible  to 
propose  a  convenient  method  of 
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generalization  of  the  experimental  data  on  haat  transfer  i..  the  initial 
section  of  a  cylindrical  pipe.1 

From  equation  (5*39)  It  follows  that 


t'p  *a  r 


(6.40) 


Further  we  have  the  relation* 8 


'  '•  (6.41) 

I  3 

**o  ~  (6.42) 


=  I  IP\ 


i  - 1 

m 


(6.43) 


Thus,  by  measuring  the  distribution  of  the  static  pressures,  the  wall 
temperatures,  and  the  heat  flux  along  the  length  of  the  pipe  by  these 
equations,  it  is  possible  to  construct  the  experimental  dependence 
of  St  on  ®s£*  . 

If  in  the  experiments  measurements  of  the  static  pressures  are 

_  •• 

not  carried  out,  then,  by  taking  into  account  the  condition  m 
«  !•**,  it  is  possible  to  calculate  the  Stanton  number  for  subsonic 
velocities  by  the  formula 


S#* 


_ _ 

(ft#,  +  5.2  i  ftrf“)  f  «>,*..  4  T 


(6.44) 


In  Fig.  48  are  shown  the  results  of  treatment,  by  the  proposed 
method,  of  the  experiments  of  B.  S.  Petukhov,  V.  L.  Lel*Chuk,  B.  V. 
Dedyakin,  V.  N.  Fedorov,  A.  I.  Leoont'ev,  I.  A.  Kozhinov,  and  S.  I. 
Kosterin.  All  these  data  are  reduced  to  the  conditions  f  *»  1  by  the 
limiting  formula  (3.55).  As  can  be  seen,  the  average  line,  drawn 
through  all  the  experimental  points,  is  described  by  formula 


1Method  was  developed  Jointly  with  V.  K.  Fedorov. 

2During  derivation  of  formula  (6.42)  it  is  assumed  that -ii-  * 
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(6.»5) 


*  = 


«■<!» 

**n*»7“' 


which  coincides  with  the  relation 
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ipf7'  * 


(6.46) 


where  cfQ  Is  calculated  by  formula  (3.8)  upon  substitution  In  It  of 
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Fig.  48.  Results  of  treatment  of  the  experimental  data  on 
heat  transfer  In  the  initial  section  of  a  cylindrical  pipe: 

1  -  turbulent  conditions  St  -  1  g*  2  "  laminar 

Vl  *  “'tot  * 

conditions  St  *  ;  •  —  experiments  of  B.  V. 

W  tor 

Dedyakin  and  V.  L.  Lel'chuk;  A  —  experiments  of  B.  S. 
Petukhov;  O  —  experiments  of  I.  A.  Kozhinov,  S.  I.  Koaterin, 
Leont'ev,  and  V.  N.  Fedorov. 


It  is  necessary  to  emphasize  that  these  experimental  data,  pre¬ 
sented  by  their  authors  in  the  usual  a.  -ensionless  number  treatment, 
diverge  noticeably. 

The  obtained  result  distinctly  confirms  the  generality  of  the 
laws  of  heat  transfer  and  friction  in  a  turbulent  boundary  layer  for 
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the  internal  and  exterior  tasks  of  aerodynamics. 


In  Fig.  49  are  shown  the  results  of  treatment,  by  the  proposed 
method,  of  the  experimental  data  on  heat  transfer  for  pipes,  nozzles, 
plates,  and  nose  cones.  These  experiments  embrace  a  wide  range  of 
changes  of  the  numbers  M  and  ty.  In  spite  of  a  slgnlficiant  variance 
of  the  points,  all  of  them  are  grouped  around  a  line,  corresponding 
to  formula  (6.46). 
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Fig.  49.  Results  of  generalization  of  experimental  data  on  con¬ 
vective  heat  transfer  in  a  turbulent  boundary  layer  of  compressible 
gas;  1  —  turbulent  conditions;  2  —  laminar  conditions;  H 
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The  equation  of  a  thermal  boundary  layer  can  be  written  in: 


4H ” 
4X 


Pt"  .  ~  (In  1 1  —»|  -  St  Pto, 


(6.47) 


where 


<• 


F#0  - 


«u 
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or  taking  Into  account  equations  (6.20),  (6.22) 


4X 


-  ^  rt  i*9a,  +M .  *t  **•), 


where 


For  the  case  i>  -  const 


WJT  Sr***  VJ  ' 


(6.48) 


For  subsonic  velocities  of  a  gas  (and  m  ■  0.25)  after  Integration  we 
obtain 


X  = 


1 

/  2  \ 

M 

S.*t 

/  la.i*)*** 


.....  ,  !  »*•*'} 

"  1  v  ti*  J 

1  +\ 
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i  ■  i 

r  **  vM 

”  -•  •  '[  ut  J 
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-r  2*rc  tf 


1  2  *1 

!«■  ' 

_ Mi _ _ 

/  * 

% _ 1 

|M  -  hr’**’1 

l  5.2  if  J 

(6.49) 


** 


From  equation  (6.49)  we  obtain  a  change  of  the  parameter  9b 
along  the  length  of  a  pipe  for  given  parameters  of  a  gas  at  the  inle'; 
of  a  pipe  and  for  a  given  wall  temperature. 

Local  values  of  the  Nusselt  number  and  heat  flux  are  determined 
by  equations. 


a(-^—Y 

"•  - -- ■  i  /»«». + 5.2*  •  /»»••! 


(6.50) 
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and 


*'l7r7r H*- 

*'m  ~  W*  r>  ■  pPA  o 

X  [Pin,  -r  5.2^*  •  Pt**). 


(6.51) 


For  the  case  of  a  given  constant  thermal  load  along  the  length 
of  a  pipe  from  equation  (6.47)  we  obtain 


Pt**.  N*x 

‘  <1  4)  * 


(6.52) 


where 


N*\  ■- 


1n-f> 

k  •  r.« 


In  the  range  of  a  change  of  V'  from  0  to  3.0  the  function  can 
be  expanded  In  a  series  and  can  be  bounded  by  the  first  term,  i.e., 
assume  (3.49): 


Then 


St  = 


2P 


Pl**W  (|  ’\)Pr"r' 


__ _ N* 1_ 

*/><«  v) 


(6.55) 


From  equations  (6.20)  and  (6.22)  It  follows  that 


P*d-P*ik  +  5,2*  •  Pe' 


(6.54) 


As  a  result  we  have  a  system  of  three  equations  (6.52),  (6.55)# 

*# 

and  (6.54)  with  three  unknowns  if,  Ft  and  Ftp. 

*# 

For  determination  of  the  change  of  Ft  along  the  length  of  the 
pipe,  from  equations  (6.52),  (6.55) ,  and  (6.54)  we  have 
X  =  3^6  Afaf'fAto,  PS*npf*  -  0,0286 /Vd,- 0,1 49  F#**- 

(6.55) 
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Knowing  local  the  values  of  fie  ,  from  equation  (6.52)  we  de¬ 
termine  the  local  values  of  the  wall  temperature: 


♦  = 


I  + 


X. 


The  proposed  method  of  calculation  can  be  also  be  extended  to 
the  case  of  an  arbitrary  law  of  supply  or  removal  of  heat  along  the 
length  of  the  channel.  Only  in  this  case  instead  of  equation  (6.52) 
we  should  use  equation  (6.47).  Further  considerations  remain  the 
same. 


6.4.  Heat  Transfer  to  a  Vapor  During  Critical  Pressures 

The  theory  of  the  limiting  laws  of  friction  and  heat  transfer 
can  be  applied  also  to  flows  of  gases,  not  obeying  the  Clapeyron- 
Mendeleev  equation  of  state.  The  problem  is  solved  most  simply  for 
the  bounded  temperature  Integrals  in  that  case,  when  there  is  possible 
a  linear  approximation  of  the  dependence  of  the  gas  density  on  the 
temperature  (or  enthalpy). 

In  this  case  for  a  uniform  medium 


(6.56) 


where,  as  in  Chapter  IV, 


♦.  - 


The  limiting  law  of  heat  transfer  at  c_  »  const  for  subsonic 


flow  is  expressed  by  formula 


9t 

J*. 


=H=L— v 


(6.57) 


For  an  ideal  gas 


»» 
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In  more  complicated  cases,  for  instance  during  a  significant 
change  of  the  heat  capacity,  it  is  possible  to  obtain  if  not  exact 
quantitative,  then  in  any  case  correct  qualitative  results. 

Let  us  consider  in  this  plan  the  character  of  the  changes  of  the 
coefficient  of  heat  transfer  during  flow  of  a  vapor  in  the  region  of 
near  critical  (in  a  thermodynamic  sense)  parameters. 


Fig.  50.  Change  of  coefficient  of  heat  trans¬ 
fer  during  flow  of  water  and  steam  in  the 
region  of  critical  parameters  from  the  experi¬ 
ments  of  A.  A.  Armand  et  al.  :  1  —  P  =  240  atm; 

2  —  P  =  260  atm;  q  =  1.7  •  1C)5  kcal/m2  •  hr; 

7w  =  650  kg/m2  .  sec. 

In  Fig.  50  are  shown  the  results  of  the  experiments  of  A.  A. 

Armand  on  the  determination  of  the  coefficients  of  heat  transfer 
during  flow  of  steam  having  supercritical  parameters  in  a  cylindrical 
pipe.  As  can  be  seen,  during  passage  through  the  critical  temperature 
there  is  observed  a  sharply  expressed  maximum  of  the  coefficient  of 
heat  transfer.  Such  a  course  of  the  curve  a(T)  is  analogous  to  the 
course  of  the  dependence  c  (T)  in  the  region  of  critical  parameters. 

r 

In  Fig.  51  is  shown  a  comparison  of  the  results  of  the  experiments 
of  Z,  L.  Miropolskiy,  M.  A.  Styrikovich,  and  M.  E.  Shitsman  with 
calculations  by  formula  (6.15).  No  unambiguous  relation  between  the 
formula  and  the  experimental  data  on  the  graph  is  observed.  An  ana¬ 
logous  result  was  also  obtained  in  the  above-mentioned  work  of 
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Fig.  51.  Comparison  of  experiments  of  Z.  L. 
Miropolskiy,  M.  A.  Styrikovich,  and  M.  E. 
Shitsman  with  Formula  (6.15). 


A.  A.  Armand.  At  the  same  time,  according  to  the  experiments  of 
V.  E.  Doroshchuk  and  V.  L.  Lel’chuk,  at  rf/  *  1  formula  (6.15)  is  also 


valid  in  the  region  of  supercritical  parameters  (Fig.  52).  Conse¬ 


quently,  deviation  of  the  experiments  with 


Fig.  52.  Heat  transfer 
for  critical  parameters 
and  quasi -isothermal 
conditions  from  the  ex¬ 
periments  of  V.  E. 
Doroshchuk  and  V.  L. 
Lel'chuk. 


calculations  by  formula  (6.19)  are  connected 
with  the  nonisothermalness  of  the  stream. 

During  simultaneous  change  of  the 
quantities  cp  and  p  the  intensity  of  tur¬ 
bulent  heat  transfer  will  be  determined  by 
the  relation 


9  =!>/’ 


iw 

*1 


(6.58) 


where  i  is  the  specific  enthalpy. 

The  limiting  relative  law  of  heat  trans¬ 


fer  is  derived  from  equation  (6.58)  by  the  already  known  method: 


St 

st. 


(6.59) 
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But  in  distinction  from  equation  (2.69)  in  the  given  equation 


St  = 


»= 


Air*  . 

itr-U  ’ 


_ Is. _ 


(6.60) 


In  the  region  of  near  critical  parameters  the  specific  volume 
of  the  vapor  is  almost  a  linear  function  of  the  enthalpy,  i.e.,  for 
definite  intervals  of  Ai  formula  (6.47)  can  be  applied,  if  in  it  we 
introduce  the  value  of  $  from  (6.60).  Then  the  integral  (6.59)  is 
expressed  by  formula  (6.57)*  if  we  compare  the  heat  flux  Qwall*  We 
have 


-5T  -  -(  /ir+1  )••  (6.6i) 

Multiplying  and  dividing  the  left  part  of  formula  (6.6l)  by  At, 
we  find  that  it  is  applicable  to  the  ratio  of  the  coefficients  of 
heat  transfer. 

In  Figs.  55  and  54  are  given  the  results  of  the  treatment  of  the 
experimental  data  on  heat  transfer  to  vapors  of  water,  carbon  dioxide, 
and  oxygen  in  the  region  of  near  critical  parameters.  The  logarithmic 
straight  line,  drawn  on  this  graph.  Is  calculated  by  the  formula 

A =  0,023  Piv4  Re;9 1 - ^ — 

\  l  *.  + 1 

Physical  properties,  entering  into  the  criterions  !G1q,  PPq  and 
Re  are  referred  to  the  average  stream  temperature  by  the  best  content. 

It  is  possible  to  ascertain  that  the  introduction  of  a  correction 
by  formula  (6.6l),  to  a  significant  degree,  approaches  the  experimental 


(6.62) 
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data  with  the  formula  for  quasi-isothermal  conditions.  Thus,  also 
in  this  case  the  influence  of  the  nonisothermalness  basically  is  ex¬ 
pressed  via  a  change  of  the  vapor  density  along  the  cross  section  of 
the  turbulent  stream. 

V.  V.  Krasnoshchekov  and  A.  I.  Protopopov  proposed  the  empirical 
dependence 


st 

J4 


-  /  *  >°‘u/  >*t  \**'i  e„ 

l  /  l  K  I  Uf.  )  ' 


(6.65) 


where 


f  _  4»‘-4 

f‘~  T„~T. 


Fig.  53.  Comparison  of  experimental 
data  on  convective  heat  transfer  in 
the  region  of  critical  parameters  with 
the  limiting  formula  (6.62),  O  -  ex¬ 
periments  of  Z.  L.  Miropolskyi  and 
M.  E.  Shitsman  (HgO);  o  —  experiments 

of  A.  A.  Armand  (H^O);  •  —  experiments 

of  Bringer  (C02);  □  —  experiments  of 

Dickinson  (HgO);  x  —  experiments  of 

V.  S.  Protopopov  (C02). 
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In  the  region  of  the  near  critical  state  of  the  vapor  the  first 
of  these  factors  is  practically  equal  to  1. 

Formula  (6.63)  gives  quali¬ 
tatively  the  same  relation  between 
the  change  of  the  coefficient  of 
heat  transfer  and  the  temperature 
difference,  as  the  limiting 
formula  (6.62).  It  is  possible 
to  note  in  connection  with  this, 
that  during  criterial  treatment 
it  is  expedient  to  introduce  the 
ratio  -HSii 

Naturally,  the  given  analysis  cannot  pretend  to  a  full  reflection 
of  the  whole  variety  of  factors,  determining  such  a  complicated  pro¬ 
cess,  as  heat  transfer  for  near  critical  parameters  of  the  vapor. 
However,  from  data  presented  it  is  clear  that  at  least  one  of  the 
important  sides  of  this  phenomenon  is  explained  by  the  theory  of  the 
limiting  relative  laws  of  heat  transfer  in  a  gaseous  medium. 


Fig.  5^.  Comparison  of  experi¬ 
mental  data  on  convective  heat 
transfer  to  oxygen  in  the  region 
of  critical  parameters  with  the 
limiting  formula  (6.53).  —  by 

formula  (6.52);  O  —  experiments 
of  Powell. 


7**  • 
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CONCLUSION 

Definitions  of  Cyrillic  Items  in  Order  of  Appearance 
ct  =  wall  =  wall 
t  =  t  =  Thermal,  Turbulence 
Kp  =  cr  =  critical 

The  proposed  theory  of  the  limiting  relative  laws  of  friction 
and  heat  transfer  in  a  turbulent  boundary  layer  of  gas  is  based  on 
two  known  theoretical  conditions,  which  express  the  relation  between 
turbulent  friction  and  the  averaged  motion  in  a  flat  boundary  layer 
in  the  form  of  equations  (2.2)  and  (2.4). 

The  problem  of  the  influence  on  the  coefficient  of  friction  of 
the  nonisothermalness,  compressibility,  diffusion,  and  pressure  grad¬ 
ient  is  reduced  to  the  integral 


i 


where 


An  analogous  equation  is  also  obtained  for  the  relation 


This  integral  at  |»  -mb  in  a  number  of  cases  has  a  limiting  ex 
pression  of  the  form 
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i 


(2.30) 


Distribution  of  tangential  stresses  along  the  cross  section  of 

T 

the  boundary  layer  t  =  ■■■■■-■■■  ■  is  approximated  by  a  polynomial,  the 

Twall 

coefficients  of  which  are  calculated  from  the  boundary  conditions, 
ensuing  from  the  determination  of  the  boundary  layer. 

The  relation  between  the  gas  density  and  the  velocity  profile 
is  established  via  the  equation  of  state  and  coefficient  of  nonsimi¬ 
larity  of  the  stagnation  temperature  and  velocity  fields. 

Thus,  problem  is  reduced  to  an  analytic  investigation  of  the  in¬ 
fluence  of  the  nonisothermalness,  mass  transfer,  and  pressure  gradient 
on  the  relative  changes  of  the  coefficients  of  friction  and  heat 
transfer. 

Values  of  the  coefficients  cfQ  and  ttQ  can  be  determined  on  the 
basis  of  the  available  reliable  experimental  data  for  isothermal 
streamining  of  a  plate. 

For  streamlining  of  an  impenetrable  plate  by  an  unbounded  stream 
of  ideal  gas  an  exact  solution  of  the  limiting  equation  (2.30),  ex¬ 
pressed  by  formula  (3.47)  for  subsonic  and  by  formula  (3.48)  for 
supersonic  flows  is  obtained.  This  solution  does  not  contain  empirical 
coefficients  (including  "constants  of  turbulence")  and  is  not  connected 
with  any  sort  of  special  type  of  semi-empirical  theories  of  turbulence. 

if  ## 

The  weak  dependency  of  the  ratios  -  and  H  on  the  R»  and 

**  fO 

Rt.  numbers  allowed  us  to  extend  the  limiting  laws  to  flows  with 

finite  values  of  the  Reynolds  numbers. 

For  an  impenetrable  plate  in  practice  a  first  approximation 

expressed  by  formula  (3.58)  appears  to  be  sufficient. 
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In  the  case  of  a  permeable  surface  the  Influence  of  the  !•  number 
on  the  relative  laws  of  resistance  and  heat  transfer  is  still  less, 
due  to  perturbations,  introduced  in  the  near-wall  region  by  a  trans¬ 
verse  stream  of  substance.  Qualitatively  correct  and  quantitatively 
satisfactory  theoretical  formulas  are  obtained  in  an  approximation, 
expressed  by  the  equation  (for  a  plate) 


where  b  =*  ..fyiSj-1  is  the  wall  permeability  factor. 

cf0  p0  w0 

At  f  =  0  this  integral  gives  the  critical  value  of  the  factor 
b. 

It  is  shown,  that  for  subsonic  flows  of  a  gas  the  exact  limiting 
solutions  are  approximated  very  well  by  formula 


Y  <r 


=  !. 


(*.50) 


For  values  of  bcr  solutions  are  given  that  take  into  account  the 

influence  of  nonisothermalness  and  heterogeneity  of  the  boundary  layer. 

For  flows  with  significant  positive  pressure  gradients  it  appeared 

to  be  possible  theoretically  to  determine  the  critical  parameters  at 

the  point  of  breakaway  of  the  boundary  layer  from  an  impenetrable  wall. 

It  is  shown  that  breaking  away  is  characterized  by  a  critical  value 

of  the  form  parameter  f,  wher®  by  the  quantity  f„  is  a  weak  function 
## 

of  the  1ft  number.  Dependencies  of  the  critical  parameters  of  break¬ 
away  on  the  temperature  factor  and  the  M  number  are  obtained. 

In  connection  with  the  fact  that  the  critical  value  of  the  form 
parameter  f  remains  practically  constant,  but  not  the  form  parameter 
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Of 

r  =  - — ,  as  was  assumed  previously,  a  new  method  of  solution  of  the 
cf0 

momentum  equation  Is  proposed.  General  solutions  of  the  equations 
of  energy  and  momentum  for  flat  and  axi symmetric  boundary  layers  on 
a  impenetrable  wall  are  taking  into  account  the  nonisothermalness, 
compressibility,  and  pressure  gradient. 

For  permeable  walls  solutions  of  the  equations  of  energy  and 
momentum  for  relatively  small  values  of  the  form  parameter  f  are  given. 
The  relative  influence  of  nonisothermalness  and  mass  transfer  on 
friction  and  heat  transfer  in  the  turbulent  flow  of  gas  is  connected 
with  a  change  of  the  density  in  the  nucleus  of  the  boundary  layer  and 
minutely  depends  on  the  molar  viscosity  and  thermal  conductivity. 

Thus,  establishment  of  the  fact  of  the  existence  of  the  limiting 
relative  laws  of  friction  and  heat  transfer  for  a  turbulent  gas  stream 
allowed  us  to  give  a  certain,  logically  consecutive  method  of  research 
of  the  considered  problem. 
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APPENDIX  I 


SUMMARY  OF  BASIC  CALCULATION  FORMULAS 


I. I.  Streamlining  of  an  Impenetrable  Plate 
Coefficient  of  friction  during  Isothermal  flow* 

<2,5  l«***  +  3.SP  * 

Coefficient  of  heat  transfer  during  quasl-lsothermal  flow1 


(1.1) 


*.= 


2  Mr~  0,71 


(2.Sl«JSr;*  +  3f8)*  f1-2) 

Limiting  relative  laws  of  resistance  and  heat  transfer  (e  »  1 
assumed) : 

a)  during  subsonic  flow 


v-  t _ 1 _ 

(/T  +  iP 


b)  during  supersonic  flow 


(1.3) 


H*  = i f 

*•-»  l 


arc  sin 


2(*«-tn-Afr 


—  arc  sin 


»  4(**  “!)(*•  +  + 

_ _M  .  .. 


V  4(*»  -  !)(*•  4  (A*) 


xDuring  use  of  the  exponential  approximation# 

<f%  «=  B  Mi**~  *  51#  -  -j-  /»/•  aT5  iC'  "  . 

where  the  coefficients  B  and  m  are  taken  from  Table  3.1. 


(1.4) 


is 
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Relative  laws  of  resistance  and  heat  transfer  taking  into  account 
the  finiteness  of  the  |0  number  (e  »  1): 


a)  during  subsonic  flow 


[V  ♦  8.2(4  -l)  V  th  ♦  l]3 
b)  during  supersonic  flow 

^  -  p=rJ  X  [arc  .in - 

-  |>(l  _  8.2  I  | 


(1.5) 


-  arc  sin - igdfl‘-7»n 

/  4(+* -!)(*• 


4(4*-  !)(+*  +  A*)  +  <**.• 

i 


(1.6) 


The  relation  of  the  average  coefficients  of  friction  and  heat 
transfer  at  Twal^  -  const  and  a  turbulent  layer,  starting  from  the 
leading  edge  of  the  plate, 

(1.7) 

Average  coefficient  of  friction  for  an  isothermal  turbulent 
layer,  developed  from  the  leading  edge  of  the  plate 


er  y 

i  - 1 

f  r>  ) 

| 

*r.Q  1 

<  ">  J 

<>o  =  0,455  i  ig  Rex) 


-JW 


(1.8) 


Coefficient  of  heat  transfer  for  a  dynamic  tui-bulent  boundary 
layer,  developed  from  the  leading  edge  of  the  plate  and  the  initial 
heat-insulated  section  with  length  xQ 


•d_  °.0U9  ,  x  -  ,v.i  \0.w*t 


(1.9) 


I.  II.  Streamlining  of  a  Permeable  Plate 
during  Subsonic  Flow  (limiting  lawsl~ 


Relative  law  of  resistance 


Am 


(ii. i) 
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Critical  value  of  the  wall  permeability  parameter: 

a)  isothermal,  uniform  boundary  layer 


(II. 8) 


(H.3) 


b)  nonisothermal,  uniform  boundary  layer  at  ii  <  1 

c)  nonisothermal,  uniform  boundary  layer  at  V'  >  1 


(ii.*) 


-7 b-( 


•re  cos 


lz±. 


)• 


(II.5) 


a  graph  of  the  function  b  ( It)  is  given  Pig.  17; 

d)  isothermal,  nonuniforra  boundary  layer  at  R1  >  Rq  (m-  <  M-0)#  where 
H  is  the  molecular  weight  of  the  gas 


(II. 6) 


e)  isothermal,  nonuniform  boundary  layer  at  R’  <  Rq  ((i<  >  Hq) 

K,  «  1 .47  +  2,53  ;  ( II .  7 ) 

f)  nonisothermal,  nonuniform  boundary  layer 


4»*~0,25  bn  ,  (11,8) 

where  bcrl  and  bcr2  are  calculated,  respectively,  by  formulas  (II.4), 
(II. 5),  and  (II. 6). 

The  relative  law  of  heat  transfer  has  the  same  form  as  (II.l), 
but  instead  of  the  parameter  b  the  thermal  parameter  of  wall 
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permeability  is  introduced: 


r»i\  ' 

CptoWtSt., 


(II.9) 


Values  of  b, ,  __  =  b„_. 

t  cr  cr 

I,  III.  Critical  Parameters  of  Breakaway  of  the  Boundary 
Layer  from  an  Impenetrable  Wall 


Critical  parameters  at  point  of  breakaway  of  isothermal  boundary 
layer: 


-/„= 0,010; 

(lu.i) 

(III.2) 

(III. 3) 

See  also  Table  5.2. 

Critical  parameters  during  nonisothermal,  subsonic  flow,  see 
Figs.  56,  37 i  and  38. 

Critical  value  of  parameters  during  supersonic  flow,  see  Figs. 


39  and  40. 

Change  of  parameter  H  = 


/f-t+ojtr'(—v-!£_  V-;.  (111. 4) 

\  «.«-r  **  / 

where 

The  relative  change  of  the  coefficient  of  friction  under  the 
influence  of  a  pressure  gradient  is  determined  from  Fig.  28. 

##  C 

The  flt  number  In  the  range  of  values  Bt  <  10^  can  be  considered 
practically  independent  of  the  pressure  gradient. 
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I.  IV.  Integrals  of  the  Energy  Equation 


Subsonic  stream  on  an  Impenetrable  surface 


#*" = -jV  f  *,£4  r1  ♦  “  • 

+  (/?«r  4  r)i*  */*■. 

The  quantity  Is  calculated  by  formula  (1.3). 
Supersonic  flow  on  an  impenetrable  plate 


dX  + 


Re 7 


= -it 


i 

f/(l  -  t/^x 


xAr,  +  -,</*  + 


<«#r  4  r)i:  -r18. 


The  quantity  ^  calculated  by  formula  (1.4). 
During  supersonic  flow 


(IV.l) 


(IV.  2) 


A  7*=  7«t  -  rtr. 


(IV.  3) 


The  local  value  of  the  Stanton  number,  determined  from  the  calcu- 


lated  value  of 


During  subsonic  flow  m-qo  *  M-q* 
Subsonic  flow  on  a  permeable  surface 


**r=-A-! 


x 


-f 


(IV.4) 


»,] A  T'*  mdX  -r  [Re? A  rji*  "}  (IV. 5) 


Distribution  of  supply  of  cooling  gas  along  the  length  of  the 
contour  (coordinate  X) 

'*  (IV.6) 
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Local  value  of  the  Stanton  number 

»  =  «',(  I  -tM'®-  (IV.  7) 

Local  heat  flux 


Vci  ~  S  f*0  ^  T  St . 


(IV.8) 


Initial  temperature  of  the  cooling  (injected  through  the  wall 
into  the  boundary  layer)  gas 


Tt  .*=  rtT  - 


fct 


/I 


During  evaporation  and  sublimation, 

9tt  =  rg/t; 

St 


St, 

K~ 


=  *,  K; 

* 

tpi*T 


')'•  ■ 

**'  =  TT  {  IrA  flR<0|  [*"■’  iff1 

x(/ 


(IV.9) 

(IV.  10) 
(IV.  11) 
(IV.12) 

(IV.  13) 


+(#*rA7)\; 


I  + 


(IV.  14) 


The  quantity  is  calculated  by  formula  (I 
In  the  case  of  gas  injection  and  given 
Tq  in  formulas  (IV.  10)  -  (IV.  15)  instead  of 
duced  the  ration 


5). 

values  of  T±,  Twall,  and 
the  criterion  K  is  intro- 


r„  -  h 

t„  -  r„ 
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I, V.  Integrals  of  the  Momentum  Equation 


Subsonic  flow  on  an  Impenetrable  surface 

x 


He"  =  5f[ — BHe.  f  «’  ' " V ,dX  -r 

1  X. 

+  (He"  i-’jV*  *f"rS i 


(V.l) 


where 


*  =  1  -r  Hn. 

The  quantity  Is  determined  by  formula  (1.3).  The  quantity 

Hcr  is  determined  from  the  graph  of  Pig.  37. 

Local  value  of  the  coefficient  of  friction,  determined  from  the 

#  * 

calculated  value  of  !•  . 


c,  =  *;.«!•,., 


/*• 


(V.2) 


where  is  determined  from  Fig.  28. 


Value  of  the  form  parameter 

/= 


l 


(V-5) 

Critical  value  form  parameter  is  determined  from  the  graph  of 
Fig.  36. 

Supersonic  flow  on  an  impenetrable  plate 

x  1 

/?#••  ==  IT 1  |  (1  -  L'*jTrT>; 

"".dx-r  (Re* * ■}  +  " .  ( V. 4 ) 

The  quantity  is  determined  by  formula  (1.4).  The  quantity 

H  is  determined  from  the  graph  of  Fig.  40. 
c  r 

The  value  of  the  form  parameter  for  the  calculated  value  of  It*1 


»« 


/  = 


Me»U*{l  -  V) 


1 

1-  » 


4U 

4X 


(V.5) 
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The  critical  value  of  the  form  parameter  is  determined  from  the 
graph  of  Fig.  39. 

The  local  value  of  the  coefficient  of  friction  for  R»  ,  deter¬ 
mined  by  (V.4), 

(v-6) 

Subsonic  flow  on  a  permeable  surface  at  small  f 

BRe.  ^  + (I  -  -—■)'  + 

X\ 

+  t\dX  4-  iRe*^yx;  (V*7) 

## 

Local  value  of  the  coefficient  of  friction  for  Re  by  (V.7) 


_L_)V  (v>8) 

I. VI.  Integrals  of  the  Equations  of  Energy 
and  Momentum  for  •Axisymmetrlc  Flow 

During  axisymmetrlc  flow  the  integrals  of  the  energy  equation 
have  the  same  form  as  in  the  case  of  flat  flow  (see  Section  I. IV), 
but  Instead  of  the  quantity  AT  the  product  RAT  enters,  where  R  is  the 
radius  of  the  surface  of  the  body. 

Correspondingly,  in  the  integrals  of  the  momentum  equation  during 
subsonic  flow  instead  of  the  quantiy  Wq,  the  quantity  R  w^  is  intro¬ 
duced. 

Local  coefficient  of  heat  transfer  in  an  axisymmetrlc  nozzle 


(VI.l) 
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where 


*#7 


r  <1 +«)**#.., 

f  D>p  mt 

A-l  j 

[  2Pt'n  ST**  m 

l  D  )  '  1 

1  *  +  l  / 

(VI. 2) 
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APPENDIX  II 

EXAMPLES  OF  CALCULATIONS 
A.  Heat  Transfer  In  a  Supersonic  Nozzle 


There  exists  a  supersonic  nozzle,  the  geometric  dimensions  of 
the  flow-through  part  of  which  are  shown  in  Fig.  4l  and  in  the  first 
two  lines  of  Table  II.  1.  The  parameters  of  drag  of  air  in  the 
nozzle:  pQ0  «  4.49;  T*  =  5.99.6°K. 

Temperature  of  the  wall  T  11  =  547°K  is  constant  along  the 
length  of  the  nozzle. 

Distribution  of  the  local  coefficients  of  heat  transfer  along 
the  length  of  the  nozzle  and  the  corresponding  distribution  of  heat 
fluxes  on  the  wall  is  found: 

q-a(r„  -  T„). 

Parameters  of  air  T  =  T^:  Fr00  =  0.725; 

|IW«3,I2  .  10_*  fcg-Me/m2!  C,  —0,251  koal/kg'degl 


4,49  •  10« 

9.81  •  29,77  •  509,6 


=  2,60 


Hence 


-*0,0381  ko&l/n'hr'deg. 


Re,»,  -  -jy~_r  V  427 . 2  •  9.81  •  0,251  •  599,6  =  1,9  •  10'. 

Take  m  =  0,25  and  B  ■  0.0129.  Results  of  further  calculations 
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are  given  in  Table  II. 1. 


B.  Cooling  of  a  Porous  Plate 


A  plate  of  length  i  m  is  streamlined  by  a  flow  of  air  with  a 
temperature  tQ  *  2000°C,  and  is  cooled  by  cold  air,  passed  through 
the  pores  of  the  plate,  flow  velocity  of  the  hot  air  wQ  ■  50  n/sec. 

It  is  necessary  to  determine  the  specific  flow  rates  of  the  coolant  air, 
needed  to  maintain  the  temperature  of  the  plate  constant,  equal  to 
‘wall  *  500°C. 

The  initial  temperature  of  the  coolant  air  is  t^  »  50°C  (diagram 
of  problem  is  given  in  Pig.  II. i). 

By  formula  (II. 4)  we  determine  the  critical  parameter  of  injec¬ 
tion 


By  formula  (IV.  15)  we  determine  the  quantity  b^: 

MO -30  7,7W» 


*,  =  7,793- 


1.1 


MCM00-5Q0) 


-») 


+  1 


-)= 


X 

2,44. 


By  formula  (5.94)  we  determine  the  parameter  f: 


V  =  2.44  • 


300—30 

*40-500 


=  0,764. 


Fig.  II. 1.  Diagram  for  calcu¬ 
lation  of  the  cooling  of  a 
porous  plate. 
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Tab±e  11,1.  Combined  Table  of  Calculation  of  the  Thermal 


itt.  •*  tmmi 


By  formula  (V.7)  we  determine 


«» 


»« =  -afi ft-  - 

fc  =  0,0I5S  — — —  ; 

=  (1.1  •  10*  ^  0,0128 . 3,204  •  l.25)M  1000 X'*. 


By  formula  (1.1)  we  determine  the  local  values  of  cf0.  After 
substitution  of  the  results  of  the  preceding  calculations  we  have 


•*  _  0,002  6 
2  X™  ‘ 


We  determine  the  flow  rate  of  coolant  airs 


To*«t  —  V* 


c/t  7,5  •  2.44  •  0.0022b 


X*' 


We  calculate  the  local  values  of  the  coefficients  of  friction. 


cf 


0.00342 

A'1’* 


Table  II. 2.  Results  of  Calculation  of  Cooling  of  a  Porous  Plate 


X 

M 

0J 

U.4 

M 

(•A 

M 

w 

M 

u> 

y  mllwt»ll 

[m2«o.  J 

0.0685 

0,0661 

0.0625 

0,0403 

0,0472 

0,0450 

0,0443 

0,0482 

0,0421 

0.0413 

e/*10» 

5,48 

4,78 

4,3* 

4.13 

3,96 

3,62 

3.7 

3.62 

3.12 

3,44 

W10*4 

pssall 
Jjtfhr  J 

2.66 

2,3 

0.13 

2.0 

1,92 

MS 

1.0 

1,75 

1.71 

1,076 

We  determine  the  specific  heat  flux: 

fct  =  C,  •  T*tW*r  (At  “A)  = 

0.34  •  7.5  •  2.44  •  0,00226  (500  -  30) 

= - P3 - 

Results  of  the  calculations  are  reduced  in  Table  II. 2  and  are 
presented  in  Fig.  II. 2. 
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Pig.  II. 2.  Results  of  calculation 
of  cooling  of  a  porous  plate. 


C.  Heat  Transfer  in  the  Inlet  Section  of  the  Pipe 


In  a  cylindrical  pipe  of  D  =  0.011  m  air  enters,  with  parameters 
at  the  inlet,  TQ  =  297. 5°K;  G  =  22.7  kg/hour,  u  =  1.85.10-6.  The 
temperature  of  the  wall  of  the  pipe  Is  constant  and  equal  to  T  ^  = 

=  572. 5°K.  It  is  required  to  determine  the  distribution  of  the 
specific  heat  flux  along  the  length  of  the  pipe. 


We  find  that 


^n.  — 


4  *  22,7  _  .  ^  jfu. 

~  3800-  9,81  •  1,83  •  10"*  •  3,14  •  0,11  ~  ' 

*>•  =  0,725;  PeD  = 

•=  Pr  =  2,95  .  10* ; 

•{*—1,25. 

Equation  (6.49),  taking  into  account  the  given  parameters,  is 


reduced  to  the  form: 


==  60 jo, €16  •  892,318  X 

_£»”  +  II.W  /■«“”  -r  67.4  f  2J(C 

Pt"°*  -  11.6 8P#  J  6/,4  67,4  -  JJ 


Assigning  the  values  Ft  =  200  ,  400,  1000  ,  2000  ,  4000,  by  this 
equation  we  determine  the  corresponding  values  of  X. 


By  formula  (6.50)  we  determine  the  specific  heat  flux.  Taking 
into  account  the  given  parameters,  equation  (6.50)  takes  the  form: 


fw  —  193 


(2.95  10* +  6,5 /»#"») 


Results  of  the  calculations  are  reduced  in  Table  II. 5.  In  Fig. 
II. 3  are  presented  the  results  of  calculations  and  experimental  data 
of  B.  S.  Petukhov,  obtained  for  the  indicated  parameters.  As  can  be 
seen  from  the  graph,  the  agreement  of  the  analytic  calculation  with 
the  experiment  on  heat  flows  as  well  as  by  values  of  are  fully 

satisfactory.  On  the  same  graph  are  presented  the  results  of  calcu¬ 
lations  for  streamlining  of  a  flat  plate  by  the  formula  of  B.  S. 
Petukhov.  From  the  graph  it  is  clear  that  the  best  agreement  with 
experiment  Is  given  by  the  proposed  method  of  calculation. 


Table  II. 3.  Results  of  calculation  of 
Heat  Transfer  in  the  Inlet  Section  of 
the  Pipe 


#»• 

D 

*«»  ' 

Hoi 

JMtt 

j  «m 

X 

1.49 

3,74 

10,45 

22.4 

47.6 

W1,H ; jr£r 
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12,52 
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13.7 
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Fig.  II. 3.  Results  of  calculation  of 
heat  transfer  in  the  inlet  section  of 
the  pipes  1  —  qwan  according  to  the 

method  of  the  authors  2  —  qwan  accord¬ 
ing  to  the  method  of  B.  S.  Petukhov,  3  — 
qwan  by  the  empirical  criterial  formula 

for  stabilized  heat  transfer  in  a  pipe, 

4  —  Ft  according  to  the  method  of  the 
authors,  O  —  experimental  data  of  B.  S. 
Petukhov  for  ®  —  experimental 

data  of  B.  S.  Petukhov  for  . 
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